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CHAPTER  1 

INTRODUCTION  AND  THE  STANDARD  MODEL 
11  The  Rjndampntftl  Particle.g 

In  the  standard  modal  all  mattar  is  composed  of  electrically  charged  fun- 
denicntal  particles  known  as  quarks  nod  leptons.  Ordinary  mattor,  Ibai  Is 
atoms  and  molecules,  consists  of  two  types  of  quarks,  the  up  and  down  quarks 
arranged  into  protons  and  neutrons  (nucleons),  and  one  type  of  lepton,  the 
electron-  The  uudeons  are  confined  in  the  niicleu.s  while  electrons  arc  found  in 
quantum  mechanical  orbits  in  and  around  the  center  of  the  atom  or  molecule. 
The  quarks  differ  from  leptons  in  two  fundamental  ways:  they  have  fractional 
electric  charge  and  possess,  in  addition,  a distinct  quantum  number  analogous 
to  the  electric  charge  known  as  color.  The  color  charges  come  in  three  different 
varieties  known  commonly  as  red  (R),  green  (G)  and  blue  (B)  and  give  rise 
to  the  strong  interactiotis.  Since  quarks  are  endowed  with  both  electric  nod 
color  cliargcs  they,  unhke  leptons,  are  infiuenced  by  both  electromagnetic  and 
strong  forces. 

The  up  and  down  quarks  and  the  electron  with  its  doublet  partnet,  the 
electron  neutrino,  make  up  the  first  fermionic  generation  from  which  all  of  the 
familiar  matter  is  constructed.  We  now  know  that  two  more  generations  of 
doublet  pairs  exist;  they  duplicate  almost  exactly  the  first  generation  but  pos- 
sess significantly  larger  masses.  These  two  new  geuerations  consist  of  doublets 
with  progressively  heavier  masses  culminating  with  the  recently  discovered  top 


Figure  1,1  The  fundanieitliil  particles. 


quark  with  a mass  of  176  GeV  (roughly  the  mass  of  a timgsteu  atom)  (1],  Ad- 
diliouai  generations  witli  neutrino  masses  leas  then  45  GeV  have  been  ruled 
out  by  LEP  meassureinent  on  the  width  of  the  Zfi.  The  origin  of  the  heavier 
geuerationa  remains  u mystery. 

Ill  addition  to  the  fermions  (spin-half  particles]  a set  of  fundamental  inte- 
ger spin  particles  known  as  vector  bosons  are  known  to  exist.  These  particles 


mediate  the  weak,  eicctromagmetic.  strong  and  gravitaliooal  forces,  f In  quan- 
tum field  theories,  particle  interactions  are  n result  of  the  exchange  of  gauge 
bosons.  For  example,  the  photon  is  the  particle  exchanged  in  eleolromagnetic 
interactions.  Its  exchange  mediates  electromagnetic  interactions  and  thus  all 
chemical  and  light  phenomena  can,  in  principle,  be  described  as  a consequence 
of  the  exchange  of  the  photon.  Eight  gauge  bosons,  known  as  gluons,  medi- 
ate the  strong  interactions  which  hind  quarks  into  hndroos  and  nucleons  into 
nuclei.  Both  gluons  and  photons  ate  massless  and  are  the  force  carriers  or 
"gauge  quanta”  of  the  theories  describing  the  dynamics  of  the  respective  in- 
teractions. In  weak  interactions  three  spin-one  fundamental  bosons,  two  with 
opposite  electric  charge,  and  one  with  no  charge  are  exchanged.  Unlike  the 
fcuce  carriers  of  the  strong  and  electromagnetic  interactions  the  weak  force 
carriers  arc  massive.  The  weak  interaction  allows  quark  and  lepton  species  to 
transform  amongst  tbcmseives  and  is  responsible  for  radioactive  beta  decay  in 
heavy  nuclei.  The  fermions  and  vector  bosons  arc  listed  in  Figure  1.1  with  their 
approximate  masses,  intrinsic  charge  electric  and  number  of  color  charges. 


Particles  with  quantum  numbers  consistent  with  those  associated  with  free 
quarks  have  never  been  observed  in  isolation.  This  has  led  physicists  to  postu- 
late that  the  strong  interactions  exhibit  a property  known  as  color  confinement. 
The  hypothesis  asserts  that  quarks  are  so  strongly  bound  within  hadrons  that 
they  can  not  escape  as  isolated  entities  but  must  lombine  with  other  quarks 

1 A spin  two  boson  is  postulated  ns  the  carrier  of  the  gravitational  force. 
Gravitational  interactions  between  subatomic  patiicles,  however,  are  negligible 
when  compared  to  the  other  three  forces  and  are  thus  excluded  from  further 
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into  final  stales  tehlcli  posses  no  net  color.  In  tbe  language  of  group  theory  a 
state  which  possess  no  net  charge  is  known  as  a coIor-singlet. 

As  already  stated  quarks  come  in  three  varieties,  red.  green  and  blue. 
These  compose  a fundamental  triplet  representation  under  tbe  color  SU(3)c 
symmetry.  To  form  hadrons  we  look  for  quark  combioations  which  traDsforui 
as  singlets  under  5£/{3)c  . The  simpiest  of  these  are  a quark  aDti>quark  pair 
aad  the  three  quark  combinatioo  of  ted,  greeo  and  blue  quarks.  These  parti- 
cle combinatbus  cmrespond  to  the  following  direct  product  decompositioo  of 
irreducible  5tf{3)  repccsentatioos 

3®3-Sei  aod  383®3  = 108S®8ei.  (1.1) 

The  quark  anti-quark  pair  forms  ao  octet  and  a singlet  and  is  known  as  the  me- 
son nonet.  The  states  in  the  three  quark  decomposition,  properly  symmelrired, 
are  known  as  baryons. 

Of  particular  ioterest  to  this  analysis  are  the  “heavy  flavor"  mesons  which 
consist  of  the  pairing  of  either  a c or  fc  quark  with  an  up  or  down  light  quark. 
Heavy  flavor  nicsoos  will  be  treated  separately  b the  folbwing  chapters  with 
emphasis  on  the  B meson  aad  its  decay.  Tbe  light  mesons  and  heavy  charmed 
mesons  are  shown  in  Figure  1.2.  The  figure  shows  the  mesoits  arranged  b 
the  badly  broken  SU(i)j  flavor  symmetry  of  the  quark  model.  (Exact  5f/(4) 
would  give  an  uumlxed  15-plet  and  a sioglet,  but  the  mass  of  the  c quark 
breaks  the  flavor  symmetry  and  allows  the  muting  between  the  5t/(3) / triplet 
and  the  5tf(4)y  singlet.)  A similar  analysis  of  the  synainetry  in  flavor  St^(3)y 
gave  physicists  suggestive  evidence  of  a deeper  structure  in  elementary  particles 
and  led  eventually  to  the  creation  of  the  modeni  quark  model. 


Figure  1,2  The  SU(i)f  flavor  17?  multiplels. 

In  (a)  the  16-amItiplel  for  the  vector  spin-one  mesons  in  (b)  the  spin  seen 
pseudosralara.  Tlie  middle  level  corresponds  to  the  S£f(3) / flavor  octet.  Tlie 
figure  is  taken  from  Ref.  (2), 


f' 


Figure  1.3  The  5C'(4)^  flavor  three  quark  multiplets. 

In  (a)  the  20-muitiplet  for  the  spin-3/2.  In  (b)  the  spin  1/2  ground  state 
baryons.  Tlie  ground  floors  in  (a)  and  (b)  correspond  to  the  SU{3)j  flavor 
decuplet  and  octet  respectively.  In  (b)  the  Brat  floor  shows  two  members  in 
the  I3  = position;  the  and  the  These  correspond  to  the  two 
possible  states  witli  the  same  third  component  oI  I'spin  Ij  but  different  total 
V'.  The  figure  is  taken  from  Ref.  [2). 


As  described  previously,  tlie  baryons  are  composed  of  color-siuglet  prod- 
ucts of  three  quark  iriplels.  UsUig  the  fiavor  5£/(4) / symmetry  we  tiod  that 
the  64  possible  qqq  combiuBtiaas  are  decomposed  into  a fully  symmetric  (S) 
2Q-plet  with  a Sl'(3)  decuplet,  Iw-o  other  20-plets  with  S£/(3)  octets  of  mixed 
symmetry  and  a 5tf(4)y  quadruplet.  Tlie  5f/(3)  octets  of  mixed  symmetry 
are  either  symmetric  (.i/5)  or  anti-symmetric  (A/^)  under  the  interchange  of 
the  Brst  two  quarks.  In  Figure  1.3(a)  the  particles  associated  with  the  totally 
symmetric  combination  of  the  product  of  the  symmetric  SU{i]j  23-plet  and 
symmetric  spin  SU(2)  multiplets  give  the  ground  state  spm  3/2  baryons.  To 
arrive  at  the  spin  1/2  ground  state  baryons.  Figure  1.3(b).  we  form  totally  sym- 
metric condtinations  from  tbe  two  20-plels  and  two  spin  multiplets  of  mixed 
symmetry.  The  symmetric  combination  is  necessary  to  form  a totally  anti- 
symmetric wavefimction  states  when  including  the  completely  antisymmetric 
color-singlet  part  into  the  total  wavefunction.  For  the  spin  3/2  baryous  the 
total  wavefunction  = (205,45,1^),  expressed  in  terms  of  the  spin,  flavor 
and  color  multiplicities  SU(4) j,  SU{2),.  SC/(3)e  is 

<'■  = 3 [(20Afs.4AJs,  U)  + (20a,..,  4a,...  U)]  - 

The  flavor  and  spin  decompositions  are 

S£'(4) / :4  a 4 ® 4 = 20s  ® 8a,j  ® 8.U.,  ® 3a 

SU  (2),  :2  ® 2 8 2 = 4s  ® 2a,5  ffl  2 a,.  . 

for  color  the  decomposition  is  given  in  by  Equation  (1.1). 

Because  of  the  large  mass  of  the  c quark  relative  to  the  Ughi  quarks,  flavor 
SV{i)f  is  not  a useful  symmetry.  However,  it  provides  a eonveoient  framework 
for  presenting  all  charmed  mesons  and  baryons.  In  fact,  the  quarts  are  iisualb' 


divided  into  the  light  quarki>,  the  u.d  and  s.  and  lieavj'  quarks,  c,  i and  t.  The 
approximate  flavor  symiiiclry  of  the  light  quarks  is  used  axlensively  in  chiral 
perturbation  theory  which  attempts  a quantitative  description  of  lots-  energy 
QCD  plieuomena.  At  the  other  end  of  the  nia.ss  spK-tra  the  large  masses  of 
the  c and  b quarks  are  exploited  in  Heavy  Quark  Elfertivc  Theory  (HQET)  to 
achieve  a model  independent  description  of  the  behavior  of  hadrons  containing 
heavy  flavored  quarks- 

1.3  St.'f31  X SI'I2)  X Uil)  : The  Standard  Modol 

The  description  of  the  snbatomic  world  given  in  the  previous  sections  con- 
sisted maiuiy  of  a discussiou  of  the  composition  of  hadrons  without  much  em- 
phasis placed  on  the  dynamics  of  the  constitueuis.  In  the  remainder  of  this 
chapter  we  focus  our  attentioo  on  a somewhat  brief  but  general  description  of 
the  fundamental  interactions  which  influence  the  microscopic  behavior  of  the 
fundamental  fermions. 

Of  primary  importance  in  describing  the  dynamics  of  the  fully  interacting 
particle  theories  is  ihe  principle  of  gauge  or  phase  invariance.  This  principle 
asserts  that  the  equatioos  which  govern  fundamental  particles  do  not  chaoge 
when  an  arbitrary  phase  is  introduced.  If  the  phase  is  a constant  then  by 
Noether’s  theorem  a conserved  quantity  results.  If  the  constant  pha.se  is  asso- 
ciated with  tile  charge  of  the  fuudaniental  particle  then  global  gauge  invariance 
implies  the  conservation  of  the  charge.  By  generalizing  the  constant  phase  to  a 
local  or  space-time  dependent  phase  and  imposing  invariance  under  the  gauge 
transformation  the  dynamics  of  the  theory  can,  in  principle,  be  derived-  Theo- 
ries which  are  gauge  invariant,  locally,  are  known  as  gauge  theories  aud  are  be- 


lievw)  10  allow  a description  of  all  the  microscopic  behavior,  that  U the  strong, 
weak  and  electromagnetic  interactions  between  the  fundamental  fermions. 

1.3.1  Quantum  Electrndvntiinir.s 

Tb  best  illustrate  the  generalisation  of  the  global  gauge  symmetry  to  a 
local  one  we  will  do  so  by  Erst  applying  the  procedure  to  the  prototypical  gauge 
theory  of  particle  imeraetioDs;  the  theory  of  quantum  electrodynamics  (QED). 
In  QED  the  dynamics  of  the  theory  can  be  completely  derived  by  demanding 
invariance  of  the  free  particle  Lagraiigian  under  a local  gauge  transformation. 
We  begin  with  the  Lagrangiau  for  a free  fermion  field 

Co  = (1.2) 

In  a local  gauge  transfonnatioii  the  constant  phase  is  generalised  by  making  it 
a space<timc  dependent  function 

= (1.3) 

This  transformation  no  longer  preserves  the  form  of  the  original  free  particle 
Lagrangian  since  the  derivative  acts  on  the  space-time  dependent  phase 

a„i/t  -*  + fc'^l'liitaao  (1,4) 

and  produces  u term  proportional  to  d„a  which  breaks  the  invariance  of  Co. 

lb  restore  covariance  to  the  Lagrangian  we  replace  the  derivative  with  the 
covariant  derivative 

This  covarinnt  derivative  transforms  like  jh{x)  i.e..  D^il>{x)  -t  e’“**)D„di(i). 
The  transformation  law  of  the  covariant  derivative  determines  the  transfor- 
mation property  of  the  vector  field  which  will  introduce  a term  that  esactly 


cancels  the  unwonted  O^q  in  Equation  (1.4).  The  \'ector  field  thus  transforms 


Ap  -»  Ap  =Ap  + -8p 

and  thus  gives  rise  to  aji  interartion  term  in  the  original  free  particle  La- 
grangian.  This  field  has  the  form  expected  for  the  electromagnetic  potential 
and  is  associated  with  the  spin-one  photon  field. 

The  Classical  QED  Lagiangian  can  now  be  written  down  from  the  free 
particle  Lagranglan  with  the  replacement  of  the  derivative  with  the  covariaat 
derivative.  Since  we  want  to  associate  the  gauge  field  A„  with  a physical  photon 
field,  we  need  to  add  a covariant  kinetic  energy  term  for  the  vector  field  to  the 
origiaal  Lagranglan 


arriving,  finally,  at  the  full  QED  Lagranglan 

Cqed  = t5(>'9  - e4)4-  - m^ti.  - (1.5) 

The  tree  level  Feynman  rules  can  he  read  directly  off  equation  Equation 
(1.5).  The  coupling  of  fermion  fields  to  the  photon  vector  field  is  of  strength 
c and  is  given  by  the  4i44i  term.  We  alTO  note  that  the  a mass  term  (ApA’‘) 
for  the  vector  field  is  absent  from  the  Lagranglan  which  implies  that  the  gauge 


ineulum  transfer  of  the  exchange  photon.  From  the  dijfjj  term  the  propagator 
for  fermion  fields  is  proportional  to  !/(;)  -m). 


Id  the  qiuuitizaclon  procedure,  using  the  path  integral  formuiatioii,  the  clas- 
sical Lagrangian  Equation  (1.5)  is  laed  as  au  input  to  the  effective  Lagrangian 


The  higher  order  terms  generate  loop  corrections  to  the  classical  Lagrangian 
which  introduce  new  interactions  and  or  finite  rorrections  to  the  basic  inter- 
actions. Evaluating  the  loop  cootrihutions  Introduces  divergences  which  can 
be  effectively  removed  by  renormalising  the  input  parameters  and  fields  of  the 
theory.  The  reucrmaliaation  procedure  results  in  the  modification  of  the  pa- 
rameters which  make  them  dependent  on  an  arbitrary  scale.  The  parameter 
can  he  evaluated  at  convenient  scale  and  then  extrapolated,  through  the  use 
of  the  renonnaliiation  group  equations,  to  the  scale  of  interest.  Fbr  example, 
the  renormalization  group  equations  give  the  running  of  the  QED  coupling 
constant  to  all  orders  in  perturbation  theory  as 


-=S?i™(S)  ■ 


(1.6) 


Here  « is  the  arbitrary  scale  and  17  is  the  mass  scale  of  the  interaction.  A 
convenient  point  for  QED  calculation  is  usually  taken  at  w 2m|  where 
is  the  well  known  value  of  1/137. 


1.3.2  Quantum  Chromodvnamics 


While  few  similarities  appear  to  exist  between  strong  and  eiectro-niaguetic 
interactions  the  same  gauge  principle  used  to  derive  the  dynamics  of  QED  can, 
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whpD  Qppli(>d  10  the  SV  (3)  color  symmetry,  be  used  to  derive  the  dynamics  of 
the  strong  interactiou  (3,4).  We  begin  as  before  with  the  free  quark  Logrengian 

do  = ?W(>^'-m)q(r)  (1.7) 

in  analogy  with  QED  Equation  (13)  a gauge  trausformatioD  takes  the  form  of 
?(r)  -» 

where  Oq  is  an  space-time  dependent  phase,  the  X°  are  the  Gell-Mann  matrices 
that  satisfy  ibe  SC/(3)  commutation  reiations  and  normalization  conditions 

].«■*■£  , »(«•].  2 J“.  (1.1) 

and  the  unitary  transformation  f/(i)  is  now  an  arbitrary  3x3  unitary  ma- 
trix. Unlike  the  U(l)  case,  the  generators  of  SU(3)  do  not  commute.  This  has 
important  consequences  to  the  possible  interactions  of  the  theory  since  it  im- 
pbes  that  the  gauge  bosons  themseives  carry  charge  and  can  thus  self  interact. 
Gauge  groups  with  non-commuting  generators  are  know  as  noo-Abeliaa  gauge 

We  now  proceed  as  in  QED  and  demand  covariance  by  replacing  tbe  deriva- 
tive 8„  with  tbe  covariant  derivative 

= ; n = l-  .8. 

where  T,  = ^ and  g,  is  the  strong  coupling  constant.  The  transformation 
of  the  vector  field  is  augmented  by  an  ndditional  term  to  cancel  the  extra 
terms  introduced  by  the  gauge  transformations  of  a non-AbeUan  group.  The 
transformation  of  the  gluon  field  is  thus  given  by 

-*  GJ  - 


(1.9) 


To  complete  the  QCD  LaBrangion  we  need  to  add  a kiueiic  energy  term  which 
is  mvarionl  under  the  transformation  similar  to  Eqiiatioo  (1.9).  The  covariant 
kinetic  energy  term  is 

•FJv  = SnC;  - 3vG2  - 9/,6cG?,CJ 

Adding  the  kinetic  energy  term  to  the  original  Lagranginn.  Equation  (1.7), 
and  replacing  the  derivative  with  the  covariant  derivative,  we  arrive  at  the 
classical  QCD  Lagranginn 

AjCfl  = + 9(1)  (n“8n  + 9(1)  - 9(i)m9(i)  (1.10) 

where  the  quark  fields  are  described  by  9(1)  and  the  gluons  by  the  gauge  fields 

03. 

Similarly  to  QED  v,«  can  qualitatively  derive  the  Feynman  niles  by  in- 
spection. Propagators  for  the  quarks  and  gluons  are  obtained  from  the  terms 
proportional  to  99  and  A®  respectively  aod  quark  gluon  couplings  of  strength 
S are  given  by  999A  term.  In  addition,  tho  aon-Abelian  fields  in  the  kinetic 
energi’  term  introduce  additional  interactions  that  involve  gluoms  coupling  to 
themselves.  These  terms.  g,A^  and  g^A*.  give  rise  to  three  and  four  gluon 
CDupIinp  of  strength  g,  and  s|  respectively  and  have  no  analog  in  QED. 

The  self-interactions  of  the  gluons  are  a result  of  the  non- Abelian  character 
of  the  gauge  group  and  imply  that  the  gauge  fields  themselves  carry  the  color 
charge.  Observationally,  this  property  of  the  theory  manifests  itself  in  the 
confinement  of  the  quarks  in  hadronic  final  slates.  The  theory  alsoexhibiM  the 
property  known  as  asymptotic  freedom  which  states  that  at  short  disiamres  the 


force  between  quarks  goes  to  zero-  This  is  illustrated  in  the  cunning  coupling 
constant 

(H-|n,)lnya/A^^j, 

where  a,  is  expressed  as  a function  of  the  energy  scale  17  and  Ajjcjj  is  the  scale 
at  which  the  coupling  constant  becomes  infiuite.  Asymptotic  freedom  follows 
since  as  Oj  -*  0,  g -»  00.  In  Equation  (1.11).  the  group  properties  of  5(/(3)c 
are  reSected  in  the  value  of  the  constant  11- 

I-3.3  Ute-Theorv  of  Electroweak  Imemciions 

The  weak  interacibn  is  not  only  responsible  for  the  transmutation  of  quarks 
and  leptons  species  but  also  changes  the  charge  of  the  quark  or  lepton  undergo 
the  transition.  The  fact  that  both  flavor  and  charge  undergoing  a change  during 
the  same  transition  suggests  that  the  weak  and  electromagnetic  interactions 
are  different  aspects  of  a single  theory.  We  now  believe  that  the  weak  and 
electromagnetir  forces  are  indeed  unified  aud  are  described  by  a gauge  theory 
commonly  referred  to  as  the  standard  model  of  the  electroweak  interactions 

[5]. 

The  masses  of  the  Intermediate  vector  bosons  mediating  the  weak  inter- 
actions have  been  measured  and  are  80.2  GeV  for  the  H'*  and  91.2  GeV  for 
the  Z®  [2]-  The  masses,  however,  cannot  be  added  to  the  Lagrangian  in  terms 
bilinear  in  the  fields  since  this  would  break  gauge  invariance  and  thus  render 
the  theory  unienormalizable.  In  the  Weinherg-Sabtn  model  the  gauge  bosons 
are  initially  massless.  They  acquire  mns-s  wheu  the  vacuum,  which  does  dot 
posses  the  original  symmetry  of  the  Lagrangian.  is  selected  from  among  the 
many  possible  vacuum  configurations.  In  this  way  masses  for  the  vector  besons 


ate  generated  in  a gauge  invariant  way  rendering  the  theory  reaotmalijable  [6]- 
In  the  electroweak  theory'  three  of  the  four  massless  vector  bosons  acquire  amec 
while  one  remains  massless.  The  details  of  how  this  is  accomplished  through 
the  Higgs  mechajiism  will  be  described  in  the  following  sections. 

h’nbrolten  elecirowcak  theory  The  electroweak  theory  is  based  on  the  non- 
Abelian  gauge  group  formed  by  tbe  product  of  a left-handed  S£/(2)i  and  tf(l)y- 
for  hyperebarge  or  S£/(2)iX  {.'(I)!-.  The  form  of  this  gauge  group  is  motivated 
by  the  need  for  au  interaction  which  violates  parity,  is  invariant  under  rotations 
in  weak  isospin  S£/(2)x  and  hypercharge  £/(l)y  space,  and  gives  at  least  four 
gauge  bosons:  two  massive  charged  vector  bosons  a massive  neutral 

boson  (Z®)  and  a massless  vector  boson  (photon).  The  simplest  group  structure 
that  esa  accommodate  these  conditions  is  the  product  of  an  Sb'(2)x  x l/(l)y. 
Before  spontaneous  symmetry  breaking  all  gauge  bosons  are  massless  and  the 
arguments  that  led  to  the  detivalbn  of  the  dynamics  of  QCD  apply.  The 
unbroken  Lagrangian,  in  analogy  with  QCD  is 

Ci=-^lFl^r‘^  + G^C'^)  (1.12) 

Fi^  = d^Wl  - S„w;,  + ge‘’‘‘wiw^.  i = 1,2,3 


The  potential  energy  term,  also  know  as  the  Weyl  Dirac  term,'  is  given  by 


' At  this  point  m the  detivalioo,  the  Weyl  Dirac  part  of  tlie  Lagrangian 
dots  not  have  a term  bilinear  in  the  fermion  fields  since  this  would  break  the 
31/(3)  X 1/(2)  symmetry.  Fermion  mass  terms  will  be  introduced  in  the  Higgs 
sector  and  are  also  a product  of  spontaneous  symmetry  breaking. 


Table  of  feruiiou  weak  isospio  and  bypercharge  assignments 
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Demand  of  local  gauge  invariance  gives  the  form  of  tbe  covariant  derivative 


Here  T — Tt  and  the  r are  tbe  Pauli  matrices  that  satisfy  the  eommutation 


and  generate  the  St/(2)  group.  The  £/(I)  group  is  generated  by  the  weak 
hypereharge  V = 2(Q  - 73). 

The  quantum  number  assignment  for  all  fermions  is  shown  in  Tkble  1.1. 
The  group  structure  allows  for  any  choice  of  hypercharge  assignment  but  the 
particular  values  chosen  allow  for  the  cnucelatiou  of  Adler-Bell-Jackiw  (ABJ) 
anomalies . With  the  quantum  number  assignments  above  tbe  covariant  deriva- 
tive acts  oQ  tbe  left-handed  lepton  fields  ns 


(ffi.Oj]  = 2«,j,.ot 


(1.13) 


and  righl'banded  skiglei  hk 


D„f„  = (a„  + ij'S,,)e,. 

SpoDianfom  i\  mmtifv  broakiin.  To  gsnernte  masses  for  three  of  the  gauge 
bosons  while  preserving  the  masslessoess  of  the  boson  associated  with  the  elec- 
tromagnetic interactions,  the  Sf.^(2)g  x symmetry  must  broken  while 

preserving  the  Cf(l)  symmetry  of  QED.  Tb  begin  we  introduce,  to  the  uriginal 
Lugrangian.  an  addiiioDaJ  Sf'(2)g  x C/(l)y  invariant  term  such  as  as 

Cs  = (i3n4')'(D'‘4>)-V($).  (1.14) 

This  term  is  knowu  as  the  Higgs  term  or  Higgs  sector.  The  Hig^  field  4>  is  a 
complex  two-component  scalar  field 

which  transforms  as  weak  iao-doublet.  To  allow  for  the  required  spontaneous 
breakdown  of  the  electroweak  symmetry,  we  let  the  potential  {the  Higgs  po- 
tential) take  the  form  of 

The  potential  describes  the  interaction  of  four  scalar  particles  (two  complex 
componeuts  for  each  member  of  the  doublet)  interacting  with  the  four  gauge 
bosons  W and  the  B.  For  positive  aud  A the  minimum  of  the  potential 
gives  the  \‘acuuro  expectation  value  (VEV). 


Wc  uow  want  to  expand  about  the  minima  of  the  potential  and  we  do  so  by 
picking  a particular  configuration  of  the  physical  \acuum  from  among  the  many' 
possibilities.  The  particular  vacuum  configuratioD  is  not  chosen  arbitrarily, 
instead  a configuration  is  picked  such  that  the  original  SU(2)  x U{1]  symmetry 
is  broken  down  to  the  U(l)m  of  QED.  The  choice  made,  expanded  about  the 


Since  4>o  has  T = J and  7^  = -J,  the  generators  for  both  the  S£'{2)  and  U(l) 
are  broken.  But  since  4'n  is  neutral  the  generator  of  QED 

Q = t5  + ^ 

remaius  unbroken  and  thus  tbe  vacuum  is  invariant  under  the  £/(!)«„.  This 
tlien  implies  that  the  gauge  boson  associated  with  this  symmetry  remains  mass- 
less and  we  get  the  massless  gauge  boson  sought. 

The  HjgM  tiiechanism,  To  see  the  explicit  mass  generation  through  the 
Higgs  mechanism,  we  introduce  four  new  fields  £(i)  and  fi(i)  and  express  tbe 
dcctroweak  vacuum  in  terms  of  the  new  fields 


This  particular  form  allow  us  to  gauge  away  the  three  massless  Nambu-Goldstone 
bosons  arising  from  the  spnntaiieously  broken  continuous  symmetry  (the  Gold- 
stone  theorem).  In  gauging  away  the  extra  degrees  of  freedom  tbe  massless 
bosons  acquire  mass  terms  when  the  fields  arc  transformed  by  a unitary  trans- 
formation of  the  form 


t/(0  = 


(1-17) 


The  new  transformed  fields  are  defined  by 

= = where  , = (1.18) 

K = /; = U 

= m (^)  £/-'«)  - il8„U(0)t/-'(€) 

Bl  = Bu.  (1.19) 

This  proeedure  is  known  as  the  Higgs  mechanism-  The  three  degrees  of  free- 
dom associated  with  the  massless  Goldstone  bosoiis  arc  transformed  into  the 
longitudinal  polarisation  associated  with  the  now  massive  vector  bceon  of  tbe 
broken  theory. 

We  DOW  show  tbe  emergence  of  the  mnss  terms  in  the  Higgs  sector  of  the 
Lagranglan  by  expressing  Equation  (l.H)  in  terms  of  the  new  fields 

^3  = (0,<$')’1D“*')-V'{$').  (1.20) 

The  kinetic  energy  term 

(a,. + • w;  +4-5  • w; + 4> 

= ^aa/i(r)9"Mj)-l-i[ff2(Wi‘  -ilV^)(H/‘'‘-hm'^) 

(1.21) 

contains  the  anticipated  mass  terms  for  the  gauge  vector  fields.  This  can  be 
shown  if  we  redefine  W'!  and  as 


(1.22) 


■4^ 


uggiug  all  of  this  ioto  the  kinctk  eoergy  terms  iu  Equation  (1.21)  gi 

i8„h(T)a^A(^)+im=Z'-Z„  + 


(1.25) 


mi,  = v/aX  = iiy/2 


£4  = Y('>£i,e5,$  + + Y'‘^Qi4$  (1.26) 


.eti.Qi  Olid  (efl,ufl.dfl)aretb 


21 

^ukawH  matrices.  The  states  so  far  discussctl  are  gauge  eigenstates  and  thus 
the  Yukawa  mass  matrices  y„  are  iu  general  not  necessarily  diagonal.  Since 
any  matrix  can  be  expressed  as  the  product  of  two  unitary'  matrices  times  a 
real  diagoiiai  matrix  (biunltary  transformation),  we  can  express  the  Yukawa 
mass  matrices  as  Equation  (1.26)  as 
Y<')  = 

Yl^'  = SjMdf'^Rrf, 

Y“-'  = sJ'Md'“>R„ 

In  the  leptoa  sector  the  biunitary  transformations  can  be  absorbed  in  the 
definition  of  the  fieids  by  redebning  them  as 

L[->SiLi.  e'„ -tRjCg-  (1.27) 

In  the  quark  sector  the  right-hand  quark  Belds  can  be  simiiarly  redefined; 
however,  the  left-haitded  quark  doubiets  couple  to  both  up  and  down  types  of 
right-handed  quarks  thus  redefining  say 

Q[~*SjQ,.  (1.28) 

wilJ  introduce  a S'  in  the  right-handed  up  type  coupling  tertu.  The  usual  pro- 
cedure is  to  redefine  the  left-handed  fields  as  in  Equation  (1.28)  and  introduce 
a matrix  to  fix  up  the  right-hand  quark  coupling-  The  matrix  V = is 
known  as  the  Cabibbo-Kobayashi-Moskawa  (CKM)  quark  mixing  matrix 


and  arises  when  the  Yukawa  mnirires  are  diagonaiised.  The  matrix  has  four 
independent  parameters  and  can  be  described  by  three  mixing  angles  and  a 


complex  phase.  A non-zero  value  for  the  phase  is  the  source  of  CP-violatioa 
in  the  Standard  Model. 

After  diagonalizing  the  mass  matrix,  the  form  of  the  \*ukawa  coupling  is 
£4  = (1.30) 

Ehtpanding  Etjuaiion  (1.30)  in  the  unitary  gauge  we  see  that  the  VEV  of  the 
Higgs  doublet  gives  masses  to  fermions  in  a universal  way, 

£«)  = [rnae  + + m.rr] 

^^4^  “ [niijdtl -I- m,8S  + mj66j  (1-31) 

^4*^  * ^1  -I-  (muuu  4-  nifCC  + mr£t| 

with  the  identifications 

^p“(()  = (me.m^.tti,) 

^V"{d)  = (m,i.m„tin)  (1.32) 

^p"(u)  = (ma,m4,m() 

The  tree  level  Feynman  rules  can  be  read  off  the  Classical  Lagranglan 
expressed  as  the  stun  of  Ci  through  £,.  The  Higgs  sector,  however,  increases 
the  number  of  couplings  and  thus  the  total  number  of  diagrams  allovred.  In 
particular,  the  newly  introduced  Higgs  scalar  couples  not  only  to  the  vector 
bosons  but  to  all  the  fermions  as  well.  Also,  additional  couplings  of  the  various 
vector  bosons  to  each  other  make  for  somewhat  more  topologically  complicated 
Fevmnan  diagrams. 


CHAPTER  2 

THE  PHYSICS  OF  THE  b QVARK 

hilrodiiction 

Experinieotal  b physics  began  in  1977  with  the  discovery  of  a narrow  res- 
onance at  9.5  GeV  in  p -t-  A'  ^ + X evenia  at  Fbrmilab's  CFS  (7j  ex- 

periment. Shortly  Ihereatier.  the  observation  was  confirmed  by  e^e”  expec- 
imenl*  including  the  CLEO  [8]  and  CUSB  [9]  collaboration  at  CESR.  These 
experiments  also  measured  the  hadronic  cross-section  in  the  region  arotind 
y/i  c=  10  GeV  and  found  a series  of  narrow  resonance  consistent  with  radial 
oxcitalbns  of  the  hi  quarkonium  state.  The  states  formed  in  e+c“  annihila- 
tion share  the  quauimn  numbers  of  the  parent  virtual  photon  and  are  thus 
produced  in  a J = 1,S  = 0 state.  These  stales  are  known  collectively  as  the 
T resonances.  This  discovery  was  followed  shortly  by  the  discovery  of  hadrons 
comainiug  the  new  6 quark  paired  with  a Ugbt  quark  [10,11].  These  states  are 
known  as  the  fl  mesons  and  measurement  of  their  decay  rates  is  the  primary 
aim  of  this  analysis. 

In  this  chapter  w#  sketch  the  theoretical  framework  which  describes  the 
deetty  mechanism  of  the  B mesons.  We  will  begin  by  briefly  describing  some 
of  the  more  important  points  in  the  physics  of  the  upsilon  system  since  all  B 
meson  used  in  this  analysis  ate  decay  products  of  the  T(4S). 


Figure  2-1  Hadrouic  cross-section  of  e+c”  iumihUstion. 

2.1  The  T Resonances 

In  Figure  2.1  the  hadronic  cross-sectiuo  is  pbtced  os  a function  of  the  c'*'c” 
center  of  mass  energy.  The  plot  shows  the  first  four  T resonances.  The  first 
three  states,  T(1S)  - T(3S),  show  the  characteristic  narrow  widths  expected 
for  quarkonium  states  with  masses  below  open  flavor  threshold.  In  fact,  the 
widths  here  are  dominated  by  the  beam-energy  resolution  rather  than  by  the 
intrinsic  widths  of  the  stales  themselves-  The  widths  and  masses  of  the  states 
aregiveo  in  Table  2.1. 

We  can  understand  the  small  widths  if  we  consider  all  possible  decay  di- 
agrams accessible  to  states  with  masses  below  BB  threshold-  The  possible 
diagrams  are  shown  in  Figure  2.2.  lu  (a)  and  (b)  the  decay  rate,  and  thus  the 
overall  width,  is  reduced,  relative  to  open  bottom  production,  by  higher  order 


Table  2.1 


Resonance 

T(1S) 

T(2S) 

T(25) 

T(4S) 

T(5S) 

T(11020) 


9.4603V  ±0.00021 
10.02330  ±0.00031 
10.3363  ±0.0003 
10.3SOO±  0.0035 


11.019±0.008 


(fev) 

1.32  ±0.03 
0.554  ±0.030 
0.45  ±0.04 
0.24  ± 0.05 
0.31  ± 0.07 
0.13  ±0.030 


26.3  ±3.5 
(23-8  ±1.8)  X 10^ 
(110±13)  X 10^ 
(79±16|  X 10^ 


2.1.1  SpectroBCODV  of  the  T System 

The  T resonances  also  decay  into  a rich  spectniin  of  states  through  strong 
and  electromagnetic  processes.  The  spectroscopj’  of  heavy  flavor  quarkonia  can 
be  described  in  close  analog!’  «'ilh  that  of  elws.romegnetically  bound  positron- 
ium.  Because  of  the  large  mass  of  the  6 quark,  the  T system  can  in  principle  by 
described  by  noaielativistic  QCD  potential  models  |12,13|-  Recently,  the  sym- 
metries of  heaxT  quark  effective  theory  (HQET)  have  been  applied  to  bottom 
and  cbarmonium  systems  [14|. 

Several  spectroscopic  tran.dtious  between  the  various  T(nS)  states  have 
been  observed.  These  include  single  photon  transitions  to  the  states  known 
as  xir  hadronic  transitions  to  and  hadronic  transitions  to  bwer  lying  radial 
excitations.  The  spectroscopy  of  the  T system  is  shown  in  Figure  2.3. 


T(I  lOiOi 


T flOS60) 


Figure  2.3  Tbe  specfrum  of  hottomouiuni. 


Thp  Bret  T resonaDce  massive  enough  to  aJIoiv  for  the  produrtion  of  BB 
mesons  is  the  T(4S).  Tlie  large  width  and  tlie  dranmlie  increase  in  the  produc- 
tion ofleptona  gave  the  first  direct  evidence  for  open  bottom  meson  production 
111].  The  lepton  spectrum  observed  was  cousistent  with  production  in  decays 
of  heavy  quark  through  the  B Xtu  channel  [lOj.  Eventually  the  B me- 
son was  reconstructed  in  exclusive  channels  and  its  mass  was  measured  to  be 
5.2786  ± 0.0008  GeV  [ID]. 

The  T(45)  Is  only  20  MeV  above  BB  threshold.  Due  to  phase-space  limi- 
tations it  can  only  decay  to  bottom  mesons  with  either  a u or  d quark.  Also, 
the  Bret  excited  stale,  the  B‘.  is  40  MeV  above  the  ground  state;  excluding  the 
possibility  of  BB  decays  other  than  to  those  which  are  back  to  back  decays. 
The  small  amount  of  energy  above  threshold  also  implies  that  the  B mesons  are 
produced  nearly  at  rest.  The  average  momentum  of  each  S is  approximately 
330  MeV,  a fact  which  will  be  exploited  when  reconstructing  the  exclusive 
decays  of  the  B. 

Tb  calculate  the  brandling  fractions  of  B mesons  we  need  to  koov/  tbe 
number  of  B mesons  produced.  The  ratio  of  charged  to  neutral  BB  produc- 
tion (/+-//00I  and  the  fractiou  of  BB  pairs  produced  (fgg)  iu  the  decay 
of  the  T(4S)  are  thus  important  quantities  to  consider.  Recently,  the  CLEO 
collaboration  has  determined  the  production  ratio  to  be  [15] 

_ T(4S)- 


fo 


= 1.04  ±0.13  ±0-16  . 


(2.1) 


The  systematic  er 
Tbe  result  is  cons 


T(45)  - 

includes  uncertainties  iu  the  ratio  of  to  B+  lifetimes. 
:nt  with  equal  production  of  B*B~  and  B°B®  pairs.  In 


this  aoals-sis  we  assume  r.bat  /m  = /+-  = Also,  CLEO  11  has  receotb'  sM 
an  uppcr-Umit  on  the  iiuiul>er  of  nou-flfl  decays  of  the  T(45)  by  comparing 
dilepton  and  single  lepton  yields.  The  value  obtained  is  less  than  0.05  at  the 
95%  confidence  level  and  Is  Ihus  consistent  with  /gg  = 1 (16]. 

2.2  The  Decay  of  tlie  /?  Meson 

Weak  decays  of  B mesons  have  been  the  subjects  of  great  theoretical  and 
experimental  interest  since  their  discovery  in  the  early  1980s.  The  large  mass 
of  the  6 quark  and  the  large  momentum  transfer  involved  in  its  decoy  provides 
an  excellent  laboratory  to  probe  both  the  weak  and  strong  interactions.  The- 
oretical models  which  describe  tbe  deeny  processes  of  badrons  can  be  better 
tested  at  the  large  energy  scales  typical  of  6 quark  decays- 

Therc  are  significant  theoretical  simplifications  that  arise  in  the  study  of 
QCD  by  considering  mesons  witli  large  masses.  In  particular,  the  large  mass 
of  tbe  heavy  quark,  when  compared  to  AqcD-  implies  that  its  actual  value 
becomes  irrelevant  to  the  dynamics  of  tbe  light  degrees  of  freedom  in  the  com- 
posite hadron.  The  heavy  quark  can  thus  be  thought  of  as  a static  source  of 
color  which  interacts  with  tbe  remaining  light  quarks  and  gluons.  In  the  infinite 
mass  limit  the  pnrtinilar  flavor  of  the  quark  and  its  spin  also  become  irrelevant 
in  the  dynamical  behavior  of  the  hadron  and  give  rise  to  SC/(2A')  spin-flavor 
symmetries.  These  ideas  form  the  genera!  framework  upon  which  the  Heavy 
Quark  EBeciive  Theory  (HQET)  is  batted  [17).  In  HQET  the  spin-flavor  sym- 
metry is  exploited  to  build  an  effective  Lagrangian  density.  The  Lagrangian 
is  then  nsed  to  derive  an  effective  theory  analogous  to  QCD,  complete  with 
its  own  set  of  Feynman  rules,  in  which  tlie  perturbative  expansion  is  also  now 
an  expansion  in  Aqcd/"1o-  These  ideas  have  generated  many  new  theoretical 


Figure  2.4  Possible  B mcsoo  weak  decay  diagrams. 

developraeats  in  the  past  decade  whicb  have  led  to  a belter  understanding  of 
the  processes  that  govern  the  decays  of  heavy  quarks. 

The  remainder  of  this  chapter  is  devoted  to  the  phenomenological  descrip- 
tion of  the  weak  decay  of  heavy  flavor  mesons.  While  all  of  these  ideas  apply 
to  charmed  mesons  as  well,  we  will  concentrate  only  on  the  decays  of  the  B 
meson.  Also,  this  analysts  deals  with  the  exclusive  reconstruction  of  two-body 
B meson  decays.  Tbe  scope  will  thus  be  limited  to  the  discussion  of  exclusive 
Don-leptouic  deci^s. 

Tb  begin  the  description  of  the  decoys  of  the  B mesons,  we  first  consider  all 
tlie  possible  weak  decay  diagrams  which  can  mediate  tlie  process.  In  Figure  2.4 
all  possible  lowest  order  weak  diagrams  which  contribute  to  the  two-body  decay 
of  the  B uiesou  are  showo.  In  (a),  the  external  spectator,  the  quarks  produced 


at  the  K‘  \'ertex  hadronije  logethtr  into  one  of  the  two  final  states-  This  forms 
the  simplest  of  the  four  decays  diagraixis  to  treat  theoretically  since  it  involves 
no  qq  popping  and  the  final  states  consist  of  matched  qnark  colors  fields.  In 
(b)  the  internal  spectator  or  "color-suppressed"  diagram  things  are  a bit  more 
complicated.  Now  tbe  original  light  spectator  quark  hadrouites  with  one  of 
the  qtiarks  produced  at  the  ll  vertex.  Hard  gluon  exchanges  can  no  longer  be 
ignored  and  must  uow  bo  included  in  the  theoretical  dcscriptbn.  There  is  also 
a requirement  that  tbe  color  degrees  of  freedom  match  within  each  composite 
hadron.  This  requirement  wouid  naiveiy  predict  that  processes  wliich  can 
proceed  only  through  internal  diagrams  ate  suppressed  relative  to  ones  which 
proceed  only  through  external  diagrams.  Since  there  are  three  colors  to  include, 
the  naive  expectation  gives  a suppression  factor  of  (l/Nj)*  = 1/9.  We  will 
see  later  that  this  situation  Is  somewhat  more  complicated  when  hard  gluons 


The  other  diagrams  in  Figure  2.4  are  usually  ignored  since  their  contri- 
bution to  the  decay  modes  considered  are  small  when  compared  to  simple 
spectator  processes  [18].  For  example,  the  contribution  from  tbe  exchange  dia- 
gram (d)  is  both  heiicity  and  color  suppressed  relative  to  the  external  spectator 
(a).  Tlie  overall  decay  rate  is  reduced  by  a factor  -^(m? -l-mj)^,  where  fg  is 
the  B meson  decay  constant,  mu  is  the  mass  of  the  light  quarks  (a  u or  a d) 
and  ms  is  the  mass  of  the  h quark.  The  expected  suppression  is  10“®  and 
thus  the  usual  procedure  of  ignoring  this  channel  Is  wdi  justified.  The  quark 
combination  possible  in  (d)  do  not  provide  a channel  for  the  decay  modes 
considered  in  this  analysis. 


32 

The  panicular  decay  modes  of  the  B meson  cousidcred  for  this  analysis  fall 
into  three  general  categories: 

S3  -*  Class  I decays 

S3  ■->  M°m'‘  Class  II  decays  (2-2) 

-t  A/°m”  Class  III  decays 

The  M+  or  M°  represcnl  a charmed  pseudoscalar  or  vector  meson  and  the 
m represents  a “light"  pseudoscalar,  vector  or  aaial-veclor  state.  The  light 
mesons  ate  all  part  of  the  family  of  isospin  singlet  or  triplet  states.  Thcrelevanl 
daughters  are: 

M*  = D*.D‘*  and  A/°  = D‘'.D‘‘' 

0 _ 0 ^0 

We  will  see  later  that  the  three  categories  essentially  separate  the  decays  into 
modes  which  can  proceed  tlirough  either  the  external  (class  I)  spectator,  the 
internal  spectator  (class  II)  or  through  a combinailon  of  both  weak  decay 
diagrams  (class  III). 

2-2-1  The  Weak  HamiltOTii.-ir 

The  starting  point  in  the  description  of  tlie  phenomenology  of  non-Ieptonic 
decays  of  heavy-flavored  mesons  is  the  Lagrangiau  density.  The  bare  rurrenl- 
current  form  of  the  Lagrangian  density  for  a CKM  allowed  6 -*  c transition  is 
given  by  the  coupUng  of  the  weak  currcots  to  the  H'*  operators.  The  charged 
current-curtent  Lagrangian  is 


whore  the  ctirreut 

Ji:  = «)'  = (u  c *)7'‘(l-1s)'cA-4,(s] 
are  Ibc  !eft-handed  currents  in  the  Sl/(2)  * J/(l)  standard  model.  The  symbol 
^CKM  denotes  the  Cabibho-Kobayashi-Maskawa  mixing  matrix  of  Eriuation 
(1.29)  and  relates  the  weak  eigenstates  to  mass  eigenstates.  The  values  of 
the  individual  entries  are  not  predicted  in  the  Standard  Model  but  must  be 
determined  hy  experiment  [2). 

Flom  the  Lagrangian  density  we  can  read  off  the  effective  weak  Hamilto- 
nian; which  to  lowest  order  in  the  weak  coupling  constant  is 

= (5^)  d(;')7e(l  - 7sl«i(«)}. 

(2.4) 

This  expression  corresponds  to  the  quark  level  spectator  diagrams  in  Figure 
2.4  where  a 6 quark  decays  to  a c quark  and  the  virtual  H’  creates  a ud  pair. 
For  mq’  > k.  the  denominator  in  the  gauge  propagator  reduces  to  a cemstaat 
and  the  mass  of  the  gauge  boson  is  absorbed  into  the  dimension-full  coupling 
constant  Gf, 

’ ^‘e4K«/{e(p')y‘(l  - 7s)l>(p)^(«')ya(l  - Ts)>'(9))-  (2.5) 

This  limit,  which  corresponds  to  a zero  range  interaction,  is  justified  for  the 
class  of  decays  considered  here  since  the  energy  scale  k of  the  reaction  is  much 
smaller  than  mvp.  the  mass  of  the  vector  boson. 

2.2.2  QCD  Corrections  to  the  Effective  Hamiltonian 

The  form  of  the  effective  Hamlitonian,  Equation  (2.5),  is  significant))'  mod- 
ified by  QCD  corrections.  To  lowest  order  in  o*  the  effects  are  included  by 
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Figure  2-5  QCD  correcteO  simple  spectator  decay  diagram. 

adding  all  topologically  distinct  weak  decay  diagrams  which  iuclude  the  ex- 
change of  a single  gluon.  At  this  point  in  the  derivation  an  assumption  is 
made  about  long  distance  QCD  effects.  The  assumption  is  simply  that  soft 
gluon  exchange  do  not  significantly  affect  the  decay  of  the  heavy  mesons  and 
are  thus  ignored  in  the  calculation.  Also,  the  energy  scale  is  assumed  to  be  suf- 
ficiently large  sucb  that  higher  order  corrections  give  uegiigibic  contributions. 

The  possible  Feynman  diagrams  which  contribute  to  first  order  in  a,  are 
shown  in  Figure  2.5.  The  diagrams  are  separated  Into  (a)  those  that  induce 
new  effective  4-Fermi  interactions  with  different  color  structure  and  (h)  those 
that  comribute  to  the  renormalization  of  the  weak  couplings,  quark  masses  and 
wavefunctions.  In  addition  there  are  penguin  type  contributions  which  are  of 


Figure  2.6  Hard  giuon  corretied  simple  spectator  diagram. 

the  same  order  in  a,  hut  lead  to  small  contributions  and  are  typically  ignored 
for  the  decays  considered  here  [18]. 

To  see  the  origin  of  the  new  effective  4-Feruii  interactions  we  first  write 
down  the  effective  Hamiltonian  with  explicit  color  indices-  The  Lorentz  indices 
are  suppressed  In  the  expression  and  summation  of  the  repealed  indices  is 


3(i 

implied.  The  effective  Hamiltonian  is 

' (*<•«■ ’'-I"}- 

[2.4) 

It  corresponds  to  the  first  of  the  four  diagrams  in  figure  2.6(a)  which  has  been 
recreated  with  explicit  Lorcnts  and  color  indices  in  Figure  2.6. 

To  evaluate  the  integral  in  Equation  (2.4)  we  first  simplify  the  expression 
by  setting  the  quark  masses  and  external  momenta  to  lero  and  collapse  the  W 
propagator  to  a point.  The  integral  can  now  be  evaluated  from  down  to 
an  arbitrary  cutoff  scale  p |19).  We  then  sum  the  four  diagrams  in  2.5(a)  and 
simplify  the  result  by  doing  a bit  of  gamma  matrix  algebra  [20) 

'“®  (^)  - 15)^Sr“c}  • 

(2.5) 

The  derivation  above  shows  how  hard  gluons,  exchanged  between  the  two 
quark  currents,  introduced  a A"  matrix  to  each  curreut  in  the  current-current 
operator.  To  show  explicitly  the  generation  of  the  new  effective-neutral  term 
we  perform  a Fiera  transformation  ou  Equation  (2.5)  by  using  the  identity 

(•*o)ij  (•^),T  = ”1  4-  26,,6j,.  (2,6) 

We  now  have  Equation  (2.5)  in  terms  of  effective  charged  and  neutral  current 

„ Gp3o,.  (Mt\  .(  0 

° ■;;^8T  I ~ I - j c,7'(l  - 75)5'  - rf.7r(l  - 75)n‘ 

-F25,7‘'(1-75)u‘  d;7e(l-7s)5’}. 

(2.7) 
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ir  we  collect  terms  with  definite  Savor  symmetries  vre  can  rewrite  Equation 
(2.7)  in  terms  of  operators  which  belong  to  different  flavor  representations  and 
thus  do  not  mix  niider  renormalization.  The  new  operators  (O*)  are  given  by 

0*  = | [ci'(l  -05)4- d-,„(l--j5)a± 5)6]  (2.8) 


Expressed  in  terms  of  0±,  Hgff  now  reads 

"e//  = ud  {C+(p)0+  + C-(p)0-} . (2.9) 


The  scale  dependent  coefficients  C±(p)  provide  a convenient  place  to  store  the 
effects  trom  renormalization  and  hard  gluon  correctioas.  Before  the  renormal- 
ization procedure  the  coefficients  C±  are  given  by 


C+fp)  = 1 - ^ log  j 

C^(p)  = l + ^log^i^j. 


(2-10) 


Use  of  the  renormalization  group  equations  gives 
approximation  (LLA)  at  the  6 quark  scale  [21). 


n the  leading  logarithmic 


where  7 27^.  ,6=11-  and  = 4 for  6 decays.  With  p = mj  =; 

5 GeV  and  \qcd  — 0.250  GeV  we  find 


C+  =0.85  C-  = 1.39  (LLA) 

(2,12) 

C*  =0.82  C-  = 1.49  (NLLA) 

The  next-to-leading  log  (NLLA)  corrections  modify  the  (LLA)  results  only 
slightly  [18]  indicating  that  at  the  6 scale  the  perturbative  expansion  gives 


reliable  results. 


2-2.3  Modela  of  Nnii-lpt^tonic  B Meson  Dpruv.q 

The  previous  section  concluded  with  the  QCD  corrected  weak  Hamiltonian. 
The  derivation  showed  how  QCD  corrections  modify  the  bare  weak  Hamilto- 
nian and  how  the  effective  neutral  current  term  arises  from  the  exchange  of 
hard  gluons  hetweeu  different  quark  currents.  Also,  the  derivation  showed  that 
two  scale  dependent  coefficients,  each  multiplymg  its  own  effective  weak  cur- 
rents, absorb  all  relevant  QCD  effects.  In  this  section  we  describe  how  the 
effective  Hamiltonian  is  used  to  find  the  amplitude  for  the  decay  of  the  heavy 


We  hegiu  by  rewriting  the  Hamiltonian,  Equntiou  (2.9).  with  ail  effective 
charged  and  neutral  currents  terms  grouped  together.  The  Hamiltonian  is 


are  known  as  the  scale  dependent  Wilson  coefficients  |22], 

Each  term  in  Equation  (2.13)  corresponds  to  one  of  the  valence  quark  di- 
^ams  shown  in  Figure  2.7.  Parts  (a)  and  (b)  rortespond  to  the  fiist  and 
second  term  in  Equation  (2.13)  respectively.  The  quark  color  structure  of  the 
first  term  forms  combioatioDs  which  are  color-singlets  and  thus  lead  directly 
to  teal  hadroas.  In  the  second  term  the  color  degrees  of  freedom  are  mixed. 
We  can  obtain  a color-singlet  flavor  combination  by  perfnrmiug  another  Fierz 
transformation  on  the  color  mixed  terra.  Applying  the  transformation.  Equa- 
tion (2.6),  to  the  second  term  in  Equation  (2.13)  projects  out  a color-singlet 
which  now  gets  multiplied  by  l/Nc  and  introduces  an  octet  term  O].  The 
effective  Hamiltonian  is  now 


^e//  = [<^->'(i  --rs)6'jp„(i  - ,5)u»] 


2 


(2.14) 


and  Oiocc>‘’(l— «)^y^Ad-,<,(l— >b)(y)u  (2.16) 


o]  = Cl  -e  — Cj 


The  preceding  dUcusaiou  involved  class  I decay  channels  were  the  external 
spectator  diagram  dominates.  In  a similar  way,  the  amplitude  for  class  II 
decays  can  be  obtained  by  relabeling  the  valence  quarks  in  Figure  2.7  (see 
Figure  2.8)  to  reBect  the  0nvor  combinations  of  an  intcroal  spectator  process. 
Rearrangement  of  the  quark  8elds  gives 

= +Oj|,  (2,27) 
which  now  a represents  the  quark  flavor  and  color  structure  in  class  II  decays. 
The  new  coefliciciit  02  and  operator  O2  are  analogous  to  those  derived  for  class 
1 decays.  They  arc  defined  as 


A-  (2.18) 


In  Equations  (2.15)  and  (2-17)  the  effective  Hamiltonians  are  expressed 
as  functions  of  two  operators,  a coIor-singlet  and  an  color-octet.  Since  we 
ultimately  want  to  evaluate  amplitudes  of  decay  processes  which  involve  treio- 
sitions  between  real  hadronic  states  the  effective  Hamiltonian  expressed  in  this 
way  lends  itself  to  this  task  it  we  And  some  way  to  deal  with  the  octet  term.  If 
the  hadrons  wliere  indeed  pure  singlets  then  the  octet  term  vanishe.s  by  color 
conservation.  However,  hadroms  are  not  pure  singlets,  they  are  instead  compli- 
cated objects  which  consLst  of  many  gluon  and  qq  pairs  which  can  interact  with 
the  octet  term  and  thus  affect  the  overall  amplitude.  In  the  usual  treatment, 
the  hadron  is  assumed  to  be  a color-singlet  and  thus  the  octet  term  can  be 
ignored.  While  this  procedure  is  somewhat  ambiguous  aud  perhaps  unjustiBed 
(1D,23,:8),  it  leads  to  simplifications  which  allow  the  calculations  of  the  matrix 
elements  from  the  effective  Harailtouian  described  above. 

The  factorization  hypothesis  Most  current  models  of  non-leptonic  decays 
employ  the  factorization  assumption  in  which  the  ovprali  amplitude  is  written 
as  the  product  of  two  matrbr  elements.  Each  matrix  element  describes  sepa- 
rately the  transition  between  the  initial  heavy  meson  and  one  of  the  daughters 
and  the  creation,  from  the  vacuum,  of  its  decay  partner.  In  this  way  the  ampli- 
tude can  be  evaluated  by  writing  the  matrix  elements  in  terms  nf  form  faclors 
and  decay  constanus.  These  can  be  either  determined  by  direct  comparison  to 
semi-leptooic  decays  nr  by  theoretical  model  calculations. 

There  arc  heuristic  arguments  that  attempt  a justification  of  the  factor- 
ization hypothesis  in  certain  kinematic  regimes.  At  large  momentum  transfer 
the  expectation  is  that  the  (ud)  quark-pair  produced  at  the  IV  vortex,  which 
at  this  stage  is  a point-like  rolot-singlet,  leaves  the  interaction  region  before 


its  dipolo  moment  becomes  sufficiently  large  to  interact  with  the  color  field 
surrounding  the  spectator  quark  system.  One  can  estimate  the  hadronizatiou 
scale  by  finding  the  time  it  takes  the  (ud)  to  be  1 fm  apart,  the  typical  strong 
iDteraotioQ  range.  While  the  a and  d are  flying  apart  they  are  also  moving 
away  from  the  spectator  quark  system.  If  by  the  time  the  u and  d arc  approtti- 
mately  1 fra  apart  the  pair  is  more  than  1 fm  away  from  the  spectator  system, 
the  ud  pair  is  thought  to  have  escaped  the  effects  of  the  color  field  surrounding 
the  spectator  system.  For  decay  modes  wliere  the  B meson  decays  to 
(class  1)  the  badronization  scale  is  reached  when  the  quark*pair  is  about  20 
ftn  away  from  the  spectator  system  [19).  For  class  II  decays  the  hadronization 
scale  is  reached  sooner,  about  1-2  fm,  so  it  is  not  clear  if  factorization  holds  in 
this  class  of  decays.  Recently  Dugan  and  Grinstein  (24)  have  argued  that  some 
basis  exists  m QCD  for  factorization  in  appropriate  kinematic  regimes. 

evaluate  the  amplitude  of  a particular  reaction  in  the  factorisation  approxima- 
tion. To  clarify  the  discussion  we  consider  class  I modes  where  the  B meson 
decays  to  plus  a “light"  partner.  These  decays  correspond  to  the  H,// 
described  by  Equation  (2.15).  Their  amplitudes  are  given  by 


Under  the  factorization  hypothesis  this  can  written  in  terms  of  the  two 
independent  matrix  elements 


' = -15H0>{DW'“1ci"(1  -15)4100)  (2.19) 


- |{P  + P')" ^;*"^(P-P')*'|fl(?^) 


{D'*(£.p')|c7'‘is6|^(p))  = +<(mfl+mi,.)  [cj  - A:{q^) 

(2-20) 

The  roeffiden^  fb(9’],f,(j^).  V(y^),  „„d  A,(q^)  are  ka<ma.  aa 


DUiabers  (Jf  = 0*  for  Fq.  1"  for  F,  aod  T,  1+  for  /I,  and  A,  and  0"  for  4o) 


The  form  factors  at  7^  ^ 0 are  also  model  dependent  and  various  methods 
have  been  used  to  calculate  them.  In  the  BSW  model  a relativistic  harmonic 
oscillator  potential  is  used  to  expand  the  wavefunction  of  the  initial  aud  BnaJ 
state  mesons.  The  weak  hadroaic  cuirenM  are  Ibeu  written  in  terms  of  creation 
and  annihilation  operators  and  matrix  elements  are  formed.  The  results  are 
then  compared  with  the  parameterization  above  atid  the  form  factors  at  7^  = D 
are  obtained  [25).  Recently  Neubert,  Rieckerl,  Stech  and  Xu  have  calculated 
the  form  factors  by  using  HQET  and  a sUghtly  modified  version  of  the  BSW 
model  |19j. 

The  creation  from  the  vacunni  of  the  pseudoscalar  follows  from  the  Lorentr 
structure  of  the  current  acting  on  the  vacuum. 

(x“|d7„7su|0)  = 

The  ouly  Lorema  4-vector  here  is  the  4-momentum  of  the  plon  and  the  matrix 
element  is  expressed  as  a the  product  of  this  4-vector  and  the  decay  constant 
It-  The  creation  of  a vector  from  the  varitum  is  given  similarly  by 

(f  “MipulO)  = /u<Smii. 

The  decay  constants  are  determined  experimentally  by  measurements  In  e+e“ 
annihilation.  leptonie  decays  or  by  extracting  them  from  other  process  such  as 
T — ♦ bur  decays. 

The  Bauer  Stech  and  Wirbel  anomsch  In  the  discussion  above  factorisa- 
tion was  assumed  to  bold  strictly  as  in  semileptonic  decays  where  the  leptons 
produced  at  the  11 ' vertex  do  not  interact  with  the  color  field.  This  assumption 
allowed  us  to  drop  the  color-octet  term  which  is  as  already  described  thought  to 


be  uot  completely  justifiable.  There  are.  in  additioD,  other  ambiguities  inher- 
ent in  the  procedure  above.  These  involve  the  exact  scale  at  which  to  evaluate 
tlie  QCD  coefficients  [19].  Also  the  extent  to  which  noii-petturbative  effects 
are  ignored  Is  another  source  of  concern  iu  the  usual  tiratmcnls  of  heavy  quark 
decay.  These  ambiguities  lead  Bauer,  Stech  and  Wrbel  [25,26]  to  propose  that 
the  coefficients  0|  and  oj  should  instead  be  ueated  as  free  patainetere  which 
can  be  determined  from  experiment.  In  this  way  processes  that  occur  at  simi- 
lar energy  scales  can  be  compared  to  one  another  and  decay  rates  which  would 
otherwise  be  difficult  to  predict  can  he  obtained.  UsuaUy  the  parameters  o, 

ai  =Ci(p)-r{C2(p) 

(2.22) 

where  ^ = 1/N^  is  the  variable  instead.  By  allowiug  these  coefficients  to  vary 
non-factorisable  contributions  are  given  a means  by  which  to  enter  the  calcula- 
tion without  the  need  to  know  how  to  evaluate  terms  which  contain  color-octet 
currents,  the  energy  scale  to  use  and  bow  to  include  non-perturbative  effects. 


CHAPTER  3 

EXPERIMENTAL  APPARATUS: 
THE  CLEO  II  EXPERIMENT 


Introduction 

The  data  used  in  this  analj-sis  were  collected  at  the  CESR  c+e”  storage 
ring  with  the  CLEO  II  detector.  The  CLEO  collaboration  coosiats  oI  over 
200  meoibere  aasocialed  with  20  institutions  in  both  Canada  and  the  United 
Sutes.  These  scientists  are  primarily  responsible  for  maintenance  of  the  CLEO 
II  detector,  data  taking  and  data  analysis.  The  accelerator  physlcisu  at  CESR 
are  primarily  Involved  in  the  development,  rnaintcnance  and  construction  of  the 
accelerator  machinery.  In  this  chapter  we  will  briefly  describe  the  machinery 
involved  in  providing  the  data  insed  in  the  analysis.  For  a more  complete 
description  refer  to  the  various  published  technical  papers  referred  to  in  the 

CESR:  The  Accelemtor 

Figure  3. 1 shows  a schematic  of  the  Cornell  Electron  Storage  Ring  (CESR) 
accelerator  complex  [27|-  CESR  is  an  electron-positron  collider  with  a circum- 
ference of  758  meters.  Currently,  only  the  Southern  iuteraclion  region  is  in  use. 
It  houses  both  the  CLEO  II  ejtperimeot  the  detector  complex  for  the  Cornell 
High  Enetgj-  Synchrotron  Source  (CHESS).  Tlie  ring  is  operated  in  the  neigh- 
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• Elecuon  Bunch  - Counter  Clockwise 


Tkblo  3.1 


CESR  operaliug  parameters 


Revolution  Frequency 
RF  Frequency 
Beam  Energy  Range 


■I.5-6.D  GeV 


50C  MHa 


Fractional  Beam  Energy  Spread 
Energy  Loss/tuin  at  5-3  GeV 


l.H  MeV 


0.062  % 


Collision  Frequency 
Bunch  Length 
Peak  Luminosity 


2.5  X lo’ 


2.7  MHz 


horhood  of  5.3  GoV  with  a peak  luminosity  of  2.9  x lO^cm'^s"-.  These  and 
other  CESR  parameters  are  listed  in  Table  3.1. 

The  positron  beam  fill  procedure  is  divided  into  a series  of  steps.  First 
electrons  are  hoilerf  off  a filament  wire  and  accelerated  in  the  linear  acceler- 
ator (LINAC).  To  create  the  positrons,  the  electrons  are  directed  towards  a 
tungsten  target  half  way  down  the  LINAC.  The  collision  between  the  electrons 


electromagnetically  separated  from  the  debris  and  accelerated  to  150  MeV  by 
the  LINAC.  The  pcsitrons  are  subsequently  introduced  in  to  the  synchrotron 
where  their  energy  Is  increased  to  the  operating  energy,  about  5.3  GeV.  Fi- 
nally, the  particles  are  injected  into  the  storage  ring  where  they  are  prepared 
for  collisbiB  with  the  electron  and  will  remain  until  the  beams  rie;.rad..|»n 
warrants  another  fill.  The  electrons  ate  accelerated  in  the  same  sequence  of 
steps  but  the  conversion  process  is  omitted.  CESR  is  operated  with  seven 
counter-rotating  bunches  of  electrons  and  positrons.  The  bunch  size  is  needle 


which  the  positrons  are 


!)baped,  0.16  by  0.011  mm  and  17  cm  long.  Each  bunch  comains  approximately 
1.4  X 10**  particles. 

During  injection  the  bunches  are  kept  apart  by  electrostatic  separators. 
The  separators  shepherd  the  electrons  and  positrons  into  separate  **prets:pr' 
like  circular  orbits  which  ate  not  allowed  to  intersect  at  any  point.  After 
the  injection  procedure  is  complete,  the  .separators  are  powered  down  and  the 
bunched  orbits  are  allowed  to  intersect  ar  the  South  interaction  region. 

The  curved  path  of  the  electrons  and  positrons  can  then  radiate  energy, 
a phenomena  known  as  synchrotron  radiation.  For  highly  relativist  electrons, 
the  energy  radiated  per  particle  per  turn  is  given  by 


where  e is  Che  charge  of  Che  electron,  p is  the  radius  of  the  ring  aud  E is  the 
ener^  of  the  electron.  At  typical  CESR  energies,  the  losses  due  Co  synchrotron 
radiation  amount  to  about  0.5  MeV  per  turn.  To  compensate  for  the  losses, 
the  storage  ring  uses  two  500  MHz  radio-ftequency  (rf)  cavities  which  provide 


The  CLEO  II  detector  is  a soleuoidal  detector  Chat  combines  excelleut 
charged  particle  detection  with  a high  resolution  electromagnetic  calorimeter. 
It  was  installed  in  the  fall  of  1988  and  superseded  the  existing  CLEO  detector 
with  improved  time<or-flighc  detectors,  muon  detectors  and  added  the  cesium 
iodide  (Csl)  crystal  calorimeter.  The  concentric  detector  design  consists  of 
three  drift  chambers,  a time-of-flight  system,  the  crystal  calorimeter  and  the 


a “kick''  Co  the  electrons  and  po 


ng  their  lost  energy. 


3.2  The  CLFO  IT  Detector  Hardware 


Figure  3,2 


II  det« 
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Figure  S.3  Endcap  view  with  several  different  position  along  the  : direction. 


muon  detector.  All  but  the  muon  detector  element  are  Immersed  In  a 1.5  tesla 
superconducting  coil  1.5  meters  in  radius. 

The  complete  detector  is  thrnugbly  described  in  Ref.  128].  lu  this  chapter 
we  will  brieiy  sketch  the  various  components.  The  elements  will  he  described 
in  order  of  radial  distance.  A schematic  dravting  of  the  whole  detector,  in 
cross-section  and  in  end  view  is  shown  in  Rgurcs  3.2  and  3.3. 

To  separate  the  detector  volume  from  the  vacuum  system  of  CESR  a thin- 
walled  beryllium  beam  pipe  surrounds  the  interaction  region.  The  beam  pipe  is 
33  cm  long  and  has  a radius  of  3.5  cm.  To  reduce  the  incidence  of  synchrotron 
radiatiou  on  the  detector  volume,  the  inner  wall  of  the  beam  pipe  is  coated 
with  a thin  layer  of  sliver.  25mi(Tons  thick  and  an  even  thinner  layer  of  nickel, 
less  than  1 micron  thick.  The  thickness  of  the  beam  pipe  is  is  500  ftm  and 
presents  a total  of  0.45t  radiation  length  of  material  in  the  radial  direction. 

3.2.2  Chaiecd  Particle  Tracking  System 

Charge  particle  trajectories  are  measured  in  CLEO  II  by  a series  of  three 
concentric  wire  drift  chambers.  The  common  axis  is  aligned  along  the  beam 
direction.  The  three  drift  chamber  achieve  different  goals  which  would  other- 
wise be  difficult  to  achieve  by  a single  detector  element.  Each  drift  chamber 
provides  precise  track  position  measurements  by  recording  the  drift-time  of 
electrons  which  arrive  at  the  sense  anode  wire.  As  a charged  particle  travels 
through  the  drift  chamber,  it  ionizes  the  atoms  in  the  gas  filled  volume.  The 
electrons  released  in  the  ioniaation  process  drift  in  the  electric  field  produced 


by  the  field  wires  nnd  are  collected  by  the  sense  wire.  Knowledge  of  the  elec- 
tron drift-time  to  distance  allows  for  a precise  determination  of  the  particle's 
path-  The  curvattire  of  the  path  and  knowledge  of  the  axial  magnetic  field 
are  used  to  find  the  momentum  of  the  particle.  The  two  outer  drift  chambers 
also  provide  inforniatiou  on  the  r position  of  the  track-  Specific  ionization 
energy  loss  measurements,  from  the  outer  drift  chamber  arc  used  for  particle 
identification. 

All  three  tracking  cbamlters  and  the  muon  detector  ere  fed  a coustant 
supply  of  argon-ethane  gas.  The  50%  argon  50%  ethane  (C2H8)  mixture  is 
delivered  by  a complicated  gas  circulation  system  design  to  maintain  a suble 
fiow  and  operating  pressure  to  tbe  various  detector  subsystems  while  providing 
continuous  monitoring  of  flow  rates  and  gas  compositious.  In  the  Spring  of 
1092.  the  argon-ethtmc  gas  mixture  iii  the  precision  tracking  layer  was  replaced 
by  di-methyl  ethane  (DME)  gas  to  improve  tracking  re.solutiou.  The  change 
occured  between  the  first  and  second  half  of  the  data  sample  used  in  this 
analysis. 

The  precision  tracking  layers  IPTLI.  The  layer  closest  to  the  beam  pipe  is 
the  precision  tracking  (PTL).  It  provides  precise  measurements  of  transverse 
panicle  position  near  the  interaction  region  to  determine  particle  directions 
and  to  separate  primary  from  secondary  vertices. 

The  PTL  is  a six-layer  straw  tube  chamber  with  64  axial  wires  per  layer. 
Each  layer  is  staggered  by  a half  cell  from  the  adjacent  layer.  The  field  cage 
for  each  tube  is  defined  by  an  aluminized  Mylar  tube  instead  of  cathode  wires 
utilized  by  the  two  outer  drift  chambers.  Within  each  Mylar  tube  lies  a gold- 
plated  tungsten  sense  wire,  15  pm  in  diameter.  The  single  hit  resolution  of 


wires  i«  50  Am  rms  with  di-methyJ  gas.  The  diameter  of  each  tube  is  determined 
by  its  distance  from  the  interaction  point  sucb  that  eaeli  tube  makes  contact 
with  all  neighboring  tubes  in  adjacent  layers.  The  tubes  are  glued  together  for 
mechanical  stability  and  internal  alignment.  The  detector  assembly,  including 
the  inner  and  outer  walls,  is  33  cm  long  and  extends  radially  from  4.5  cm,  just 
outside  the  beam  pipe,  to  8.1  cm.  A schematic  of  the  PTL  is  show  iu  figure 
3-4.  Neither  polar  angle  nor  track  r position  measurements  are  provided  by 
this  detector  element. 

The  vertex  detector  fVDl.  The  intermediate  drift  ejiambcr.  originally  in- 
stalled as  a the  vertex  detector  in  the  CLEO  I experiment,  provides  for  both 
azimuthal  position  measurements  and  position  measurements  along  the  r di- 
rection. The  spacing  and  cell  widths  in  this  detector  element  are  70%  smaller 
than  the  outer  drift  chamber  providing  better  granularity  for  track  separation. 
The  detector  also  provides  for  position  measurements  along  the  : direction 
with  two  layers  of  segmented  cathodes  readouts  and  charge  division  on  all 
sense  wires- 

TVansverse  particle  momenta  are  determined  by  measurements  of  track 
azimuthal  position  collected  with  a total  of  800  sense  wires.  The  sense  wires 
are  interspersed  with  2272  field  wires  which  define  hexagonal  cell  boundaries 
around  each  field  wire.  The  cells  are  arranged  into  10  axial  layers  of  various 
sizes  to  completely  fill  the  volume.  The  layers  are  divided  into  two  groups.  The 
Inner  group,  layers  1 through  5,  hove  64  cells  per  layer  while  the  outer,  layers 
6 through  10,  have  96  cells  per  l^er.  All  cells  between  layera  are  staggered 
by  half  a cell  within  each  group.  The  wire  configuration  of  the  VD  is  sliowii 
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Bcbematically  Id  Figure  3.4.  The  single  hit  resolution  is  90  /im  in  the  (r  - 
plane,  and  has  a single  hit  efficiency  of  95%. 

The  inner  and  outer  surfaces  consist  of  segmeuted  caCiiodes  wbicli  define 
the  inner  and  outer  field  cage.  The  segmentation  is  5.85  {6.85)  mm  along  the 
beam  direction  on  the  inner  and  (outer)  cathode  surface.  The  sense  wires  are 
made  of  a nickel-chromium  alloy  witli  about  three  limes  the  resistivity  of  gold 
plated  tungsten  and  arc  instrumented  for  charge  division  mcasuremeuls.  The 
segmented  cathodes  and  charge  divisiou  provide  for  polar  track  measurements. 
The  resolution  of  the  cathodes  is  750  ;tm. 

The  drift  chamber  tPR.1.  The  outer  drift  chamber  (DR)  provides  track 
position  measurements  in  both  the  azimutli  and  : directions  (29).  These  niea- 
suremeuts,  together  with  measurements  obtained  with  both  the  VD  and  PTL, 
are  used  to  determine  the  particle's  transverse  momentum  (jn).  the  radial  dis- 
tance of  closest  approach  of  the  track  projected  back  to  the  interaction  point 
and  the  asimuthai  («)  direction.  Tlie  outer  drift  diamber  alone  is  used  to 
measure  the  specific  ionisation  energy  loss  {dE/dx}  for  particle  identification. 

The  outer  drift  chamber  consists  of  small  rectangular  cells,  all  of  nearly 
equal  size,  that  fill  the  entire  volume.  A total  of  12240  sense  wires  are  inter- 
spersed with  36240  field  wires.  The  cells  are  arranged  In  a pattern  of  51  layers 
with  three  field  wires  for  eocli  sense  wire.  The  layers  are  grouped  in  fives  or 
threes,  with  equal  number  of  wires  per  layer  within  each  group.  To  resolve  left- 
right  track  sign  ambiguity,  the  cells  are  staggered  by  half  a cell  with  respect  to 
Its  uearesi  neighbor  in  tbe  radial  direction  within  each  group.  The  groups  of 
five  or  three  layers  are  separated  by  a single  layer  of  stereo  sense  wires.  Eleven 
of  the  51  wire  layers  are  stereo  layers. 


iS 

Table  3.2 

Material  cocirihution  of  the  various  detector  elements 


Detector  Element 

Amount  of  Material  (%  R.L.) 

Beam  Pipe 

0.40 

B.P./PTL  interface 

0.06 

PTL 

0.29 

PTL/VD  interface 

0.63 

\'D 

0.04 

VD/DR  interface 

1.09 

DR 

0.42 

Tbe  stereo  wires  are  positioned  at  various  small  angles  relative  to  tbe  axial 
direction.  Tbe  sign  and  magnitude  of  the  angle  vary  as  a function  of  radii. 
Tbe  wires  In  the  first  stereo  layer  are  positioned  at  3.77°,  the  second  U at 
-4.22“  the  last  two  at  -6.45°  and  6.89°.  respectively  (30).  All  sense  wires  are 
20  am  diameter  gold-plated  tungsten,  the  field  wires  are  made  of  gold-plated 
aluminum  (layers  1-40)  and  gold-plated  copper-beryllium  layers  (41-51). 

The  outer  and  inner  volume  are  defined  by  segment  cathodes  surfaces. 
These,  as  in  the  VD,  form  the  field  cage  for  the  inner  and  outer  surface  wire 
layera.  The  segmented  anodes,  together  with  the  stereo  sense  wires  provide  for 
polar  angle  measurements.  A schematic  of  the  DR  is  shown  in  Figure  3.5. 

Momeptum  and  angular  resoliulon.  The  momentum  resolution  cau  be  pa- 
rameteriaed  by  two  components.  One  component,  due  to  distortion  of  the 
track's  true  helix,  is  caused  by  multiple  scattering  as  the  particles  traverses 
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Figure  3.6  A view  of  a quadrant  of  the  CLEO  II  detector. 


the  various  detector  elements.  A substantial  portion  of  the  scattering  mate- 
rial. 1.08%  radiatiou  lengths,  lies  at  the  interface  between  the  intermediate  and 
outer  drift  chambers.  To  minimise  the  effect  of  the  concentration  of  material, 
a kink  is  allowed  in  crack  fitting  at  the  radius  of  this  interface.  A second  com- 
ponent is  due  to  the  error  in  measuring  the  track  curvature  due  to  iadividual 
measurement  errors  in  the  drift  distance.  Both  of  these  romponents  contribute 
to  the  transverse  momentum  resoiiitiou  as  [30j 


where  s (in  meters)  is  the  accuracy  of  the  individual  position  measurements, 
Pt  (in  Ge\  ) is  the  transverse  momeutum,  B is  tlie  magnetic  held  in  tesla.  L 
is  the  length  (in  meters  over  which  measuremenU  ore  made),  n is  number  of 
position  measurements  and  ( is  the  thickness  of  material  in  radiatiou  lengths. 


resohitioii  for  a particle  with  sufficiept  momcnluuj 


The  expected  piomcnium 
to  traverse  the  entire  tracking  chamber  is 

(^)  + 

with  S = 1.5  T,  a = 150  pm.  n = dO,  L = 0,85  ni,  and  t = 0.D25  r.l..  The 
scattering  material,  in  radiation  lengths,  of  the  detector  elements  are  listed  in 
Thble  3,2.  An  expanded  side  view  of  the  CLEO  II  detector  is  illustrated  in 
figure  3.6shnwing  the  relative  positions  of  the  various  detector  elements. 

The  angular  resolution  of  the  combined  detector  elements  is  obtained  by 
studying  a sample  of  e+e"  -»  p+p-  events.  These  events  provide  an  estimate 
for  the  resolution  at  high  momentum.  The  rms  measured  resolutions  are 

5©  = 1 mrad;  dff  = 4 mrad 

wliere  8©  and  SB  are  the  resolutions  in  azimuth  and  polar  angles,  respectively. 
The  larger  resolution  in  polar  angle  is  due  to  the  smaller  number,  a maximum 
of  IS  measurements,  and  diminished  accuracy  of  polar  measuremenis  made. 


Particle  identification.  The  identity  of  charged  tracks  are  obtained  through 
specific  energy  loss  mcasnrements  performed  by  the  51  layer  outer  drift  cham- 
ber. The  energy  loss  at  each  hit  point  is  determined  by  measuring  the  charge 
collected  on  the  sense  wire.  The  raw  measuremenM  are  then  corrected  for  dip 
angle,  drift  distance,  entrance  angle  and  proximity  to  the  stereo  layers.  A 
five-dimensiouaJ  1600  bln  mapping  is  calibrated  and  then  mied  to  correct  the 
charge  associated  with  each  hit.  For  tracks  associated  with  40  or  more  good 
him.  resolutions  of  6.2%  for  Bhabhas  and  T.1%  for  minimum  ionizing  pions 


Figure  3.7  Energy  loss  (d£/iir)  as  a funrtion  of  particle  momenta, 

The  top  figure  shows  the  two  dimensional  scatter  plot  of  real  data.  For  clarity 
the  bottom  figure  depicts  the  centroid  and  la  bauds  for  charged  pions  and 


Figure  3.8  Barrel  time-of*flight  rounter. 

are  obtaiued.  The  dEjdx  system  provides  /ftr  separation,  at  the  2<r  level,  for 
tracks  of  momeatuin  up  to  700  MeV. 

The  specific  ionization  energy  losses  arc  a function  of  the  velocity  of  the 
particle.  By  plotting  the  dE/dx  as  function  of  niotueutum  the  particles  sepa- 
rate into  distinct  bands.  Figure  3.7  shows  the  particle  identification  capabilities 
of  the  CLFO  II  detector.  Present  In  the  plot  are  the  bands  formed  by  elec- 
trons, charged  pions  and  kaons,  protons  and  deuterons.  The  latter  are  formed 
predominantly  by  beam  wall  coilisions. 


The  time-of-fiight  (TOF)  system  is  used  primarily  as  a fast  primary  trigger 
in  data  recording  and  as  a tool  for  particle  identification.  Particle  identification 
is  obtaiued  by  measuring  the  time  it  takes  a particle  to  reach  the  TOF  scin- 
tillation counter  referenced  to  CESR  timing.  The  TOF  system  consists  of  two 
major  components,  the  barrel  and  endcap  time-of-flight  systems.  The  barrel 
counters  cover  the  polar  angle  region  from  36°  to  144°.  The  endcap  counters 
cover  the  region  from  15°  to  36°  and  144°  to  165°.  The  solid  angle  subtended 
by  the  combination  of  barrel  and  endcap  counters  is  97%  of  4s. 
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Figure  3.9  Endcap  timMf-fliglit  counteri. 

Barrel  Iime-of-6iebt  counters.  Figure  3.8  is  a drawing  of  a single  barrel 
time-of-flight  counter.  There  ate  s total  of  64  barrel  counters  strapped  to 
the  inside  surface  of  the  crystal  calorinieter.  The  counter  assembly  consist 
of  a organic  scintillation  material  (Bicron  BC-408)  optically  connected  to  a 
UVT  Incite  light  pipe.  The  scintillation  material  is  5 cm  thick.  The  size  is 
chosen  to  maximize  the  thickness  without  degrading  the  performance  of  the 
Cs!  cabrimeter.  The  scinriUatbn  signals  are  cnllected  by  photomultipliers 
glued  to  both  ends  of  the  light  pipe;  each  barrel  TOF  counter  is  thus  readout 
by  two  photomultipliers. 
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Endtip  time-of-flieht  rniintgrs.  A drawing  of  tho  eodcap  iime-o(-flight 
counlors  is  shown  iii  Figuie  3.9.  The  figurr  shows  the  eudeap  TOF  counter 
configuration  and  the  relative  location  of  the  barrel  counters.  The  numbers 
on  the  barrel  counters  represent  the  phototube  addresses  for  the  east-end  pho- 
totubes. Since  each  barrel  counters  is  readout  out  by  two  phototubes,  the 
west-end  phototube  addresses  are  even  numbered. 

There  are  28  wedge-shaped  counters  mounted  in  a circle  on  each  endcap 
calorimeter.  The  thickness  and  scintillation  material  are  identical  to  the  bar- 
rel counters  but  the  geometrical  design  differs  since  the  endcap  counters  are 
designed  to  record  particle  as  they  travel  through  the  ends  of  the  detector. 
The  narrow  end  of  the  scintillators  are  shaped  into  a 45'  prism  upon  which  is 
attached  a photomultipEer  tube.  This  configuration  reduces  the  need  lot  bng 
light  guides  and  additional  glue  joint  which  would  degrade  the  endcap  TOF 
performance. 


L To  obtain  timing  resolutions 
in  the  time-of-fiigbt  system,  a sample  of  Bhabhas  events  was  used  in  the  cal- 
ibration procedure.  The  measured  time-of-flighls  are  calculated  6om  events 
for  which  there  is  only  one  hit  per  counter  and  the  momentum  and  trajectory 
of  the  electron  are  known.  The  timing  measurement  is  then  calculated  from 
phototube  TDC  counts  and  electronic  calibration  constants.  The  phototube 
measurement  is  then  compared  to  the  expected  lime,  T,j.  The  expected 
time  is  parameterizes  by  various  quantities  which  relate  the  physical  properties 
of  the  counters  and  detector  system.  They  include  the  time  require  for  an  elec- 
tron to  move  from  the  interaction  point  (IP)  to  the  scintillator,  the  distance 


Momentum  (GeV/c) 

Figure  3.10  Barrel  tiuie-of-aigbt  \ jB  versus  trank  momcDta. 

between  the  point  where  eieetroa  crosses  the  scintillator  to  the  end  of  the  scin- 
tillator. the  pulse  height  and  shape  of  the  signal,  the  offset  constant  for  each 
chanuel  and  the  velocity  of  the  signal  propagation-  Three  parameters,  the  time 
of  flight  from  the  IP,  the  velocity  of  signal  propagation  and  the  pulse  shapes 
are  determined  by  minimizing  the  difference  between  the  measured  and  ex- 
pected time-of-flight  with  the  Bbabha  data.  AUowing  each  phototube  to  have 
independent  velocity  parameters  results  iu  a resolution  of  139  ps  for  Bhabhas 
in  the  barrel  TOP.  For  pious  from  hadronic  events,  with  momentum  greater 
than  700  MeV,  the  resolution  obtained  was  154  ps  rms.  The  timing  resolution 


provides  a 2a  separation  for  pions  and  kaons  at  a inomenta  of  I.l  GeV.  The 
separation  of  hadrons  as  a function  of  momenta  is  shown  in  Figure  3.10.  The 
plot  derives  from  liadronic  events  at  the  T(45). 

Witha  similar  calibration  procedure,  the  resolution  of  the  endcap  thne-of- 
flighl  counters  ore  also  determined  from  Bhabhas  events.  The  timing  resolution 
in  tbe  endcaps  is  272  ps.  Caiibrating  tbe  endcaps  counters  proves  to  be  a more 
difficult  task  since  advantage  of  tbe  position  dependent  correiations  cannot  be 
taken  in  counters  readout  by  only  one  phototube.  Also,  tbe  use  of  Bhabhas 
compUcates  tbe  calibrations  procedure  due  to  the  spread  in  timing  from  the 
conversion  of  secondary  electrons  in  the  drift  chamber  cudplates. 

3.2.4  Electromaenetic  Calnrlmeter 

The  electromagnetic  calorimeter  consists  of  7,800  thaUiuin-doped  cesium 
iodide  (Csl)  crystals,  and  like  the  TOF  system,  is  divided  into  barrel  and 
endcap  portions.  The  large  number  of  ciystals  provides  for  an  angular  seg- 
mentation finer,  by  an  order  of  magnitude,  than  previous  generations  of  of 
crystal  detectors  at  e*e"  fadlitios  (30). 


95%  of  the  solid  angle.  The  barrel  coverage  starts  at  a polar  angle  ol  32“. 
overlapping  with  the  endcap  which  ends  at  36“.  Each  barrel  crystal  is  shaped 
into  a tapered  rectangular  block  approximately  5 cm  square  by  30  cm  (16  r.l.) 
long.  There  are  0144  crystals  in  the  barrel  portion  of  the  calorimeter.  These 
are  arranged  into  a nearly  vertex-pointing  geometry  of  48  rovrs  along  the  s 
direction  with  each  row  segment  containing  128  block.s  along  the  azimuth.  The 
;-rows  ate  located  symmetrically 
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Figure  3-11  Barrel  crystaJ  calorimeter  holder. 

about  the  axial  direction  and  have  identical  Upered  shapes.  The  shapes 
and  gaps  between  blocks  were  selected  such  that  photons  originating  at  the 
interaction  re^on  strike  the  barrel  at  nearly  nonnal  incidence.  The  crystals 
arc  held  in  place  by  the  crystal  holder  shown  inFiguie3.il.  This  strurturc  not 
only  bears  the  considerable  weight  of  the  entire  crystal  calorimeter  assembly, 
approximately27x  11)3  leg,  hut  also  provides  a clean  and  dust  free  environment 
for  each  individual  ervstaj  bbek. 
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Figure  3.12  Single  quadrant  endcap  crystal  geometry. 

The  enjcgp  cryalal  calorimeier.  Each  endcap  portion  of  the  cTystalcalorime- 

ter  consists  of  828  Csl  rectangular  blacks  5 cm  square  by  3D  cm  long.  These 
are  arranged  as  shown  in  Figure  3.12.  The  pattern  of  207  blocks  is  repeated 
every  90'  forming  a four-fold  symmetric  pattern.  The  aluminum  holder  has  an 
inner  radius  of  0.321  and  0.914  meters,  respectively,  and  is  capped  by  circuloi 
cover  plates.  Each  front  plate  is  located  1.248  m from  the  interaction  point. 

The  scintillation  light  from  the  Csl  crystals  is  converted  into  electrical  sig- 
nals by  four  sUicon  photodiodes  mounted  on  6 nim  thick  CVT  Incite  window 


on  the  back  of  each  ciysta]  block.  Eadi  photodiode  is  connected  to  an  in- 
dependent nearby  preamplifier.  The  four  preamplifiers  are  connected  to  the 
mixer/shaper  card,  wliich  sums  the  shapes  and  signals  from  each  crystal  for 
input  to  the  ADC  circuits. 

Shaa’er  reconstnirtion.  To  reconstruct  showers,  crystals  must  be  grouped 
together  to  form  clusters  of  hits  which  define  a single  shower.  The  cluster 
finding  algorithm  first  locates  the  most  energetic  crystal  in  a cluster  that  ex- 
ceeds 10  MeV  aud  has  an  energy  higher  than  any  of  its  immediate  neighbors. 
Members  of  the  same  cluster  can  not  be  more  than  two  bbcks  away  from  each 

The  energy  aud  iiosition  of  each  cluster  is  then  computed  by  taking  the  N 
most  energetic  crystals  iu  the  cluster.  The  value  of  N varirrs  logarithmically  as 
a function  of  energy  from  4,  at  25  MeV,  to  17,  at  4 GeV,  These  values  of  N 
are  chosen  to  minimize  the  energy  resolution.  This  algorithm  proves  to  be  the 
most  effective  possible  [30].  If  the  crystal  belongs  to  two  separate  cluster,  its 
energy  is  apportioned  between  the  two. 

The  position  vector  of  each  cluster  is  computed  by  first  finding  the  centroid 
of  the  shower.  The  centroid  is  calculated  as  the  energy-weighted  sum  of  the 
coordinates  of  each  crystal's  geometric  center.  The  vector  is  then  corrected  for 
shower  location  in  the  plane  perpendicular  to  the  incident  particle's  direction 
and  depth  within  the  cesium  iodide  where  the  an  electromagnetic  shower  would 
have  deposiled  its  mean  energy.  The  location  of  the  shower  within  the  crystal 
is  important  when  matchiug  showers  to  charged  tracks  aud  in  assigning  correct 
anguJar  pasitions  wheo  showers  fall  in  the  endcap. 
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are  calculated  using  Bhablia  events.  The  constants  are  computed  by  minimie- 
ing  the  rms  widths  of  the  Bbabha  electron  shower  energy  distrihutiou  and 
constraining  it  to  peak  at  the  beam  energy.  To  cover  tlie  entire  pliotou  energy 
range,  various  techniques  are  used  to  calibrate  the  absolute  energy  normaliza- 
tion for  photon  dusters.  They  include  w®  -n.  radiative  Bhabhas  and  two 
photon  events.  Photon  pairs  from  ir®  decays,  mass  constraint  imposed,  yield 
a calibration  accurate  to  about  ±0.5%  below  2 GeV.  To  sample  higher  energy 
photons,  radiative  Bhabhas  are  used,  again  resulting  in  calibration  accuracies 
of  ±0.5%. 

With  these  calibrations,  tbe  photon  energy  resolution  in  the  barrel  region  is 
1.5%  at  5.0  GeV.  In  theeudeap  region  the  energy  resointion  degrades  to  2.6%. 
At  5.0  GeV  the  anguJar  resolution,  in  the  azimuthal  diiectba,  is  3 mrad  in  the 
barrel  and  D mrad  in  tbe  endcap  regions.  The  angular  and  energy  resolutions 
can  parameterized  in  the  barrel  region  by  the  expression 


<r^(miad)  = -^  + 2.5  irslmrad)  = O.Sir^sin^ 
where  the  photon  energy  £ is  in  GeV.  In  the  cndcap  region  the  parameteriza- 


(T*(mrad)  = ^ + 7.3  ff*(mrad)  = ^ + 5.6. 

These  parameterizatioiis  come  from  Monte  Carlo  simulations  of  showers  in  the 
CLEO  II  detector.  Not  included  in  the  simulation  are  effects  duo  to  other 
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tracks  in  the  event.  The  effects  var>'  depending  on  the  event  topology  and 
tend  to  slightly  degrade  the  resolution. 

3.2.5  The  CLEO  11  Maenet 

The  CLEO  II  magnet  consists  of  a 3.5  m long  supercouducting  cylindrical 
coll  with  a bore  of  2dl  meters  and  a large  steel  return  yoke.  The  coll  is  cooled  by 
liquid  helium  contained  in  a 700  1 dewar  mounted  on  top  of  the  detector.  The 
helium  is  delivered  through  uaiural  flow  circulation  know  as  a thermosyphon 
system.  The  system  is  self  regulating  and  avoids  the  need  for  liquid  helium 
pumps  which  Insures  that  refrigeration  malfunctions  will  not  result  in  a magnet 
quench. 

The  return  yoke  has  three  layers  of  steel  interspersed  with  muon  detection 
equipment  and  thus  serves  not  only  as  a medium  lor  the  magnetic  flux  return 
but  as  a hadronic  absorber.  The  entire  assembly  weighs  about  8 x 10*  kg,  most 
of  it  in  the  slabs  of  steel  arranged  in  on  octagonal  geometry  around  the  coil. 

The  conductors  used  in  the  coil  windings  arc  flat  high  purity  aluminum 
coextruded  over  Cu-NbTi  wires.  The  flat  wires  measure  5 x 16  mm  and  were 
wound  in  two  layers,  one  above  the  other.  Both  layers  are  mounted  to  the 
inner  surface  of  the  coil  shell  with  insulation  material  and  a radiation  screen 
separating  the  coil  assembly  from  the  vacuum  vessel  surfaces.  The  wires  carry 
a current  of  3300  A which  provide  the  1.5  tesla  field  used  to  bend  the  paths 
of  charged  particles.  The  magnetic  field  is  uniform  to  0-.2%  throughout  the 
drift  chamber  volume.  During  data  taking  tlic  magnetir  field  is  monitored 
by  NMR  probes  located  a tew  centimeters  beyond  the  end  of  drift  chamber. 
The  effective  average  magnetic  field  is  calibrated  by  using  event  reconstruction 
involving  p pairs  aud  known  masses  of  particles  such  as  D\  ds  and  A‘°s.  This 


Figure  3.13  Schematic  of  muon  detector  superiayer. 


yields  a correciion  of  O.U%  at  the  ceoter  of  the  magnet  comparable  to  the 
overall  accuracy  of  field  measurements. 

3.2.1!  The  Muon  rthamhers 

The  muon  detector  [31)  is  the  last  of  the  detectors  component  which  will 
be  described.  It  covers  the  polar  angle  range  from  30“  to  150“,  about  85%  of 
the  solid  angle.  The  barrel  muon  chambers  are  embedded  in  the  iron  yoke  of 
the  magnet  at  depths  of  36,  72  and  108  cm  (see  Figures  3.2  and  3.3).  The  two 
ends  of  the  detector  are  also  covered  by  muon  chambers.  Depending  of  the 
direction  of  the  track,  the  total  equivalent  thickness  of  iron  absorber  a particle 
sees  varies  from  7.2  to  10  radiation  lengths. 

The  muon  chambers  are  streamer  counters  operating  in  the  proportional 
mode,  at  2500  V.  Within  each  iron  gap  there  are  three  layers  of  counters,  the 
middle  layer  suggered  with  respect  to  the  other  tvro.  ^ schematic  of  the  three 
layers  (known  as  a superiayer)  is  shown  in  Figure  3.13.  Each  counter  is  about  5 
m long  by  8.3  cm  wide  and  is  divided  into  eight  cells  by  graphite  coated  plastic 
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«-alls  wluch  define  the  electric:  field  cage.  Within  each  cell  lies  SO  ym  diameter 
sliver-plated  Cu-Be  anode  wire.  The  eight  anodes  arc  gauged  together  giving  a 
spatial  resolution  of  2.4  cm,  better  than  the  spatial  uncertainty  due  to  multiple 
scattering  for  panicles  reaching  the  muon  chambers.  Measurements  along  the 
c direction  ore  provided  by  external  copper  strips  positioned  along  the  counter 
and  by  charge  division  readout  of  the  anodes. 

3.3  The  Trigger  Decisinr  l.evele 


figuration  of  CESR  consisted  of  seven  electrons  positrons  beams  each  consisting 
of  seven  bunches  per  beam.  The  eollisiou  rate  between  bunches  in  this  con- 
figuration was  3C4  usee  for  all  but  the  first  and  last  bunch  where  the  interval 
was  increased  to  378  nsec.  The  beam  crossing  rate  thus  averages  about  2.73 

MHz.  Considering  the  large  number  of  channels  involved  and  the  non-oegiigible 

amount  of  time  required  to  decide  whether  an  particular  event  warrants  a write 
to  tape,  reading-out  the  detector  at  this  rate  would  prove  to  be  extremely  dif- 
ficult. Also,  interesting  physics  does  not  occur  at  every  beam  crossing.  In  fact, 
given  the  current  CESR  luminosity  and  the  value  of  the  hadronic  CTTWs-scction 
suggests  a more  reasonable  rate  of  a few  tens  of  Hz.  The  CLEO  II  trigger 
system  [32]  was  thus  designed  to  take  the  2.73  MHz  bunch  crossing  time  and 
read  out  the  detector  only  at  the  occurrence  of  “interesting  events". 

3.3.1  IVigger  Txigir 

To  accomplish  these  goals  a three  tiered  trigger  system  was  devised.  Tlie 
different  decision  levels  are  referred  to  as  LO.  U and  L2.  Each  successive 
trigger  level  requires  that  the  previous  level  be  satisfied.  In  the  abscuce  of 


Figure  3.14  Flow  of  the  tiered  trigger  in  CLEO  II. 

triggers,  the  detector  is  actively  collecting  time  and  pulse  height  inforraaiioii 
and  gated  every  beam  crossing.  Satisfying  the  LO  trigger  disables  the  gate  ant] 
a search  for  all  LI  combinations  defined  is  initiated.  If  all  triggers  are  satisfied 
the  detector  is  read  out  and  the  gate  is  enabled.  The  trigger  logic  is  sketched  in 
Figure  3.1J.  The  various  elements  which  feed  information  into  the  trigger  logic 
circuits  arc  depicted  on  the  left-hand  side,  included  are  the  TSP.  BLT  and  PD 
logic  algorithms  which  combine  various  detector  element  to  do  rudimentary 
track  finding. 
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The  lewl  0 (LO)  trigger  Is  very  fast,  very  simple  and  quite  discriminatory. 
It  reduces  the  rate  from  the  2.73  MHz  crossing  frequency  to  20  kHz.  It  receives 
input  from  the  time-of-flight  scintillators,  from  the  vertex  detector  and  from 
the  Csl  crystal  calorimeter.  The  experimenter  can  define  a set  of  criteria  tor 
LO,  all  of  which  are  logically  OR'd.  The  uumber  of  ciossinp  “lost"  when  a 
valid  LO  combination  fails  to  find  a valid  LI  combination  is  currently  six. 

The  LI  level  trigger  takes  the  20  kHz  rate  from  LO  and  reduces  it  to  25 
Hz.  It  takes  its  input  from  scintillators,  vertex  detector,  central  drift  chamber 
and  electromagnetic  calorimeter.  The  time  it  takes  for  all  of  these  devices  to 
be  ready  for  LI  integratioD  is  1.0  psec.  With  an  LO  rate  of  20  kHz  this  implies 
a dead-time  of  2%.  Each  set  of  Ll  criteria  has  its  own  set  of  L2  requirements 
which  are  examined  once  the  particular  Ll  trigger  is  met.  Currently  the  L2 
trigger  only  takes  input  from  the  vertex  detector  and  drift  chamber.  Additional 
readiness  time  from  L2  only  adds  about  0-25%  to  the  overall  detector  dead- 


3-3-2  TVieger  Innms 

Briefly  outlined  in  the  following  section  are  the  logic  elements  that  con- 
tribute to  the  trigger  system.  A more  detailed  description  of  the  various  inputs 
can  be  found  in  Ref.  [32). 

TimT-of-fl'gtl'  trigger  logic.  The  time-of-0ight  system  leads  iuelf  well  to 
the  trigger  system  due  to  its  fast  response  time.  Inputs  to  the  trigger  sys- 
tem from  the  TOP  scintillators  consi-sls  of  ganging  four  adjacent  counters  into 
sectors  which  reduces  the  granularity  of  the  detector.  This  coarsening  effec- 
tively reduces  both  the  lime  and  number  of  bits  required  to  readout  the  data. 


Six  flags  are  set  involving  various  combination  of  endcap  and  barrel  sector 
hits.  These  can  be  modifled  by  the  experimenter  to  select  for  different  physics 


CrvMal  calorimeter  trigger  loeic.  The  inputs  from  the  crystal  calorimeter 
(CC)  follows  a similar  potteru  hs  with  the  TOP  system.  The  barrel  segmenta- 
tion consists  of  384  segments  of  16  crystals  per  segment.  The  828  crystals  in 
the  endcap  are  also  grouped  to  form  segments  but  here  the  number  of  crystal 
pet  segment  vary  ffoui  6 to  16.  The  signals  from  the  grouped  segments  are 
summed  and  Che  result  is  discriminated  by  both  a high  and  low  level  com- 
parator. The  high  being  approximately  twice  the  minimum  ionizing  wliile  the 
low  is  half  minimum  ionizing.  While  all  segments  are  readout  for  potential 
software  triggers  the  LO  and  LI  hardware  triggers  group  these  segments  iuto 
larger  settore.  In  the  barrel  region  the  entire  detector  is  divided  into  16  sectors, 
eight  along  the  azimuth  and  two  aloug  polar  direction.  The  endcap  segments 
are  grouped  into  eight  sectors  each.  All  segments  within  o sector  are  logically 
OR'd.  The  choice  of  16-sector  topology  is  made  to  match  the  address  space  of 
commercially  available  fast  sutic  random  access  memory  (RAM). 

Drift  rliamber  trieeer  loeic.  The  drift  chamber  and  vertex  detector  are 
used  as  input  devices  to  the  trigger  system.  Signals  from  these  devices  are 
incorporated  into  the  trigger  logic  by  using  several  algorithms.  These  perform 
crude  track  finding  and  also  allow  the  experimenter  to  specify  the  charge  par- 
ticle multiplicities  and  momenta  necessary  for  the  acceptance  of  events.  The 
Tfack  Segment  Processor  (TSP)  uses  analog  signals  generated  by  charge  col- 
iected  on  the  anode  wires  in  all  layers  of  the  vertex  detector  and  12  layers  of 


the  drift  chamber  at  various  radii  from  0.21  to  0.87  meters.  The  TSP  algo- 
rithm groups  the  anodes  from  tire  10  'V'D  and  12  DR  layers  into  a total  200 
biu  of  information,  These  serve  as  inputs  to  32  PROMs  (programmable  reail- 
only  memory);  each  memory  ha.s  a.s  its  address  10  bits,  one  from  each  layer 
aligned  In  azimuth.  To  find  a track,  each  PROM  compares  the  bit  pattern 
of  the  event  with  the  radial  pattern  expected  from  charged  particles.  Using 
a pattern  recognition  scheme,  the  TSP  searches  for  and  counts  the  number  of 
short,  medium  and  long  tracks  in  the  event.  The  definition  of  what  coitstitutes 
a long,  medium  or  short  track  defined  at  run-time.  For  long  tracks,  nominally 


and  time-of-flight  counters  are  used  in  LI. 

A second  processor  used  in  the  LI  and  L2  trigger  decision  level  is  the  Binary 
Link  Thacker  (BLT)  (33j.  This  processor  groups  the  wires  in  the  DR  into  gangs 
of  wires  to  form  segments  consisting  of  17  layers  at  various  radii.  The  algorithm 
uses  correlations  between  the  bit  pattern  set  by  the  various  segments  to  find 
trad:  momenta  and  multiplicity  of  the  event.  Combining  the  TSP  and  BLT 
processors  increases  the  LI  efficiency  by  S-10%  over  using  the  TSP  alone.  The 
BLT  is  also  employed  in  the  L2  decision  level  where  more  time  is  allowed  for 
calculations.  At  this  stage  the  BLT  maybe  used  to  establish  the  sign  of  a track 
thus  enabling  tbe  rejection  of  beam  wall  events  which  give  a preponderance  nf 
positive  tracks. 

A third  processor,  the  Precision  IVackiiig  Device  (PD),  uses  the  hit  infor- 
mation exclusively  from  the  vertex  detector-  By  examining  hit  patterns  which 


are  consistent  with  tracks  originntiiig  at  the  vertex,  it  helps  in  rejecting  cos- 
mic ray  and  beam  related  events.  The  algorithm  also  discriminates  against  a 
large  number  of  hits  clustered  io  a given  azimuth,  thins  further  biasing  against 
beam-wall  and  beam-gas  events.  The  PD  is  used  in  the  level  2 decision. 


3.4  Data  Acquisition 

A new  data  acquisition  systems  was  installed  in  the  Fall  of  1991  to  readout 
the  CLEO  II  detector.  The  change  was  required  to  take  advantage  of  the 
steadily  increasing  CESR  luminosity  and  also  to  allow  for  the  broadening  of  the 
search  for  interesting  physics.  Since  the  previous  data  acquisition  system  was 
used  for  less  than  a quarter  of  the  data  sample  and  has  been  described  elsewhere 
[34]  we  will  only  present  a brief  description  of  the  new  data  acquisition  system, 
known  as  the  DAQ90  system. 

The  DAQ90  system  [35)  is  responsible  for  digitizing  and  collecting  data 
from  the  40  front-end  detector  creates  reading  out  the  various  detector  ele- 
ments. It  also  sparsifies  the  data  to  remove  unwanted  information,  formats 
the  fragmenus  into  complete  events  and  distributes  them  to  the  online  work- 
station. The  performance  goals  of  the  system  is  to  handle  a 50  Hz  trigger  rate 
with  a dead  time  of  no  more  than  10%.  This  corresponds  to  a bandwidth  of 
250  kBytes/s,  for  an  average  sparsified  event  size  of  5 kBytes. 

The  D.AQ90  system  hardware  is  illustrated  In  Figure  3-15.  The  process 
of  writing  data  to  tape  begins  at  the  froot-end  crates  where  the  pulse  heights, 
timing  information  and  position  bits  are  readout  and  digitized  from  the  various 
detector  elements.  As  they  ore  digitized,  the  data  from  each  front-end  crate  is 
sent  via  a point  to  point  connection  to  a remote  dual-port  RAM,  link  board, 
where  they  are  buffered  for  subsequent  transfer.  After  the  transfer  of  data  is 
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Figure  3.15  A block  dingram  of  the  DAQ90  hardware. 

complete,  each  crate  controller  sends  a “done’  signal  to  its  dedicated  link.  The 
data  arc  then  transfered  from  the  link-board  memory  to  ITvlE  memory  once 
all  link-board  “done"  signals  have  been  received  from  its  respective  branch. 
The  four  \'ME  CPU  boards  further  process  and  sparsily  the  event,  in  parallel, 
while  communicating  with  the  task  supervisor  which  sends  signals  to  the  trigger 
system  and  examines  error  logs.  Once  the  sparcihcation  of  all  event  fragments 
is  complete  thesparsifier  task  (SPA)  signals  the  event  builder  task  (EVBD)  that 
all  event  fragments  are  available.  The  EV'BD  then  gets  the  data  from  buffers  on 
the  sparsifier  nodes  and  merges  the  fragments  into  a ZEBRA-formatted  event. 
The  event  consumer  requests  events  frnm  the  event  distributor  and  allows  for 
the  application  of  software  filters,  in  the  currently  configuration  a single  event 


Figure  3-ia  Au  illustration  of  the  iafoiroalion  flow  in  Mie  DAQSO  system. 

builder  feeds  a single  event  distributor,  a DECstntion  5000/240,  which  also 
runs  the  LVL3  software  filter  and  delivers  the  data  to  another  DECstatiou 
5000  for  on-line  processing.  The  data  flow,  described  above,  is  illustrated  in 
Figure  3.16.  The  illustration  shows  how  event  fragments  propagate  through 
the  senes  of  buffers  to  arrive  at  the  event  distribntor  as  a whole  formatted, 
sparsified  event. 

3.5  The  CLEO  II  .Software 

As  with  any  High  Energy  experiment,  the  CLEO  II  software  is  comprised 
of  many  layers  of  utility  programs,  online  and  offline  analysis  routines  and  a 


multitudp  of  user  apfdications.  It  would  tfquire  at  iaasl  au  additional  200 
pages  just  to  begin  to  describe  all  of  the  software  compoDeots  used  to  com- 


plete an  analysis.  Instead,  iu  the  concluding  sections  we  will  concentrate  on 
describing  the  two  main  track  fitting  processors  used  in  the  online  and  offline 
reconstruction  programs,  sketch  the  procedure  used  to  generate  the  Monte 
Carlo  simulations  of  events  and  finajly,  outline  the  high  level  script  lauguage 
(CABS)  u-sed  to  fully  reconstruct  the  decay  modes  considered  in  this  auaiysis- 

3-5-1  PASSl:  The  Ouliue  Ansivek 

The  on-line  analysis  program-  PASSl,  runs  on  a DECstation  5000  con- 
nected through  FDDI  to  the  DAQOO  data  acquisition  system.  It  is  the  first 
pass  at  transforming  the  raw  hii  iuformntion.  recorded  at  runtime,  to  fully  re- 
constructed events.  It  is  also  used  to  monitor  the  operation  of  the  experiment 
and  to  create  the  first  set  of  calibration  constants  later  used  in  the  off-line 
analysis  program,  PASS2.  PASSl  also  provides  event  filtering,  rejeciuig  beam- 
wall,  beam-gas  and  cosmic  ray  events  and  classifies  the  remaming  events  before 
storage  outo  4 mm  DAT  tapes. 


PASS2  is  the  off-line  analysis  program  it  also  transforming  raw  hit  infor- 
mation into  fully  reconstructed  evenU  consisting  of  charged  and  neutral  tracks. 
The  algorithms  used  are  more  sophbticated  than  those  used  iu  PASSl  and  thus 
require  more  time.  The  main  tracking  algorithms  used  in  PASS2  are  TRIO 
and  DbET  [36].  These  use  the  raw  hit  information  and  find  tracks  which  fit  a 
series  of  wire  hits  in  a given  event.  TRIO  is  a fast  system  which  lakes  advan- 
tage of  the  triplet  geometry  of  tlie  grouping  of  layers  and  the  staggering  of  the 


trackiug  chamber  cells  to  fiod  tracks.  The  serood  program  uses  a tree  finding 
algorithm  that  provides  for  more  accurate  track  fitting  and  higher  efficiency  in 
find  tracks.  The  TRIO  algorithm,  being  much  faster,  is  used  in  PASSl,  Both 
DUET  and  TRIO  are  used  in  PASS2. 

l&JO-  The  TRIO  program  searches  for  groups  of  three  hits  in  consecutive 
layers  forming  units  called  tracklets.  TRIO  t hen  forms  combinations  of  inner 
and  outer  tracklets  to  aeate  candidate  tracks.  Additional  tracklets  are  included 
if  they  fall  witliiu  a given  distance  and  lie  Itetween  the  initial  two  tracklets. 
The  average  number  of  tracklets  in  a candidate  track  is  four.  The  candidate 
track  is  then  fit,  and  if  of  sufficient  quality,  added  to  the  track  iist- 

PVET-  The  DUET  program  starts  off  with  seed  tracks  from  TRIO.  It 
looks  for  two  hits  segments  within  a narrow  region  of  the  candidate  TRIO 
track.  Tlie  segments,  known  as  links,  are  then  combined  into  larger  and  larger 
segments  called  chains.  Candidate  chains  are  fit  to  a circle  and  a straight  line 
m the  r - © plane.  Chains  which  satisfy  the  requirement  are  refit,  in  three 
dimensions,  to  a helix.  Once  DUET  exhausts  all  TRIO  seed  tracks,  it  begins 
its  own  search  for  tracks  from  hits  unused  in  the  initial  search, 

QtJier  processor  in  PASS2.  Both  PASSl  and  PASS2  mu  under  the  CLEVER 
[37)  (CLEO  EVEut  Reconstruction)  program.  It  provides  a fully  interactive 
user  interface  to  mauipnlate  ZEBRA-formatted  data.  The  user  can  select  from 
among  a series  of  analysis  routines  called  PROCESSORS.  In  addition  to  TRIO 
and  DUET,  other  proreasors  are  included  in  the  PASS2  executable.  These  in- 
clude the  three  processors,  CCFC,  CDCC  and  XBAL,  which  provide  informa- 


tion  on  the  raw  crystal  hits  that  aJloa’s  for  the  formation  of  shower  candidates. 
The  CCFC  routines  find  the  position  and  angles  of  individual  crystals  {see  Sec- 
tion 3.2.4)  and  the  CDCC  provides  the  routines  to  match  showers  to  charged 
tracks.  The  XBAL  routines  provide  an  alternative  way  to  access  the  crystal 
information.  They  arc  based  on  the  clustering  algorithms  developed  for  the 
Crystal  Ball  experiment.  Other  processors  do  vertex  finding  (VFND).  event 
filtering  (LVL3  and  BMWL),  particle  identification  (TFAN  and  DEDR),  and 
calculation  of  global  event  properties  (GLOB).  In  addition,  a processor  called 
MAKE  is  used  to  compress  the  large,  .fax,  ZEBRA-formatted  files  produced  by 
CLEVER.  It  makes  a DST  (Data  Summary  “Tape")  in  ROAR  format  achiev- 
ing lager  compression  factors  by  minimising  the  storage  by  excluding  raw  hit 
information.  The  ROAR  data  format  can  then  be  unpacked  by  the  analysis 
programs  into  FORTRAN  common  blocks. 

3-5.3  Event  .SimulatioTis 

To  understand  the  recon-slruction  efficiencies  of  the  CLEO  II  detector  and 
the  eBects  of  the  various  selectiou  requirements  imposed  in  the  reconstruction 
algorithm,  we  use  Monte  Carlo  simulations.  The  Monte  Carlo  is  calibrated 
against  world  averages  of  previously  well  measured  results.  In  this  analysis 
the  Monte  Carlo  simulations  are  used  to  find  the  detector  acceptance  neces- 
sary in  determining  absolute  branching  fractions.  However,  Using  the  Monte 
Carlo  simulations  to  find  these  acceptances  introduces  systematic  tracking  er- 
rors which  depends  on  the  multiplicity  and  hadronic  content  of  the  event.  The 
tracking  error  introduces  a 2%  (595)  error  per  charged  track  above  (below)  250 
MeV  and  5.095  error  per  neutral  pion.  The  tracking  errors  are  also  determined 
by  calibrating  the  Monte  Carlo  to  weU  measured  processes. 
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The  generation  of  Monte  Carlo  events  begins  with  the  event  geiterator  QQ, 
[38|.  The  program  generates  the  initial  four-vectors,  event  vertices  and  the 
general  topology  of  the  event.  It  was  orlgiaaily  written  to  generate  qq  . BB  and 
ggg  jets  and  has  evolved  into  a shell  structure  which  can  call  a variety  of  other 
event  generators  such  KORALB  to  generate  tan  pairs  and  LUND/JETSET 
applied  in  hadronising  qq  and  ggg  events.  Once  the  events  are  created  a lookup 
table  with  information  on  all  known  decay  inodes  and  branching  fractions  is 
used  to  decay  the  generated  particles.  A user  con  change  the  decay  information 
of  any  particle.  The  angular  distributions  of  all  decays  are  also  included  in  the 
lookup  tables.  These  can  also  be  specified  or  changed  by  the  user. 

The  second  step  in  the  Monte  Carlo  chain  is  the  detector  simulation.  The 
event  history  from  QQ  are  saved  in  ZEBRA  format  and  used  as  input  to  the 
CLEOG  [39]  program.  CLEOG  is  the  CLEO  II  interface  to  the  GEANT  [40] 
detector  simulation  routines.  The  job  of  CLEOG  is  to  propagate  each  gener- 
ated particle  through  the  various  detector  element  keeping  track  of  multiple 
scattering,  interactions  with  detector  material,  electron  conversions  and  prop- 
erly decaying  long  lived  particles  such  the  pions  and  kaons  when  needed.  The 
response  of  each  detector  element  is  recored  just  as  in  a real  event,  digitized 
and  written  to  banks  in  the  same  manner  as  the  data.  These  data  are  then  fed 
to  the  P.ASS2  routines  used  in  track  finding  and  fitting. 

For  this  analysis  millions  of  simulated  events  were  used.  These  include 
two  million  events  employed  in  the  study  of  BB  backgrounds  and  hundreds 
of  thousands  of  events  used  in  the  study  of  detector  acceptances  and  angular 
distributious  of  the  various  B decay  modes  investigated. 


The  large  number  of  modes  considered  in  this  analysis,  60  in  all,  would  tra- 
ditionally require  the  user  to  write  many  thousands  of  lines  of  FORTRAN  code 
to  impletnenl  fully  all  the  reconstructed  events.  Nobu  Kalayama  has  written 
a simple  scripting  language  for  CLEO  II  called  CABS  141).  The  language  re- 
duces the  number  of  lines  a programmer  writes  to  a few  lines  of  easily  readable 
code.  The  small  number  of  Hues  also  make  it  possible  to  quickly  and  reliably 
implement  changes  in  reconstructing  and  analyiing  data. 

CABS  takes  advantage  of  an  object  oriented  approach  to  programming.  All 
particle  definitions  are  considered  as  objects  which  are  endowed  with  quantities, 
classes,  tagged  to  the  particular  instance  of  the  partide-decay  definition.  The 
tagged  quantities  are  internally  defined  kinematic  variables  or  cap  defined  by 
the  user.  Furthermore,  objects  can  be  used  to  construct  other  objects  while 
retaining  all  previously  defined  properties.  Another  feature  of  CABS  is  the 
elimination  of  user  written  DO  loops  since  this  is  handled  automatically  by 
the  high-level  design  phiiosophy.  Al  that  is  requited  of  the  user  is  to  define  the 
properties  of  the  charged  tracks  and  showers  and  combines  the  objects  with 
other  objects  to  form  composite  particles. 


CHAPTER  -1 

ANALYSIS  PROCEDVRE 
liilrodiirlini) 

The  objective  of  this  analysis  U to  measure  various  properties  of  B mesoos 
by  exclusive  reconstruction  of  their  decays.  To  do  so.  we  first  defined  the 
selection  criteria  for  each  of  the  final  decay  products,  kaons,  pions  and  pho- 
toDS,  combined  these  into  candidates,  selected  from  among  these  by  Imposing 
additional  requiremenls,  defined  further  candidates  until  we  obtained  a My 
reconstructed  B meson  candidate.  Further  selection  criteria  are  then  imposed 
OD  the  B candidate  and  their  properties  are  examined  by  measuring  event 
yields  in  distributions  of  particular  interest.  In  this  chapter,  we  will  describe 
the  data  sample  used,  the  set  of  global  event  requirements  imposed  and  the 
selection  criteria  and  methods  ttsed  in  reconstructing  the  B mesons.  Also,  the 
performance  of  the  selection  criteria  will  be  discussed  in  some  detail,  with  em- 
phasis on  Hie  efficiency  and  amount  of  background  removed  1^  imposing  the 
various  selection  requiremeDts. 

4.1  Data  Sample 

The  data  used  In  this  analysis  were  collected  by  the  CLEO  II  detector  at  the 
CESR  e‘*'e“  storage  ring  beginning  in  November  of  1900  and  ending  in  January 
of  199.3.  During  this  period  various  changes  were  made  to  the  CESR  accelerator 
and  CLEO  II  detector.  These  included  improvements  to  data  acquisition  bj- 
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the  instaUaCiou  0/  the  DAQ90  aystem  and  LVL3  software  triggers  during  the 
Summer  of  1991,  replacing  the  gas  in  the  PT  with  di-methyl  ether  iu  April  of 
1992  and  replating  the  synchrntrou  RF  cavities  with  new  MARK  III  cavities 
in  the  Fall  of  1992.  These  changes  affected  detector  acceptance  and  resulted 
in  an  increase  to  the  total  number  of  events  in  oiir  data  sample.  The  changes 
did  not  significantly  affect  the  quality  of  the  data  and  did  not  require  separate 
analysis  procedure  for  different  data  sew. 

The  data  sample  used  consisted  of  a total  integrated  luminosity  of  2.04 
on  the  T(4a)  resonance  and  0.97  /fi"'  off  resonance.  The  entire  sample  was 
divide  into  seven  data  sets,  labeled  the  4s2  through  4s8,  each  corresponds  to 
several  months  of  data-taking.  The  luminosity  times  BB  cross-section  for  this 
sample  gives  a total  of  2,2  x 10«  BB  events.  The  off  resonance  sample  was 
taken  approximately  26  MeV  below  the  T(45)  and  was  used  for  continuum 
background  studies. 

The  dau  processed  by  the  off-line  event  reconstruction  programs  were  di- 
vided into  various  data  type  eJassificatious.  These  ineJuded;  QED  (Bhabhos. 
p pairs  and  radiative  Bhabhas),  e^e"  ->  t*t~,  two  photon,  hadronic  and 
various  classes  of  junk  events,  i.e.,  beam-wall  and  beam-gas  interactions.  The 
QED  sample  was  used  primarily  for  detector  calibration  since  these  processes 
are  well  understood  theoretically.  The  hadronic  sample  used  was  extracted 
from  the  full  data  set  by  imposing  a set  of  selection  criteria  designed  to  retain 
hadronic  events  wMle  eliminating  QED,  junk  end  other  non  BB  events. 

4.2  Global  Event  Requirements 

The  hadronic  sample  was  dehned  by  events  which  passed  the  following 
requiremenw:  the  number  of  charged  tracks  in  the  event  greater  than  or  equal 


to  four,  the  total  euergj-  dapositcd  io  the  talorimeter  less  than  65%  of  the  ceuter 
of  mass  energy  (Ecm  ) or  the  highest  energy  shower  must  have  an  energy  less 
than  60%  of  Ecm  , the  visible  energy  (computed  from  showers  and  chatged 
tracks)  greater  than  20%  of  E^n,  and  the  vertex  position  of  the  event  was 
tequired  to  be  within  2 cm  (5  cm)  in  the  direction  perpendicular  (parallel) 
to  the  beam  with  respect  to  the  nominal  run  averaged  colliston  point.  These 
requirements  elimiunted  more  than  99%  of  all  Bhabhas  while  retaining  more 
than  99%  of  the  BB  and  95%.  of  generic  udsc  and  c?  hadronic  continuum 

The  topological  characteristic  or  shape  of  the  event  was  also  used  in  the 
selection  CTiteria.  We  imposed  two  global  event  shape  cuts  to  help  separate 
continuum  from  BB  events,  the  variable  R^.  or  the  ratio  of  the  2"^  to  O"" 
Fox* Wolfram  moments,  and  the  event  sphericity  angle.  These  two  cuts  are 
correlated  and  together  eliminated  about  35%  of  the  (ontinuum  while  retaining 
96%  of  the  BB  events. 

The  R2  |42]  variable  is  defined  as  the  ratio  R2  = H2IH0  where 

= (4,1) 

The  indices  ij  run  over  alJ  tracks  in  the  event,  is  the  angle  between  the 
tracks!  andj,  f)  is  the  l"‘  Legendre  polynomial  aud  £is  the  total  energy  of  the 
event.  The  mimetator  in  the  ratio  H2/R0  is  proportional  to  J(3cos^(^j)  - 1) 
while  the  denominator  is  used  to  normalize  Hj.  A value  close  to  I is  obtained 
when  the  event  is  jet-like.  For  events  that  are  distributed  isotropically  the  value 
tends  too.  Since  we  were  interested  in  BB  decays,  which  decay  isotropically,  we 
placed  a cut  on  R%  which  accepted  evenuwith  a value  of  0.5  or  less.  In  Figure 


Figure  4.1  Rj  plotted  for  continuum  and  BB  Monte  Carlo. 

4.1  a plot  of  for  continuum  and  B6  eventt  is  shown.  The  requirement  on  ^2 
alone  is  seen  to  eliminates  a substantial  amount  of  continuum  while  retaining 
most  of  the  BB  events. 

A second  global  event  shape  requirement  made  was  on  the  cosine  of  the 
sphericity  angle.  The  sphericity  angle  {&/)  is  defined  as  the  angle  between  the 
major  sphericity  axis  calculated  using  all  tracks  consistent  with  a B candidate 
and  the  major  sphericity  axis  calculated  using  the  remaining  tracks  in  the 
event.  The  major  sphericity  axis,  or  jet  axis,  is  calculated  by  diagonalizing  the 
momentum  tensor  Ty  defined  as 
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Figure  4.2  Diagram  of  the  sphericity  axis  superimposed  on  typical  events. 

= (4-2) 

where  Pi(n)  and  Pj(n)  are  the  and  component  of  the  momentum  vector 
and  the  sum  runs  over  all  applicable  tracks.  The  major  axis  corresponds  to 
the  smallest  eigenvalue. 

A pictorial  description  of  the  sphericity  angle  is  shown  in  Figure  4.2  with 
the  major  sphericity  axis  superimposed  on  a typical  continuum  aad  a low 
multiplicity  BS  event.  The  figure  illustrates  the  tendeaev  of  continuum  events 
to  duster  nearcoslSy)  = ±1,  while  all  values  of  cos(9y)  are  equally  likely  for 
BB  events  since  the  angle  &/  can  assume  any  value  from  0 to  In.  The  major 
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COS(0.) 


Figure  4.3  The  cos(6/)  distribution  in  a continuum  and  a BS  sample. 

sphericity  axes  are  indicated  by  dashed  lines,  tlic  heavy  solid  lines  are  the 
charged  tracks  which  combine  to  form  a B meson  and  the  light  soUd  lines  are 
the  tracks  that  make  up  the  rest  of  the  event. 

In  this  analysis  a cut  on  the  |cos(9/)|  was  imposed  on  each  event  after 
it  pa-ssed  all  B candidate  requirements.  In  modes  with  a prominent  signal 
the  cut  was  determined  by  maximizing  the  signal  to  noise.  In  modes  where 
a signal  was  not  evident  the  value  of  the  cut  used  was  taken  from  a mode 
with  similar  kinematics  and  multiplicity.  The  cut  thus  varied  mode  by  mode 
with  the  tightest  restrictiou,  |cos(e/)|  < 0.7,  imposed  on  B decays  to  a a, 
plus  a charmed  meson  and  the  loosest  cut,  |cos(6/)l  < 0.9,  on  B decays  to 
a D'  plus  a single  pion.  In  Figure  4,3  a plot  of  the  cos(e;)  distribution  for 
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Figure  4-4  Diggram  of  signed  impact  parameter  with  respect  to  beam-spot. 

continuum  and  BB  Monte  Carlo  events  illustrates  the  discrimination  power  of 
this  variable.  The  plot  shows  that  a cut  of  lco5(6/)|  < 0.8  eliminates  80%  of 
continuum  events  and  retains  80%  of  the  B6  events. 

4-3  Charged  ftack  Selectlnn 

The  varions  requirements  imposed  on  charged  track  candidates  not  only 
provided  good  quality  tracks  originating  at  the  primary  vertex  but  also  helped 
in  determining  whether  the  particular  track  wa.s  a kaon  or  a pion.  The  latter 
was  accomplished  by  using  the  d£/dx  particle  ID  system  described  in  Chapter 
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momentum  (GeV) 


Figure  4.5  DBCD  plotted  vs.  momcDCa  for  4s2-4s8  and  Monte  Carlo  data. 

To  Insure  that  tracks  originated  from  the  primary  vertex  a cut  on  the 
momentum  dependent  impact  parameter  with  respect  to  beam-spot  (DBCD) 
was  imposed  (The  diagram  in  Figure  4.4  defines  the  variable  DBCD).  The  cut 
used  accepted  tracks  with  an  impact  parameter  less  than  40i„i,p  for  track.s  with 
momentum  above  250  MeV.  fbr  tracks  with  momentum  less  than  250  MeV  tbe 
requirement  was  changed  to  The  quantity  is  tbe  spread  in  the 

measured  impact  parameter  with  respect  to  beam  spot  measured  for  a large 
sample  of  charged  tracks,  The  spread  depends  on  the  momentum  of  the  tracks 
as  weU  as  the  angle  d they  make  with  the  j-  axis.  The  momentum  dependence 


of  the  spread,  for  each  of  the  data  sets  and  a Monte  Carlo  sample,  is  shown  in 
Figure  4,5,  A detailed  analysis  Is  presented  in  Appendix  A. 

The  iEjdx  particle  identification  system  was  used  to  establish  whether  a 
particular  track  was  a kaou  or  a pion.  First,  tracks  were  required  to  have  more 
than  ten  central  drift  chamber  hits  to  Insure  that  energy  loss  information  was 
available.  The  A'/tr  candidate  was  then  selected  if  its  energy  loss  measurement 
differed  from  that  expected  for  a A'/tt  hypothesis  by  less  than  3ir,*  Since  low 
momentum  tracks  do  not  travel  far  into  the  central  drift  chamber  all  particle  ID 
requirements  were  suspended  if  the  track's  momentum  was  less  than  250  MeV, 
These  tracks  were  only  used  as  slow  pion  candidates  in  D"*  reconstruction. 


Photons  were  selected  from  clusters  of  cells  in  the  barrel  region  of  the  crys- 
tal calorimeter  with  an  energj’  greater  than  30  MeV,  Because  of  the  degradation 
in  energy  resolution  in  other  regions  of  the  detector,  the  energy  requirement 
for  the  cluster  in  the  endcaps  was  increased  to  50  MeV  and  dusters  which  fell 
in  the  endcap-bairel  overlap  region,  between  32'  and  36“,  were  rejected  alto- 
gether, If  the  daster  matched  a charged  track's  projection  into  the  calorimeter 
the  duster  was  rejected  as  a photon  candidate.  To  form  a match,  the  track’s 
projection  was  required  to  be  within  8 cm  of  at  least  one  ceil  in  the  duster. 

To  distinguish  electromagnetic  showers  tom  hadronic  showers  we  used  the 
shower’s  lateral  shape  distribucion.  Electromagnetic  showers  deposit  their  en- 
ergy in  an  narrow  region  around  the  central  crystal  and  tend  not  to  spread  as 
hadronic  showers  do.  A quantity  Eg/Ejs  is  defined,  for  each  candidate  duster, 

• cr  is  the  rms  resolution  of  the  expected  dElix  energy  loss  for  the  A'/ir 


hypothesis. 


as  the  suni  of  energy  of  the  block  of  9 crystals  divided  by  the  block  of  25  oys- 
tals  both  centered  on  tbe  central  crystal.  Tbe  quantity  depends  on  the  energy 
and  polar  angle  of  tbe  cluster  and  is  close  to  one  for  electromagnetic  showers. 
This  quantity  was  evaluated  for  each  candidate  shower  and  compared  with  the 
value  expected  for  99%  of  all  electromagnetic  showers.  Since  candidate  shower 
were  used  only  to  reconstruct  neutral  plom  or  etas,  we  used  a slightly  modified 
vereion  of  E0/E25.  referred  to  as  E9/E25  unfolded.  This  quantity  allows  for 
the  possibility  of  merged  sliowers  in  high  momentum  two  photou  candidate  by 
using  not  only  the  crystals  assigned  to  the  shower  but  also  a fraction  of  the 
shared  crystals. 


4-5  Particle  RcconstnirtinTi 

We  now  turn  to  describing  the  criteria  used  to  reconstruct  tbe  rest  of 
tbe  particles  which  will  ultimately  be  used  in  reconstructing  the  B meson 

candidates.  The  recoustruction  method  used  is  rather  simple:  a chosen  particle 

candidate  is  made  by  combining  the  4-momenla  of  its  decay  products.  A set 
of  quantities  are  then  calculated  for  the  composite  particle,  constraints  are 
Imposed,  and  the  selected  particle  candidate  is  stored  for  use  in  reconstructing 
the  B or  other  particles. 

4.5.1  Light  Neutral  and  CharppH  Mesnits 

Neutral  pion  candidates  were  first  formed  by  combining  two  photons  and 
reqmrmg  that  the  invariant  mass  of  the  combination  be  within  2.5o  of  the 
fitted  neutral  pion  mass.  The  neutral  pion  mass  distributions  are  asymmetric 
and  tbe  asymmetry,  the  resolution  of  the  asymmetric  Gaussian  and  the  fitted 


Figure  4.6  Two  photon  invariant  ma&s  distributions  at  different  momenta. 


The  invariant  mass  for  the 
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ir”  candidates  with  momenta  between  (a)  500  and 
and  1000  MeV,  (c)  between  I.O  and  1.5  and  (d) 


lible  4.] 
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centra]  value  vary  slightly  as  a function  of  momenta,  (the  momentum  depen- 
dence of  these  quantities  are  shown  in  Figure  4.6  where  the  invariant  mass  is 
plotted  ID  four  different  momentum  bins).  To  compensate  for  this  we  varied 
the  mean,  n,  and  ir„  as  a function  of  momenta  as  shown  in  Table  4.1.  Candi- 
dates with  momentum  less  than  250  MeV  vrere  rejected  if  one  or  both  of  their 
photon*  came  from  the  endcap  or  overlap  region.  Finally,  the  neutral  pions 
were  kinematically  fit  with  their  masses  constrained  to  the  known  pion  mass 
and  a cut  of  12.0  was  imposed. 

Etas  decay  to  two  photons  about  39%  of  the  time.  The  two  remaining 
channels  with  relatively  large  branching  fraction,  the  3^°  and  Tr+ir'a®,  modes 
were  not  used.  The  Sir®  mode  contained  too  many  poaible  2-f  combinations 
to  be  useful,  while  the  ir+tr-ff®  mode  provided  little  additional  information 
in  determining  0^7  or  8“  -»  D’°r,  branching  fractions  due  to  large 

background  levels. 

The  I)  -)•  77  candidates  were  selected  from  two  photon  combinations  that 
had  at  least  one  of  the  photons  in  the  barrel  region.  The  candidates  were 
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Decay  modes  and  resolutions  for  Ihe  D and  D‘  candidates 


Deray  Mode 
D“-.K-7t+ 

D° 

D<'  -»  A'-ir+jr-7T+ 


Resoultion  (MeV) 
<rmoo  = 90 
i7mo«  = 13.0 
<Tmo»  = e.i 


D*  K-ir-'-tr+ 
D'“  -*  Dltr® 
D-+  D®7r+ 


required  to  have  an  imariant  mass  within  30  MeV  of  the  PDG  value  (2)-  Can- 
didates with  asymmetric  decay,  i.e.,  |cos«^|  > 0-85  were  rejected.  The  decay 
angled^  is  defined  as  the  angle  betvreen  the  direction  of  one  of  the  photons  and 
the  direction  of  the  eta  in  the  lab  frame  as  measured  in  the  rest  frame  of  the 
eta  Real  etas  have  a fiat  cosd^  distribution  tvhile  random  27  combinations, 
that  happen  to  fall  within  the  eta  mass  cut,  are  strongly  peaked  in  the  region 
near  cos9,j  = 1.  This  cut  eliminates  background  from  fake  etas  but  only  re- 
dtices  the  sample  of  true  etas  by  15%.  Once  candidates  passed  all  requirements 
they  were  kinematically  fit  with  their  masses  constrained  to  the  PDG  q mass 
value  and  a < 12  was  impesed. 

The  rest  of  the  light  mesons  vrere  reconstructed  in  the  following  decay 


The  mass  was  constrained  to  be  within  one  full  decay  width 
(ISO  MeV)  of  Ibe  nominal  p mass. 


required  to  be  within  30  MeV  of  tlie  nominal 


The  mass  was 


The  mass  of  each  candidate  i/  was  required  to  be  within  30 
Me\'  of  the  nominal  i)'  mass.  The  j)'  mode  was  ex- 

amined but  yielded  no  additional  iafonnation  in  determining 
fl  branching  fractions  was  obtained  due  to  increase  in  back- 
grounds without  providing  a signal  in  tlie  event  j-ield.  Only 
the  r/  -*  i|jr+7r'  mode  with  i?  -*  r>  was  used. 

The  0)  meson  has  a very  large  and  poorly  determined  decay 
width.  Its  mass  was  constrained  to  lie  between  1-0  GeV  and 
1.6  GeV.  The  mass  of  the  p was  constrained  in  all  O]  candi- 
dates as  described- 

4.5.2  Charmed  Meson  SelecHnn 

Unlike  the  light  mesons,  pseudoscalar  charmed  mesons  can  decay  through 
a large  number  of  channels  each  with  a relatively  small  branching  fractions. 
A large  number  of  these  modes  have  multiple  neutral  pions  and  or  neutral 
kaons,  which  have  low  reconstruction  efficiencies.  The  losses  ineuned  when 
fully  reconstructing  a 0 makes  most  of  these  modes  unavailable  hence  only 
three  modes  were  used  in  reconstructing  ifis  and  one  mode  for  the  D+s.*  The 
relevant  mode.s  used  are  listed  in  Table  4.2  together  with  their  measured  rms 
' Charge  conjugate  states  are  implied  if  not  indicated  explicitly. 
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Figure  4.7  iDvariant  mass  distributions  of  A* mr  combinations. 

In  (a)  the  invariant  mass  distribution  of  the  A'“ir+  combination,  in  (b)  the 
^ (n)  the  A'”ir*7r”ir‘‘'  combination  and  in  (d)  the 

K*ir~ir  combination. 


101 


/no— m^-6m  (OeV) 


Figure  4.8  mass  dUTereace  plot. 


mass  resolutions,  Cmo  1^3].  The  invariant  mass  distributions  for  the  four  decay 
modes  used  are  shown  in  Figure  4.7.  Each  plot  cousist  of  a subsample  of  the 
dataset,  the  4s6,  and  is  plotted  for  candidates  whose  momenta  lie  between  2.0 
GeV  and  2.5  GeV.  The  arrows  in  the  figure  show  the  bounds  of  the  2.5omu 
mass  requirement  imposed  ou  all  candidates. 

Ebr  B modes  which  decay  to  a D\  we  took  advantage  of  the  excellent  res- 
olution on  the  mfl.  - mp  mass  difference.  In  D'  decays  the  small  amount  of 
available  phase  space  leads  to  a sharp,  narrow  peak  in  the  tup.  -mp  distri- 
bution. The  sharp  peak  was  exploited  in  reconstructing  the  B ->  D’  modes 
by  eliminating  non  D*  combinations  that  were  outside  the  mp.  -mp  mass 
region.  An  example  of  the  mp.  — mp  ma^  resolution  is  shown  in  Figure  4.8 


with  the  2.5(T  cut  indicated  by  the  arrows.  The  D°  wa-s  rctoastructed  in  the 
D®  -►  decay  mode. 

For  D'  candidates  we  first  required  the  roaiis  o!  tlie  to  be  within  2.5tr  of 
iM  measured  resolution  and  then  imposed  the  mo-  - mp  requirement.  As  with 
pscudoscalar  decays  we  only  considered  charmed  wetor  meson  decay  modes 
which  had  both  larse  branching  fractions  and  low  backgrounds.  The  decay 
modes  and  their  measured  rms  resolutions  ate  also  listed  in  Table  4.2  [43j. 

4.6  B Mesons  Srlectinn 

All  the  necessary  ingredients  are  now  iu  place  to  f 
mesons.  The  uon  color-suppressed  decay  modes  used  w 


B"  -t  D“fl- 
fl-  D‘°p- 


the  color-suppressed  modes  wt 


B'-tD-V 

(4-4) 

For  eacli  mode  various  physical  properties  of  the  fully  reconstructed  B mesons 
were  analyzed  to  determine  brauching  fraction-s  aud  polarization  states.  The 
results  and  couelusions  will  be  discussed  in  Chapter  6. 

To  calculate  the  branching  fraction  we  need  to  determine  the  number  of 
fully  reconstructed  B mesons  in  our  sample.  To  do  so  we  used  the  beam 
constrained  mass  (Mgc  ) defined  as 


Mbc  = , I Eita 


(4.5) 
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where  p,  are  the  reconstructed  momenta  of  the  decay  products  of  the  B and 
^beam  the  energ)’  of  the  e'*'  e~  beam.  The  Mac  distribution  has  a resolution 
an  order  of  magnitude  better  than  the  invariant  B mass  distribution  and  was 
used  to  determine  event  yields.  The  resolution  of  Mils  variable  is  dominated  by 
the  spread  in  the  beam  energy  which  Is  approxiuialely  3 MeV  at  the  T(4S)  . 
The  small  width  gives  a sharply  peaked  Mac  distribution,  at  the  B mass  and 
was  thus  easy  to  fit  with  a signal  plus  backpouad  shape. 

The  criteria  used  to  seleet  the  B meson  candidates  varied  with  the  decay 
mode  used.  Two  quantities,  liowever,  were  used  in  selecting  all  candidates,  the 
B meson's  polar  angle  9g  and  the  difference  between  its  reconstructed  energy 
and  the  energy  of  the  beam, 

- Erccanitruacd-  (4,6) 

Since  other  selection  aiteria  are  mode  dependent  we  defer  their  discussion  to 
later  sections. 

The  cut  imposed  on  the  polar  angle  was  |cos(fiB)|  < 0.95,  it  takes  advan- 
tage of  the  spatial  distribution  resulting  from  angular  momenta  conservation 
in  the  decay  of  T(45).  The  T(4S)  is  produced  in  c"*"  e~  collisions  through  a 
vutual  photon  with  J = 1,  Af  = ±1.  Its  decay  to  two  spin  zero  objects  is  de- 
termined by  angular  momentum  conservation  forring  the  orbital  distribution 
to  follows  a sin^lSg)  form. 

The  AE  distribution  is  centered  on  zero  with  a width  determined 
by  the  energy  resolution  of  the  individual  particles  in  the  final  state.  Tl;e  AE 
widths  were  measured  for  each  exclusive  mode  and  found  to  vary  from  14  to 
44  Me\‘.  These  measurements  were  performed  on  tagged  Monte  Carlo  data 
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Figure  4.9  distributions  of  tbe  S" -*  Z>“ir“  modes. 

In  (a)  the  D®  ->  fc'"jr+  mode,  in  (b)  the  D®  -i  A'"ir+7r®  mode  and  in  (c)  tlie 
D°  -•  K“ir‘*'7r”ir‘*'  mode. 


which  was  found  to  be  consistent  with  measurements  ijerfbrmed  on  the  data 
sample.  In  Figure  4.9  the  iE  distributions  for  all  three  B~  -*  D^tr"  modes  ate 
shown.  Each  iE  distribution  was  fit  to  a Gaussian  plus  a straight  line.  These 
measured  AE  resolutions  were  then  compared  with  fits  to  Monte  Carlo  data 
and  were  found  to  differ  by  no  more  than  where  cr,iE  ^ ^he  statistical 

error  of  the  resolution  measurements. 

A mode  dependent  AE  cut  of  2.5o^e  w“  imposed  on  each  candidate  to 
discriminate  against  background  events  resulting  from  poorly  reconstnicted  B 
and  or  other  B decays  which  differ  by  one  or  more  pion.'  In  B decays  where  a 
p is  produced  the  resolution  of  AE  is  affected  by  the  momentum  of  the  neutral 
pion  produced  from  the  decay  of  the  p.  This  dependence  is  parameterized  bv 
the  helidty  angle  0/,  (defined  below)  and  the  cut  on  AE  varied  as  a function 
of  0/|  for  these  modes. 

The  largo  number  of  particles  in  each  event  sometimes  leads  to  multiple  B 
candidates  per  mode  that  pass  all  the  selection  criteria.  In  these  instances,  the 

* This  is  discussed  in  greater  detail  in  Chapter  5. 
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Table  4.3 

Selection  criteria  for  non  color-suppressed  B -»  decays 


6 Channe. 

D Sub-channe: 

|coa(e,) 

liEI  (GeV) 

D°ir- 

K-;r+ 

K-ir+ffO 

<0.8 

<0.055 
< 0.065 

D+rr- 

A'-ir*T+ 

<0.8 

<0.051 

K-n* 

K-n*n° 

A-)T+7r-!r+ 

<0.9 

< 0.050 

< 0.068 
< 0.047 

D*  + 7T- 

K-ir+ir'' 

K-T+JT-JT+ 

<0.9 

< 0.050 

< 0.062 
<0.045 

Table  4.4 

Selection  criteria  for  non  folor-suppressed  B decays 


B Chatme! 

D Sub-chaunel 

|cos(0,)| 

|AE|  (GeV) 

D°p- 

K-jr* 

K-x+7r-ir+ 

<0.8 

< (0.025  •cos(ea)-b  0.070) 

< (0.025  cosOftl-P  0.079) 

< (0.025  •cos(eft)-t- 0.067) 

D+p- 

K-X+X+ 

<0.8 

< (0.025  • cosOa)  + 0.0701 

4D‘“p- 

A-x+7r» 

A-x*x-x+ 

.... 

< (0.025  •cos(ea)-f  0.078) 

< (0-025  •cos(eft)-P  0-089) 

< (0.028 -cos(0,,) -1-0.072) 

D-^p- 

<0.8 

<(0.025  C05(eft)-t- 0.079) 

< (0.025  cc6(eA)-i- 0.083) 

< (0.025  •cos(ej)  + 0.076) 
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candidate  with  the  "best"  smallest)  AE  value  was  chosen-  This  procedure 
of  pidtins  the  candidate  with  the  best  AE  picks  the  correct  candidate  86.2%  of 
the  time.  The  success  rate  was  detenniiied  b>'  comparing  lagged  and  untagged 
Monte  Carlo  yields  in  the  -*  D“7r“  with  D°  -*  A-x+  . This  procedure 
also  tends  to  bias  the  sample  of  selected  B to  those  with  AE  values  close  to 
zero.  However,  the  results  are  unaffected  since  the  Monte  Carlo  simnlatious 
correctly  predict  the  loss  iu  efficiency  for  this  choice  of  the  "best"  B candidate. 

The  cuts  used  vary  depending  not  only  on  the  B decay  channel  but  also  the 
D subchannel  considered.  In  order  to  Est  the  differences  in  a systematic  way 
we  divide  tlie  modes  into  color-suppressed  and  non  color-suppressed  modes  and 
then  further  divide  the  non  coior-suppresaed  modes  into  modes  containing  one. 
two  and  three  pion  light  mesons.  The  various  AE  and  global  sphericity  angle 
cuts  are  discussed  iu  the  following  suteections. 


4.6.1  Non  Color-Sunorp-Ked  Mndos 

B ->  modes.  In  B modes  that  decay  to  a charmed  meson  plus  a 
single  pioD  we  used  the  cuts  described  above  with  no  additional  requirements 
imposed.  The  values  of  the  sphericity  angle  and  AE  cuts  are  listed  in  Table 
4.3.  The  values  listed  in  the  table  correspond  to  2.5  times  the  measured 
resolutions. 

P ippdes.  In  these  modes,  we  employed  a cut  designed  to  take 

advantage  of  the  spatial  distribution  of  the  p's  decay  products.  The  spatial 
distribution  is  best  parameterized  by  a variable  which  we  refer  to  as  the  helicity 
angle  6* . In  p"  -»  x®x-  decays,  it  is  defined  as  tlie  angle  bctvreen  the  direction 


Figure  4.10  Diagrsim  of  helirity  angle  definitions. 

The  helitity  angle  9/,  (A  = p.ui)  is  shown  (or  the  decay  modes  (a)  B~  -> 
and  (b)  -*  . 


of  thc:r®,  in  the  o“’srest  frame,  and  the  direction  of  the  p~  in  the  fl  candidate’s 
rest  frame  (see  Figure  4.10).  As  an  example,  we  take  the  decay  B"  -» 

The  p-  is  produced  longitudinally  polarised  along  its  direction  of  motion  in  the 
0 's  rest  frame.  Simple  angular  momeutum  conservation  shows  that  the  decay 
of  a |1.0)  object  to  two  spinless  pious  is  governed  by  a cos^(©a)  distribution. 
For  B decays  to  a o plus  another  vector  particle,  a D',  the  effectiveness  of  a cut 
based  on  these  arguments  is  somewhat  compromised  since  now  the  initial  total 
angular  momentum  state  of  the  p~  is  undetermined.  We  could  require  that 
both  the  D‘  and  the  p have  identical  helicity  distributions,  but  this  reduces 
the  acceptance  without  much  improvement  in  the  signal  to  noise  in  the  already 
clean  modes  where  the  D'  - D mass  difference  cut  cau  be  used. 

The  helicily  angle  cuts  used  depended  not  only  on  the  B decay  channel 
but  also  on  the  D subchannel  considered.  In  the  decay  B~  ->  a cut  of 

l'=05(S/.)l  >0.4  vms  used  tor  all  three  subchannels.  For  the  D” -»  K" jr+jr“ 
mode  a large  amount  of  background  was  eliminated  by  discarding  0 candi- 
dates with  a soft  ffO  produced  in  the  p~  decay.  This  was  accomplished  by 
excluding  candidates  with  neutral  pions  in  the  backward  region  of  the  p decay 
plane  (plane  A in  Figure  4.1D  (a))  thus  a cut  cos(efc)  > 0.4  was  imposed  as 
a consistency  check.  For  the  B®  -*  D*p-  channel  a |cos(eft)|  > 0.4  cut  was 
also  employed. 

The  AE  and  Bf  angle  cuts  used  for  the  0 ^ modes  are  listed 

in  Table  4.4.  The  helicity  angle  dependence  of  the  AE  resolution  discussed 
earlier  is  included  in  the  AE  cut  column.  The  helicity  angle  dependence  of 
the  AE  resolution  was  found  by  measuring  the  AE  resolutions  at  different 
'“(©a)  values  in  Monte  Carlo.  The  plots  where  then  fit  to  straight  lines.  The 


<i,^(MeV) 
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(b) 
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Figure  4.11  Plots  of  belicity  as  a functiou  of  AE. 

AE  resolution  as  a function  of  the  o~  coafes)  for  the  decay  mode  B~ 
with  (a)  Cfi  -*  K-n*.  (b)  D“  -r  K-jr+Tr®  and  (c)  D®  ^ 
fits  for  all  three  modes  are  shown  in  Figure  (d). 


^D®p- 
TT*.  The 
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Table  4.S 


fitted  cos(e/,)  vs.  distributions  are  shown  in  Figure  4.11  for  the  three  D® 
subchannels  in  B~  -»  decays. 

g_->  modes.  As  in  6 decays  to  D plus  a single  plon.  no  additional 

cuts  were  imposed  ou  modes  with  a nj  in  the  final  state-  The  dominant  decay 
modes  of  the  oj,  the  fm  mode,  does  not  allow  the  use  of  a helicity  type  out 
since  the  initial  angular  momentum  state  of  the  g is  indeterminate  and  thus  no 
useful  spatial  angular  distribution  is  available.  The  broad  oj  resonance  together 
with  the  large  number  of  particles  in  the  final  also  results  in  particularly  large 
combinatoric  background  which  hampered  the  ejctraotion  of  signal  events  yields 
in  the  B~  —*  D^n“  and  S®  — * D'^’nJ'  charmels. 


4-6-2  Colur-SunpresapH 

The  foior-suppressed  modes  were  treated  in  a similar  way  to  the  modes 
discussed  above-  The  difference  mainly  lay  in  our  inabilil.y  to  verify  that  the 
Monte  Carlo  predicted  calues  used  to  select  aud  refine  the  cuts  matched  the 
data.  For  example,  we  used  the  AE  resolutions  obtained  from  Monte  Carlo 
only  for  the  D®  mode  in  ail  color-suppressed  modes.  For  the  two 

remaining  Cfi  decay  modes,  we  took  the  values  obtained  for  D®  -*  I<~n*  mode 
and  scaled  that  by  a factor  obtained  from  a similar  non  color-suppressed  mode 
where  the  AE  resolution  could  be  checked  In  data.  This  Is  done  for  all  but  the 
fi®  -♦  D“ir®  mode  since  the  resolution  of  the  D®  -»  fr— jr+  mode  is  already  44 
MeV. 

Eseudo,scjil»r  light  mesoia,  As  with  the  non  color-suppressed  light  pseu- 
doscalar modes,  the  candidates  in  these  modes  were  selected  without  additional 
cuts.  Only  the  standard  AE  and  sphericity  cuts  were  applied.  The  values  used 
for  each  mode  are  listed  in  Table  4.6.  The  |AE|  cut  entries  in  the  table  reflects 
the  2.5  times  the  AE  resolution. 

Vector  light  mesons.  Four  color-suppressed  modes  contain  a light  vector 
mesou  in  the  B decay  product.  Two  of  these,  the  -*  D’°e°  aad  -f 
are  vector-vector  decays  in  which  complementary  helicity  requirements  on  both 
the  D-°  and  its  vector  partner  were  imposed.  The  heUcity  angle  requirements, 
however,  did  not  lead  to  any  improvement  in  obtaining  a signal  in  the  event 
yield  and  were  tlius  not  used  in  the  final  analysis.  The  remaining  two,  the 
B“  -*  J>6p0  and  fl°  -*  D°w,  have  a spin  0 plus  a spin  one  object  in  the  final 
state  so  advantage  was  taken  of  the  helicity  requirement  described  earlier.  In 


Table  4,6 


Sclectioa  ctiteria  fot  folor  suppressed  B D^’’’P  decays 


B Chaime 

D Sub-chamie 

|cos(S,) 

|AE|  (GeV) 

K'-^+ 

< 0.090 

DO,r° 

K-ir+ff» 

<0.8 

< 0.090 

il‘-jr+jr-:r+ 

< 0.090 

A'-s+ 

< 0.074 

D“r/ 

A—it+t'’ 

<0.S 

< 0.084 

K-:r*<r-s-+ 

< 0.074 

K-ir-^ 

< 0.049 

A-ff+rr® 

<0.8 

< 0.059 

K-ir+iT-7r+ 

<0.049 

K-7T  + 

< 0.090 

D’°ir“ 

<0.9 

< 0.090 

A***;r-ir+ 

<0.090 

< 0.060 

D‘^r) 

ii-ir+s® 

<0.9 

< 0.070 

K->r+e-rr+ 

< 0.060 

A"7r+ 

<0.045 

A-rr+ir® 

<0.9 

< 0.055 

A--s*^-x-^ 

< 0.045 
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Tlible  4.7 

Selection  criteria  for  color  supprra.«ed  B -*  decays 


B Channel 

D Sub-diannel 

|cos(6f)| 

|AE|  (GeV) 

D°pO 

K-jt* 

fv"7r+ir® 

K-jr+tr-s* 

<0.8 

<0.044 

< 0.054 

< 0.044 

K-1T+ 

A"tr+ir° 

fC-7T+x-ir+ 

<0.8 

< 0.050 
<0.060 

< 0.050 

K-j* 

K'x+x® 

ff-x+x-x* 

<0-0 

< 0.041 
<0.051 

< 0.041 

ff-x+ 

K-x+x^ 

A"x’’-x-x+ 

<0.9 

< 0.054 

< 0.064 
<0.054 

fl®  -»  c^p'>,  a cut  on  | cos(S;,)|  > 0.4  was  imposed.  Ip  -*  D®_,,  the  helicity 
angle  definition  is  complicated  because  the  u decays  to  three  spinless  pioiu. 
The  9;,  angle  is  now  defined  as  the  angle  between  the  normal  to  the  w decay 
plane  and  the  u direction  in  the  13  meson’s  rest  frame  (see  Figure  4.10  (h)). 
Angular  momentum  conservation  of  the  u decaying  to  three  spinless  particles 
requires  that  the  norma]  to  the  plane  be  distributed  as  a sin^(0|,).  A cut  on 
|sin(0ft)|  <0.6  was  imposed  on  ail  candidates.  The  cos(6;)  and  iE  ruts  used 
are  listed  in  Table  4.7. 


CHAPTER  5 

BACKGROUND  STUDIES 


Introduction 

In  this  chapter  we  investigate  the  backgrounds  which  contribute  to  the 
Mac  spectra  for  the  decay  modes  considered  in  the  previous  chapter.  The  Mac 
background  contributions  can  be  separated  into  two  components:  continuum 
component  and  BB  component.  The  continnum  background  is  due  to  the  large 
non-resonant  fraction  of  the  hadronic  cross-section,  approximately  75%,  from 
direct  e+e‘  -*  qq  production  imderneath  the  T(4S).  The  BB  background 
component  is  a result  of  misreconstructing  other  BB  decays.  The  relative 
contributions  and  the  overall  amount  of  background  varies  decay  mode  by 
decay  mode  depending  primarily  on  the  multiplicity  of  the  B decay. 

The  background  contributions  to  each  Mgc  distribution  were  not  sub- 
tracted when  establishing  the  event  yield.  Instead  the  backgroimds  were  pa- 
rameterized by  a function  that  fit  tbe  general  characteristics  of  the  overall  back- 
ground shape.  Determination  of  the  yield  by  fitting  Mbc  the  spectrum  to  a 
background  function  plus  a signal  Gaussian  assumes  that  the  background  does 
not  peak  significantly  in  the  siguai  region  and  that  the  background  function 
used  adequately  models  the  true  background.  To  insure  that  these  assump- 
tions are  valid  we  petfonned  a detailed  and  careful  study  of  the  backgrounds, 
as  discussed  below. 
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5.1  Backgrounds  iii  .\on  Cnlor-SuDurpMpd  Mode.^ 

To  simplify  thp  analysis  and  the  discussion  thal  follows  wc  divided  the 
study  into  a study  on  the  non-color  suppressed  modes  and  a study  on  the  color- 
suppressed  modes.  The  non  color-suijprc.ssed  mode  were  examined  in  greater 
detail  since  signals  were  obtained  in  all  mode  investigated.  Determining  the 
event  yield  m modes  with  prominent  signals  requires  an  accurate  measure  of  the 
amouut  of  background  so  that  the  systematic  errors  associated  with  the  back- 
ground  shape  could  be  established  accurately.  In  the  inodes  where  signals  are 
not  observed,  sensitivity  to  background  systematics  do  not  significantly  affect 
the  resulM  since  the  errors  here  are  dominated  by  uncertainties  in  statistics. 

5.1.1  BB  Backerniinds 

Background  contributions  from  other  BS  decays  are  impossible  to  com- 
pletely isolate  in  real  data  since  tbey  arc  products  of  the  same  physical  pro- 
cesses which  produce  signal  events.  To  estimate  this  component  we  used  Moute 
Carlo  simulations,  where  complete  control  can  be  mainiainod  over  the  produc- 
tion and  decay  mechanism.  The  Monte  Carlo  data  set  cousLated  of  a sample  of 
1.95  X 10'=  BB  events  in  which  each  B meson  in  the  event  was  allowed  to  decay 
generically.  The  data  set  was  separated  into  sub-samples  where  a particular  B 
meson  decay  chain  was  removed.  The  removal  of  the  decay  chain  here  implies 
the  removal  of  the  complete  B decay  chain.  This  included  aU  possible  reso- 
nant substructure  of  the  D mesou.  For  example,  in  removing  the  B~  -t  £flir~ 
with  the  D®  ^ K-a+ir'’  decay,  the  D°  ^ A-OfO,  the  A— tt*  and  the  l<-p* 
substructure  decays  were  removed  as  well.  Each  siibsample,  referred  to  as  the 
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“aatiskim  sample,  thas  coutained  all  known  B meson  dways  except  the  decay- 
chain  of  the  particular  mode  whose  background  we  investigated. 

The  B reconstruction  algorithm  described  in  Chapter  4 was  then  applied 
to  each  aatiskim  sample  and  the  Mgc  ploU  were  made-  Figures  5.1,  5.2  and 
5.3  show  the  results  for  non  color-suppressed  modes  with  one,  two  and  three 
pion  light  mesou  states  respectively.  Most  of  these  plots  show  no  significant 
enhancement  in  the  signal  region  which  cannot  be  accommodated  by  the  back- 
ground shape  used. 

For  Bd  backgrooud-s  to  peak  in  the  signal  region  two  conditions  must  be 
met:  the  number  and  types  of  final  state  particles  must  match  that  of  the  B 
decay  mode  considered  and  the  total  reconstructed  energy  must  satisfy  the  AE 
requirement.  In  Figure  5.2  (ImI)  an  excess  is  found  in  the  signal  region  of  the 
Mac  spectra.  The  background  candidate  B decay  wliich  meets  these  criteria 
is  the  B-  -t  D"®*-  decay  mode  where  the  D*°  -*  D“-y.  This  mode  can  easily 
fake  a B -*  D®p“  if  the  (T  is  misreconstructed  by  combining  the  fast  rr" 
from  the  B decay  with  a low  momentum  picked  randomly  from  the  large 
number  of  two  photon  combinations  possible  in  the  event.  The  requirement 
that  the  energy  be  that  of  the  fl  is  met  from  the  toss  of  the  photon  in  the 
decay  This  scenario  was  Investigated  by  taking  a sample  of  fl“  ->  D‘°x- 
Monte  Carlo  decays  with  the  flOy  and  reconstructing  it  as  D°p--  The 
Mbc  distributions  are  shown  in  Figure  5-4.  The  loose  constraint  on  the  p~ 
invariant  mass  can  also  conspire  with  X)*''  -t  VOy  decay  to  produce  the  broad 
enhancement  in  the  B~  Lfip-  signal  region. 
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Figure  S.l  SB  Monte  Carlo  Mbc  Spectrum,  fl -*  Dl'ljr. 

In  (a)  the  B«  -*  p*ir-  iu  the  D*  decay  mode  O'*  K-ir*x*.  In  (b)  through 
(d)  the  B“  ^ D^ir"  in  the  D°  decay  modes  A'-tt+tt''  and  A''7r*ir“)r+ 

rwpectively.  In  (e)  through  (g)  the  B"  -t  D-Oir'  in  the  E|0  decay  mode^ 

A IT  .A— ir  T°andA'  ir‘^>r‘’r+.  In  (h)  through  (j)  the  B“ -t  in  the 

O'*  decay  modes  A'  jr*.  K~:r*tr'^  and  A'~7r+ir“!r‘*’. 


Mac  (GeV) 


Figure  S.2  BB  Monte  Carlo  Mgc  spectrum.  B -* 

In  (a)  thefl°-»p+p-  in  tbeD+  decay  mode  D+  -»  A'-ir+Tr-^.  In  (h)  through 
(d)  the  B~  -*  D°p~  in  the  decay  modts  A'"ir‘*’,  A'”a*7rh  aud  A'“ir‘*'ff~7r+ 
respectiwly-,  In  (e)  through  (g)  the  B'  -f  D”>p~  ia  the  Lfi  decay  mode^ 
A IT  .A  ir  a andA'-a+e  In  (h)  through  (j)  the  B° -t  D'+p' in  the 
D“  decay  modes  K'n*,  aud  A'^ir+tr'a"''. 
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Figure  5.3  flS  Monte  Carlo  Mgc  spKtrum,  B -» 

In  (a)  theB‘>-.D+uf  in  the  D+  decay  mode  D+  -*  K'-tt+it*.  Id  (b)  through 
(d)theB  in  the  decay  modes  A''n+,  A—7r+ir<' and  K'ir‘‘'!r-ir+, 

respectively.  In  (e)  through  (g)  the  B~  ->  D">of  in  the  eP  decay  modes 
A IT''',  A'  ir‘*'ff'’  and  A'“ir+ir"T+.  In  (h)  through  (j)  the  B°  ->D**or  in  the 
rP  decay  modes  A'^nte®  and  K“ir‘''iT'7r+. 
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Figure  5.4  Background  contribution  from  £)*“  -►  D°->. 

In  (a)  the  D“  ^ in  (b)  the  rfl  -*  and  in  (c)  the  D°  -» 

It“ir'*'7T  JT"*  Er  subdeeay  modes. 


To  remove  this  background  from  the  B~  -t  D^p~  Mac  distribution  a cut 
on  cos(eft)  > 0.4  WM  on  imposed  on  ail  B candidates.'  This  requirement 
effectively  vetos  events  where  the  neutral  pion  produced  in  the  decay  of  the  p' 
IS  directed  in  the  backwards  directiou  with  respect  to  the  direction  of  the  p~. 
Eliminating  decays  which  contain  a slow  neutral  pion  agnlficantly  reduces  the 
background  contributions  from  D'®  ->  D®7  decays.  Figure  5.4  illustrates  the 
dramatic  drop  in  the  backgrounds  achived  by  vetoing  slow  ir®s.  The  plots  were 
made  by  taking  a Monte  Carlo  sample  of  20(10  B~  -t  -»  D°t 

events  and  reconstructing  them  as  D°p~.  The  solid  (dashed)  histograms  in 
(a)-(c)  show  the  background  levels  before  and  (after)  the  appbeation  of  the 
veto.  A reduction  of  approximately  90%  was  obuined  in  all  three  modes. 

WTiile  this  background  contributes  significantly  in  the  signal  region,  the 
background  shape  used  to  fit  the  distribution  accommodates  the  broad  peaks 

' The  cosine  of  the  helicity  angle  describes  the  poUriaation  of  the  p~  and 
forms  a cos^  distribution  (see  Section  4.0.1). 
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Erom  the  -*  baekgrouad  when  the  large  contiimum  component  is  in- 
cluded in  the  overall  background  shape.  As  a consistency-  check,  the  branching 
fractioas  for  these  S decay  mode  were  recomputed  using  the  yields  and  ac- 
ceptances corresponding  to  the  slow  jr®  vetoed  sample.  The  results  agreed  by 
better  than  one  standard  deviation  (sec  Chapter  6). 

In  Figure  5.3  (e  and  j)  some  excess  in  the  signal  region  is  also  found  for  the 
two  B -*  D’ai  modes.  By  looking  at  the  details  of  the  generated  events  in  this 
region  of  tlie  Mac  spectrum  (5.270  < Mac  < 5.286)  we  find  that  most  of  the 
events,  about  70%,  are  due  to  fl  ^ D*p®7r“  or  B -»  D'!r“»+ir'  decays.  This 
background  is  consistent  with  the  upper  limits  on  the  amount  of  non-rcsonant 
contribution  to  the  event  yield  in  B”  and  B®  -4  D'+af  discussed 

in  Section  5.2. 

5-1-2  Background  Contribution  from  Continiiiiin  Events 

The  continuum  background  contributions  is  the  largest  single  component 
of  the  overall  background  in  the  Mac  spectra-  The  fractional  contribution  of 
continuum  backgrounds  ranges  from  58%  in  the  B~  -*  D*®p“  decay  mode  to 
91%  m the  B"  -*  D®fr'  mode.  The  fraction  also  varies  u a function  of  the 
D®  subchannel,  typically  larger  in  the  cleaner  D®  -*  iv-ir+  modes  where  the 
contribution  from  backgrounds  combinatorics  is  small. 

The  Mac  spectra  for  off  resonance  data  are  shown  in  Figure  5-5  through 
5.7.  The  off  resonance  data  was  taken  at  a center  of  mass  energy  30  MeV  below 
the  T(4S).  Also,  the  integrated  luminosity  is  a factor  off  two  smaller  than  the 
on-resonance  sample.  To  compensate  for  these  differences,  the  distributions 
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Figure  5.5  Scaled  coutinuum  Mgc  spectrum,  B -*  Z)''*7r. 

In  (a)  the  B®  -*  D+Jt"  la  the  D*  decay  mode  D*  ->■  K-r+ir*.  In  (b)  through 
(d)  thefl"  -♦  D®ir  in  the  D®  decay  modes  K'jr'’',  and  A'”r+ir“ir‘'', 

respective^'.  In  (e)  through  (g)  the  B~  -»  D’®it-  in  the  D®  decay  modes 

K iT‘^iT®andA'“!r‘*’7r'a*.  In  (h)  through  (j)  the  B® D"*’a“  in  the 

D®  decay  modes  A*”tr'*'ff®  and  A”tr'*’it“7r'^, 
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Figure  5.6  Scaled  coutinuum  Mj^.  spectrum,  B -*  D^'h 


In  (a)  tUe  B°  -*  D+p"  in  the  D+  decay  mode  D*  ^ A-a+tr 
(d)  the  B~  ->  Z)®p-  in  the  decay  modes  A“ir‘*',  A'a+ir® 
respectively.  In  (e)  through  (g)  the  B~  -»  D‘“p-  in  the 
K a"*",  and  A''a‘‘'a“a*.  in  (h)  through  (j)  the 

[fi  decay  modes  A'a*.  A"a*a°  and  A“a*a“a‘^. 


T"*"-  In  (b)  througli 

f CP  decay  modes 
“ ^ D-*p-  in  the 
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MB^(GeV) 

Figure  5.7  Scaled  continuum  Mgc  spectrum,  S -*  £)<’)at. 

In  (a)  the  B"  -t  D+af  in  the  D*  decay  mode  D*  -t  K" ir+ir+.  In  (h)  tlirough 
(d)  the  B~  D®o,  in  the  £1®  decay  modes  K~tr*.  it  “jr+j®  and  A‘“jr''’ir“7r+ 
rrapectively.  in  (e)  through  (g)  the  B'  ^ D'°a7  in  the  D°  decay  modes 
^ and  A'“ff‘^7r'ir+.  In  (h)  through  (j)  the  B®  -*  D'‘*'ar  in  the 

D®  decay  modes  A'^ir*.  A‘”a*a®  and  A''ir‘‘'jr”>r'*^. 


Figure  5.8  AE  distributions  in  Monte  Carlo. 

The  Monie  Carlo  AE  spettra  for  (a)  B"  -*  D°)r-  and  (b)  B~  ->  D°p“  modes, 
were  scaled  by 


and  shifted  by  30  MeV. 

5.1.3  Mnr-  Background  Spectra  from  AS  SiHehaiid.; 

To  model  the  Mbc  backgrounds  in  on-resonance  data,  B mesou  candidates 
were  formed  using  the  AE  sidebands.  The  AE  sidebands  were  taken  on  oppo- 
site sides  of  the  AE  interval  centered  at  aero.  Events  selected  in  this  way  have 
the  same  number  and  type  of  final  state  particles,  with  similar  kinematics,  as  do 
signal  events  but  cannot  be  a product  of  real  B decays  since  their  reconstructed 
energy  does  not  equal  the  energy  of  the  beam.  While  the  similarities  indicate 
that  we  have  a relatively  good  background  model,  two  imiiortant  facts  make 


this  method  less  than  ideal  for  determining  the  overall  background  contribu- 
tion. First,  if  we  exclude  the  AE  regions  less  than  S.CtojiE  from  and 
the  region  of  AE  where  contributions  from  other  decay  inodes,  one  pion  mass 
removed  are  expected,  only  a small  range  of  AE  is  available  for  comparison. 
In  Figure  5.8  the  available  AE  regions  for  B~  -*  D^jr"  and  B~  -»  D°p'  that 
forms  the  AE  sideband  spectra  are  shown.  The  available  range  is  the  region 
between  the  large  peaks.  Secondly,  the  AE  distributions  are  biased  towards 
zero  by  the  procedure  of  picking  the  candidate  with  the  smallest  AE  a.s  the 
best  candidate  per  event  per  mode-  Both  of  the.se  conditions  lead  to  a Mjc 


background  levels  in  signal  Mac  distributions. 
The  plots  in  Figure  5.8  also  reveal  the  bocl 


AE  spectrum  from  decay  modes  other  than  the  one  being  reconstructed.  The 
plots  were  made  by  reconstructing  Monte  Carlo  simulations  of  various  B decay 
channels.  The  various  decay  modes  are  indicated  within  each  plots  ond  were 
reconstructed  as  either  a (a)  B~  dH,-  or  (b)  a S' -t  D°p-.  The  peaks 
centered  at  0 are  the  results  of  reconstructing  the  correct  Monte  Carlo  sample 
while  the  peaks  on  either  side  are  the  result  of  reconstnictiug  the  other  BB 
decay  modes.  ' 


malized  to  the  number  of  events  in  either  the  B~  ->  D^ir-  or  the  B®  -*  D*p~ 
Monte  Carlo  sample  to  insure  that  the  combined  plots  represented  the  expected 
backgrounds  In  data  without  the  continuum  component. 

In  Figures  5.9,  5.10  and  5.11  the  AE  sideband  Mgc  spectra  are  plotted 
for  all  non  color-suppressed  modes.  The  AE  sidebands  are  deBned  in  the  AE 
region  greater  than  or  equal  to  and  less  than  ±100  MeV.  The  100  MeV 
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Mjt.(GeV) 

Figure  5.9  AE  sideband  Mgc  spectrum,  B -*  D("*rr. 

In  (a)  the  6®  -►  D*n~  in  the  D+  decay  modeD*  -»  In  {b)  through 

(d)  thcB“  ->  D®ir“  in  the  D°  decay  modes  A’^rr+xOand  A'“x‘*‘x“x+. 
respectively.  In  (e)  through  (g)  the  B~  -t  D‘\~  in  the  Cfi  decay  modes 
K x-^,  I<  x*x“and  A— x*x*x+.  In  (h)  through  (j)  the  8° -*  D'+x"  in  the 
D®  decay  modes  A it*.  A“x^x®  and  A”x'*'x'x'''- 


,h.n.i  mMnfln 
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Figure  5.10  AE  sideband  Mbc  spettrum,  fl -»  D*')/ 


!d  (a)  the  B®  -t  D*p~  in  the  D+  decay  mode  D*  K~jt*z' 
(d)  theS“  -t  in  tbcD®  decay  modes  A"  tt'*',  K“jr+ir®  ! 
respectively.  In  (e)  through  (g)  the  B~  -*  D'°p*  in  the 
^ ir+,  and  A-*+ir-ir-^,  In  (h)  through  (j)  the  B“ 

D“  decay  modes  K ir*.  A'“7r‘'‘ir®  and  A‘”ir*7r“7r‘''. 


■*.  In  (b)  through 

rP  decay  modes 
->  D’+p-  in  the 
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Figure  5.11  AE  sideband  Mbc  spectrum,  B -* 

In  (a)  the  B®  -»  D‘*‘a'  in  the  D*  decay  mode  D*  -t  In  (b)  through 

(d)  the  B~  -»  D®oj"  in  the  Cfi  decay  modes  K^ir*.  K'jr‘‘‘!r®  and  A'“ir‘''ir”ir'*'. 
respectively.  In  (e)  through  (g)  the  B"  -»  in  the  D°  decay  modes 

I''”*'  iv^ir+it®  and  K“ir‘*'ir"jr+.  In  (h)  through  (j)  the  B®  -*  in  the 

EP  decay  inodes  K~w*.  K~ir*ji°  and  K~n*T~:r*. 
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cutoff  is  used  to  avoid  contamiuatiou  from  B mesou  derays  which  differ  by  one 
pion.  The  lOO  MeV  value  was  chosen  lo  accommodate  the  width  of  the  AE 
peak  in  the  region  one  pion  mass  removed. 

5.1.4  Mnr  Background  Snerrm  from  Wrone  Si.n  -H- 

The  wrong  sign  combination  Mgc  spectrum  was  also  used  to  model  the 
Mbc  backgrounds,  in  on-resonance  data  The  Mbc  distributions  of  wrong  sign 
“B  "candidates  were  formed  by  applying  the  algorithm  used  in  forming  real  fl 
candidates  but  with  the  light  meson  swapped  with  its  charge  conjugate.  This 

rately  model  contributions  to  the  background  from  decays  of  other  B mesons. 
This  technique  thus  models  background  contributions  from  random  combina- 
tions of  particle  but  suffere  in  its  ability  to  predict  contributions  from  real  3S 
backgrounds. 

shown  in  Figures  5.12,  5.13  and  5.14.  In  modes  where  the  decay  chain  pro- 
duces  a D®  an  additional  veto  was  introduced  to  eUminate  the  possibility  of 
contaminating  the  wrong  sign  Mac  spectrum  with  real  charge-conjugate  (B) 

if  both  the  K~  and  ?r+  from  the  D*’  decay  are  misidentified  as  a tr”  and 
respectively.  The  double  misidentification  would  lead  to  an  in  the 
number  of  events  in  the  signal  region  from  the  misidentified  (conjugate)  B de- 
cay. The  veto  used  eliminated  events  if  the  D“  reflection  mass,  defined  as  the 
invariant  mass  calculated  with  the  mass  assignments  reversed,  falls  within  30 
MeV  from  the  nominal  D°  mass-  The  Cfi  reflection  mass  veto  was  used  for  both 
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Figure  5.12  Wrong  sign  Mbc  spectrum,  B -* 

In  (a)  the  B*’  -t  D+ir"  In  the  D*  decay  mode  D*  -*  /v'ir‘*‘!r+.  In  (b)  through 
(d)  the  B"  -*  D^s~  in  the  Cfi  6ea.y  modfs  K~  t* . Ifir+fr'l  and  R‘“jr+ir-jr+, 
respectively.  In  (e)  through  (g)  the  fl'  ->  in  the  D“  decay  modes 

K IT+,  K'ir"*-7r®  and  K'ir‘''7r“ir+.  In  (h)  through  (j)  tbeB®  -t  in  the 

Cr  decay  modes  K'r+ir'^  and  K'Tr'tjr-jr*. 
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Figure  5-13  Wrong  sign  Mgc  spwtrum,  B -i 

In  (a)  the  B®  -*  D*p~  in  fhe  D*  detay  modeD'*"  -»  A'"ir*ir‘*'.  In  (b)  through 
(d)  the  B~  -*  lfip~  in  the  D®  decay  modes K“tr+rr®  and  K“ff'‘'rr”7r''' 
respectively.  In  (e)  through  (g)  the  B'  D‘“p-  in  the  D®  decay  modeti 
A ir+,  A-ir+ir®  and  A'-rr+tr-tr*.  In  (h)  through  0)  the  6®  -*  D‘*p-  in  the 
D®  decay  modes  K~n*.  A''ir+7r®  and  A'"7r‘*‘7r“ir'*'. 
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Figure  5.14  Wrong  sign  Mbc  spectrum,  fl -t  DWoi- 
In  (a)  tbe  B°  D+nf  in  the  D*  deray  mode  D*  ->  it  In  (b)  through 

(d)  theB”  -»  D^a”  in  theD°  decay  modes  K“jr'''.  i»“jr‘’'!r'’ and  it  ";r‘'‘:r“ff'''. 
respectively,  [u  (e)  through  (g)  the  B~  D'^af  in  the  D“  decay  modes 
JC-7T+¥®  and  i<-n+ir-T+.  In  (h)  through  (j)  the  B®  -»  D'+of  in  the 
D"  decay  modes  K“<r'''.  it“ir‘‘’iT°  and  ii'“ir+>r'ir"''. 


Figure  5.15  Reflection  D®  mass  spectrums  for  wrong  sign  “S"  candidate. 
Misidentificd  D®  mass  spert.ra  for  the  (a)  K'  n*  (b)  and  D°  -»  K~7r+  modes 
in  the  wrong  sign  B~  -»  D°ir-  decay  mode. 
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modes.  Id  the  latter,  the  momentum  spectrum  of  the  kaons  and  piooe  are 
softer  and  thus  panicle  ID  is  improved.  The  larger  number  of  daughters  also 
reduces  thelikeUhood  of  misidentification.  In  Figure  5.15  (a]  aud  (b)  the  peaks 
in  the  Zfi  reflection  mass  spectrum  at  the  D°  rnttss  for  the  D“  -t  K" )r+  and 
wrong  sign  ‘B  " candidates  are  shown-  The  plots  were  formed 
by  analyzing  on-resonance  wrong  sign  data,  requiring  events  to  satisfy  all  B 
candidate  requirements  and  to  fall  within  fi  MeV  of  tlie  Mbc-  I"  Figure  5.15  (c) 
and  (d)  the  Mbc  dislributious  for  WTOng  sign  data  with  (hatched]  and  without 
(unhatchcd)  the  reflection  moss  veto  are  shown. 


In  the  preceding  discussion  we  have  presented  the  Mgc  backgrounds  spec- 
tra separated  into  continuum  and  BB  components  and  presented  two  ways  of 
modeling  the  background  in  data,  in  this  section  we  discuss  the  conclusions 
which  were  drawn  from  the  analysis  of  these  background  Mgc  distributions- 
To  test  the  vahdity  of  the  assumption  that  the  Mbc  backgrounds  ronsist 
of  a continuum  and  a BB  component  we  compared  tbe  sum  of  the  scaled 
continuum  plus  the  BB  background  model  to  signal  distributions  by  superim- 
posing the  ptots.  The  continuum  component  is  taken  from  the  off  resonance 
data  scaled  a.s  described  in  Section  5.1.2.  The  BB  component  was  taken  from 
the  BB  Monte  Carlo  (antiskim)  sample  and  was  normalized  to  tlie  number  of 
on-resonance  BB  events.  Figure  5-16  reveals  that  the  sum  accounts  tor  all  of 
the  background  present  and  that  both  tbe  shape  and  the  the  normalization  are 
adequate  modeled  by  this  background  estimate. 

To  insure  that  the  iE  sidebands  model  the  shape  of  the  true  background 
the  Mbc  siiertrum  was  formed  with  the  iE  sideband  tedmique  on  the  full  flS 


Mb^(Gc;V) 


M„c(OeV) 


Figiire  S.lfi  SuporpositioD  of  contmuum  and  55  background  componpnts- 
The  dashed  histogram  shoa-s  the  sum  of  the  scaled  continuum  plus  BB  back- 
ground components,  the  hatched  histogram  shows  the  continuum  componeut 
superimposed  on  the  signal  Mgc  spectra.  The  B decay  modes  shown  arc  the 
(a)  B~  -*  D®ir-,  (b)  B~  (e)  B*  -►  D'°ir'  and  (d)  S'  D’°p' . 

All  three  of  the  CP  decay  modes  are  included. 
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Figure  5.17  Bfl  compared  to  AE  sideband  barkgrounds. 

A Monte  Carlo  simulation  of  the  Mgc  distributiora  of  flS  background  (solid 
^Btogram),  compared  with  the  AE  sidebands  from  the  same  Monte  Carlo 
(points).  All  events  in  the  signal  modes  have  been  removed  from  the  BB  Monte 
Carlo  simulation,  (a)  B~  -*  (bl  B'  -»  D°n-  (d  B'  O’Oir-  —H 

(d)  a-  -e  ' ' 


Mgf.(CieV) 


Mac(GeV) 


Figure  5.18  BB  plus  tominuum  compared  to  iE  sideband  backgrounds. 
Mbc  distributions  for  iE  sidebands  in  on-resonanre  data  (solid  histogram) 
compared  to  tlie  sum  of  continuum  data  aud  a Monte  Carlo  simulation  of  BB 
background  (poinw)  (a)  B~  -*  D°w-.  fb)  B~  ->  rfip-.  (c)  B~  -*  D’^tr" 
and  (d)  B~  -i  C^V- 
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Monte  Carlo  sample.  These  distributions  were  then  compared  to  our  moilcl 
of  the  true  background,  the  Bb  aniiskim  Mbc  spectra.  The  comparison  was 
made  bj-  superimposing  the  antiskim  aud  the  iE  sideband  distributions,  see 
Figure  5.17.  The  figures  reveal  that  the  shapes  of  the  distributions  match 
rather  well  and  thu-s  AE  sideband  technique  accurately  modek  the  shape  of 
the  true  background.  As  an  additionai  check  the  AE  sidebands  distributions 
for  on-resonance  data  were  compared  to  the  sum  of  the  scaled  continuum  plus 
the  BB  Moute  Carlo  background  Mbc  distribution  (Figure  5.18).  lb  allow  for 
direct  comparison  of  tbe  background  shapes,  the  overall  normalization  of  the 
combined  BB  and  continuum  distributions  was  allowed  to  float.  The  shape  of 
the  AE  sideband  spectra  also  matches  the  shape  of  the  sum. 

Several  functions  were  used  to  represent  the  background  dato,  including  a 


energy  used  by  tbe  ARGUS  [44]  experiment  and  tbe  CLEO  background  shape. 
The  CLEO  background  shape  was  found  to  fit  the  background  distribution 
best.  The  CLEO  background  shape  is  defined  as  a sttaighi  line  with  a parabolic 
roll-off  at  the  Mac  endpoint.  Tlie  functioual  from  is  described  by 
fmu  + 4. 

k(x)  = \ ^ (5,2) 

[ mio{l  - (i  - xi,]^  -Hi,}  Xo<x<E^ 

where  Ej  is  the  beam  energy  and  i,  = i*  the  point  where  the 

parabolic  roU-ofl  begins.  Figures  5.19  and  5.20  show  the  CLEO  background 
shape  fit  to  AE  sideband  aud  wrong  sign  background  distributions.  They  show 
that  both  distributions  are  fit  adequately  by  the  shape.  The  offset  was  fixed 
at  8 MeV  on  all  nou  color-suppressed  modes.  Ip  color-suppressed  modes  the 
D&et  was  varied  and  the  value  which  gave  the  most  conservative  upper  limit 


Figure  5.19  CLEO  background  function  fit  to  AS  sidebands  Mac  spectra. 
In  (a)  the  fit  to  the  B~  , (b)  fl-  _ [fip-  . (c)  B~  D-Otr'  and  (d) 

^ eath  plof  ^ distributions.  AU  three  ZP  decay  modes  are  included 


Figure  5.20  CLEO  background  tunction  fit  to  the  wrong  sign  Mgc  spectra. 
In  (a)  the  fit  to  the  B~  , (b)  B"  -»  D^p-  . {c)  B“  and  (d) 

fl-  D'®p“  Mbc  wrong  sign  distribution.?.  .411  three  of  the  [fi  decay  inodes 
ore  included  in  earh  plot. 
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'I^ble  5>1 


5.1.6  Systematic  Errors  from  Background  Shanes 


Systematic  errors  due  to  the  background  .shape  were  determined  by  com- 
paring the  Mbc  ytolds  obtained  with  diSerem  CLEO  background  function  pa- 
rameterizations.  The  yields  were  determined  four  times  for  each  mode.  Each 
time  the  Mgc  distribution  was  6t  using  different  values  for  the  background  pa- 
rameters. First,  all  sTiriable  background  parameters  were  allowed  to  float.  The 
parameters  were  then  fixed  to  three  different  sets  of  values,  each  obtained  by 
fitting  n particular  background  model.  The  background  models  used  were;  the 
sum  of  BB  Monte  Carlo  and  scaled  continuum  with  both  normalizations  fixed, 
the  AE  sideband  and  the  wrong  sign  combination  oii-resonance  data.  The  er- 
rors from  the  background  contribution  were  obtained  for  each  mode  by  finding 
the  standard  deviation  of  all  four  yield  determinations.  The  standard  devia- 
tion and  the  yields  determined  with  all  the  background  functions  are  listed  in 
Tbble  5.1. 


In  decays  to  multiple  pion  final  states  the  non-resonant  contribution  to  the 
number  of  observed  events  were  determined  by  comparing  the  invariant  moss 
spectrum  of  the  composite  final  state  to  the  expected  distributions  for  reso- 
nant and  non-resonant  decays  in  Monte  Carlo  simulations.  The  invariant  muss 
distributions  were  formed  by  first  reconstructing  the  events  as  described  in 
Chapter  4 without  the  invariant  mass  cut  ou  the  “light"  meson  daughter.  The 
invariant  mass  plots  are  then  formed  by  taking  events  that  fall  within  the  Msc 
signal  and  subtracting  tbe  background  taken  from  the  Mgc  sidebands.  The 
sidebands  were  defined  to  be  between  5.20  GeV  and  5.26  GeV.  The  Mgc 


Figure  S.21  Resonanl  and  non-ccsonant  einmiatioa  ol  B ->  Da\  decays. 
The  iDvniiant  njrif  mass  in  Monte  Carlo  simulations  of  (a)  B~  ->  0®of  (ps) 
and  (b)  B®  -*  D^a”  (pn)  modes.  The  solid  histogranis  shows  the  distribution 
when  the  rmr  are  the  decay  products  of  an  resonance.  The  hatched  his- 
togram shows  the  resulting  trirs  distribution  when  the  particles  are  first  formed 


sideband  invariant  mass  spectrum  was  scaled  before  the  subtraction.  The  bach- 
ground  subtracted  pint  was  then  fit  to  the  sum  of  resonant  and  non-resonanl 
invariant  mass  distributions  whose  shapes  were  determined  from  Monte  Carlo 
with  Che  normalizations  taken  as  free  parameters  in  the  fit. 

To  illustrate  the  method  we  show  the  procedure  as  applied  to  the  extraction 
of  the  relative  ooo-resouant  (p^ir”)  contribution  to  the  B~  -*  and  S®  —* 
decays.^  The  n^  invariant  mass  distributions  for  resonant  and  non- 
resonant Monte  Carlo  are  shown  in  Figure  5.21.  These  histograms  were  used 
to  fit  the  background  subtracted  oj  invariant  mass  plots  (Figure  5.22  (a)  and 

' Tbe  shape  of  the  ±x'*‘ir”  iiivarieiit  mass  distribution  Is  very  similar  to  the 
p^jr~  distribution. 


Mgt,(GeV) 


Figure  5.22  Nou-resoDaot  contribution  to  £ -t  Doj  decays. 

In  (a)  and  (b)  the  background  subtracted  jr”jr'*'tr”  invariant  mrcsa  distributions 
for  the  B"  -*  D^af,  D®  -t  f»-(r+  and  the  B®  ^ D+o^  , K'- x+x+ 

modes.  The  black  squares  represent  the  background  subtracted  data,  the  solid 
and  batched  histogrants  are  the  resuit  of  the  fit  from  the  Monte  Carlo  distri* 
butions.  In  (c)  and  (d)  the  Mgc  distributions  are  shown  for  tlie  B~  -e  D®n^ 
and  B®  -*  0*0^  with  a cut  in  the  fii  poor  region  of  the  n ?r'*'s”  spectrum. 


(b)).  Imegr&tmg  tbe  DQn<resou>uit  shape  (hatehed  histograiD  iq  Tigure  5.31 
(a-b)  we  find  thai  5.9±3.G  events  tall  between  I.O  Ge\'  and  1.6  GeV  in  the 
B~  -»  mode  and  19.6±d.O  events  in  the  S"  -*  D*aj  mode.  This  implies 
that  non-resonant  decays  contribution  to  the  B~  — t and  -f 

are  less  then  11.55(  and  34.35(  at  the  905t  confidence  level. 

As  a constancy  check  the  Mpc  distributions  were  rcplotted  but  now  tbeui 
invariant  mass  was  required  to  lie  in  the  aj  poor  region  defined  to  be  between 
2.0  GeV  and  2.6  GeV.  The  Mbc  distributions  when  selecting  C13S  from  tbe 
poor  region  are  shown  in  Figure  5.22  (c-d)  for  l)Otb  the  B~  ->  D°aJ"  and 
8*’  -i  D'''a^  respectively-  Tbe  number  of  events  found  in  the  Mbc  signal 
region  were  1.3±6.2  for  the  8"  -»  D^aJ"  and  18.8  ±10.1  for  the  8®  -•  D*oj". 
Tbe  number  of  non-tesonant  contribution  using  this  method  gives  less  than 
18%  aud  36%,  at  tbe  90%  confidence  level  for  8"  — * and  8®  — » D'*’ap 

These  values  are  consistent  with  tbe  upper  limits  found  by  fitting  tbe  three 
pion  invariant  mass  distributious. 

A similar  analysis  was  performed  on  tbe  other  resonant  B meson  decay 
modes  the  B~  -t  8“  -*  D’+oJ'  and  D<’>+p*  by  M-S- 

Alnin  e(  of.  [45].  The  teclinique  used  differed  only  ia  that  the  background 
subtracted  resonant  mass  plots  were  fit  with  a Breit-Wigner  and  polynomial 
parameterizations  determined  from  Monte  Catio  simulations.  Tbe  results  ob- 
tained showed  that  the  non-resonaut  component  10  tbe  B — » modes  are  less 

than  9.0%.  and  leas  than  2.5%  for  the  3 — * Dfl  modes.  For  the  8®  — ♦ 
and  J3”  — t tbe  contributions  from  non-rcsonant  7r~p^  and  D’“/i  were 

found  to  be  less  9.4%  and  10.6%  at  tbe  90%  confidence  level. 


5.3  Backgrounds  in  Color-Supptesscd  Modes 


5.3.1  BB  Backgrounds 

In  order  to  examine  tbe  BB  background  contributions  to  color  suppressed 
modes,  no  signal  mode  subtraction  scheme  was  employed  since  the  Monte  Carlo 
generator  decay  Cables  do  not  contain  appreciable  branching  Iraction  to  these 
modes.  The  generator  docs  however,  uses  a branching  fractions  of  1 x 10"’  in 
some  of  the  lower  multiplicity  modes,  such  as  B**  ->  D®ir®  and  -*  D°p®. 
With  a sample  of  2 x 10®  events  and  low  x B^g  these  small  rates  do 
not  appear  in  the  plots.  However,  in  the  B®  — • mode  the  Monte  Carlo 

efficiency  Is  somewhat  larger  and  a few  events  do  appear  in  the  signal  region. 
We  verified  that  these  events  did  come  from  the  small  generator  -t  D®ir® 
branching  fraction  by  inspecting  a dump  of  the  Monte  Carlo  events  in  the  signal 
region.  In  Figure  5.23  and  5.24  the  Mgc  distributions  for  each  of  the  color- 
suppressed  BB  decay  modes  are  shown.  The  figures  reveal  that  no  significant 
contributions,  other  than  that  already  discussed,  are  present  from  other  known 
BB  decays. 


S.3,2  Continuum  Backgrounds 

The  continuum  background  contributions  are  shown  in  Figure  5.25  and 
5.26.  These  plots  were  created  as  previously  by  using  off-resonance  data  and 
scaJing  the  Msc;  distributions  by  the  appropriate  factor,  Equation  (5.1).  As 
with  the  BB  background  no  significant  contributions  in  the  »gnal  region  ore 
observed. 


Mji-  (GeV) 

Figure  5.33  BB  backgrounds  in  color-suppressed  decsj'S  B -t 
Generic  BB  Monte  Carlo  data  used  in  reconstructing  the  B -*  decay 
ntodes  (a<)  0“^“,  (d-f)  D\  (g-i)  (j-l)  and  (m-o)  Cfiu.  The  first 

columns  corresponds  to  the  K'”?:'*'?:®  and  decay  modes  of 

the  C“  respectively. 


Figure  5-24  BB  bockgrounds  in  color-Buppressed  decays  B — * 

Generic  BB  Monte  Carlo  data  used  in  reconstructing  the  B -»  £>"®h®  de- 
cay inodes  (a-c)  D'"ir'>.  (d-f)  D'%,  (g-i)  Q-l]  D'^i>°  and  (m-o) 

The  first  columns  corresponds  to  the  K~w'*'.  K”j'*"ir®  and  7C”jr+n")r+  decay 
modes  of  the  D**  respectively. 


130 


Mgc(GeV) 

Figure  6.25  Continuum  background'!  in  color-suppressed,  B -* 

Scaled  and  shifted  continuum  data  recoiBtnicted  in  B®  -►  6®b®  decay  modes 
(a-c)  £)®)r®,  (d-f)  D®ti,  (g-i)  ifiri.  (}-l)  D°p®  and  (m-o)  D®oi.  The  first  columns 
corresponds  to  the  and  K~s*s~s*  decay  modes  of  the  LP 

respectively. 


Mj(.  (GeV) 

Figure  5.26  Cooiiuuum  backgrounds  in  color^suppressed.  B -¥ 

Scaled  and  shifted  coDtiDuiim  data  reconstructed  in  — r D*^bP  decay  modes 

(a-c)  (d-f)  (g-i)  0-1)  D’“p°  and  (m-o)  D'®-;.  The  first 

rolumoa  correspoDcU  to  theK'ir"'',  /v'7r’'7r°  and  K'7r'''jr“ir'''  decay  modes  of 
the  rP  respectively. 


CHAPTER  6 

RESULTS  AND  DISCUSSION 


Iniroductiop 

In  this  chapter  we  provide  a detailed  description  of  the  experimental  mea- 
surements and  discuss  these  results  in  the  context  of  the  theoretical  framework 
introduced  in  Chapter  2.  The  measuremems  consist  of  branching  fraction  and 
polarisation  measurements,  both  pertormcd  b>'  counting  the  number  of  events 
found  in  particular  distributions  and  correcting  the  yields  by  appropriate  ef- 
ficiencies and  decay  daughter  branching  fractions.  The  results  are  then  used 
in  tests  of  the  factorization  hypothesis,  color  suppression  and  to  determine  the 
values  of  the  BSW  parameters  0|,  nj  and  the  relative  sign  of  ns/n]. 


6.1  Determination  of  Branching  Fractions 
The  fi  branching  fractious  were  determined  by  using  the  event  yields  for 
each  6 decay  mode,  the  product  of  the  total  number  of  BB  events  in  the 
data  sample,  the  branching  fraction  of  the  daughters  and  the  reconstruction 
efficiency  (cmcs)'  determined  form  Monte  Carlo  studies.  The  explicit  formula 
used  to  calculate  the  branching  fractions  was 


(6.1) 


The  index  b specifier  the  particular  B decay  mode  and  the  index  dia  the  product 
over  aii  the  reievojit  hranchiiig  fraction  of  the  B daughters  in  the  particular 
decay  chain. 

The  D and  D“  branching  fractions  are  taken  from  recent  CLEO  mea- 
surements. The  latest  CLEO  results  represent  the  most  recent  and  accurate 
measurements  made  to  date.  Also  by  asing  the  CLEO  measurements,  the  to- 
tal systematic  error  can  be  reduced  by  removing  tracking  errors  common  to 
this  and  other  CLEO  analyses.  The  branching  fractious  for  the  D‘*  and  D‘^ 
modes  are  [-)6| 

6(D’‘'  -»  = 63.6  ± 2.3  ± 3.3% 

B(D’+  ^0°:?*)  = 68.1  ±1-0  ±1.3% 

The  branching  fractions  for  the  D*’  and  D*  decay  modes  are  normalized  to  the 
-*  K~ir*  mode.  The  value  of  the  uormalizatioo  mode  and  the  ratios  are 
B(£P  K"ir‘*')  = 3-91  ± 0.08  ± 0.17%  |47| 

B(D°  -»  K'ir+ir‘')/B(D‘’  -*  K~s*)  = 3.78  ± 0.07  ± 0.49  [48| 

-»  fC-ir^ir-ir+j/SfD®  ^ ff“7r+)  = 2.05  ± 0.11  (2) 

B(D*  -»  K'!T*tr+)/B(D°  ^ /v*ir+)  = 2.35  ± 0.16  ± 0.16  [49] 

We  use  the  ratios  of  D branching  fractions  to  determine  the  B branching  frac- 
tions since  these  are  known  to  greater  accuracy  than  are  the  absolute  values. 

6.1.1  Decays  to  Nou  Color-Suppressed  Modes 

The  event  yields  for  the  non  color-suppressed  decays  were  determined  by 
frtcing  the  Mbc  distribution  for  events  that  passed  ail  the  B meson  recon- 
struction requirements  discussed  in  Chapter  4.  The  distributions  were  fit  to  a 
Gaussian  of  fixed  width  plus  a background  shape.  The  width  of  the  Gaussian 
was  fixed  to  2.64  MeV  for  all  modes.  This  value  was  determined  from  Monte 


significaDtly 


3de.  The 


background  shape  used  was  the  straight  line  with  parabolic  roll-off  discussed 
in  Section  S.l.S. 

In  Figures  6.1,  6.2  and  6.3  the  on-resonance  Mgc  dieuibutions  for  B decay 
to  one,  two  and  three  pion  light  meson  states  are  shown.  In  each  plot,  theffsed 
width  Gaussian  and  the  background  function  are  superimposed  on  the  Mbc 
distributions.  The  amount  of  background  observed  in  each  mode  is  commen- 

exatnple,  the  ^gniffcant  reduction  of  background  levels  fur  modes  which  decay 
to  a vector  charm  meson  is  due  to  the  mp.  — cut  which  takes  advantage 
of  the  small  in  D‘  -*  De  decay. 

The  branching  traction  results  for  all  the  non  color-suppressed  modes  are 
given  in  Tables  6. 1 through  6.3.  Three  values  were  determined  tor  each  decay 
mode  that  included  a in  the  decay  chain,  one  for  each  of  the  D^'  subchannels. 
For  these  modes  the  values  were  averaged,  weighted  by  tbeir  statistical  errors. 
Two  systematic  errors  are  quoted  for  tbe  averaged  results.  The  second  error, 
tbe  first  systematic  error,  was  calculated  by  adding  in  quadrature  the  following 
quantities:  tbe  background  errors  established  as  described  in  section  5.1.6,  tbe 
error  due  to  Monte  Carlo  statistics,  typically  less  than  296,  and  the  uncertainty 
in  detection  eSciencies  or  tracking  errors.  Uncertainties  in  detection  effidendes 
for  charged  tracks  with  momenta  greater  than  225  MeV  are  within  ±2.0%  of 
the  measured  efHdency.  For  charged  tracks  with  momenta  below  225  MeV’,  the 
uncertainty  increases  to  ±5.0%.  For  neutral  pions  tbe  uncertainty  is  ±5.0%. 
These  errors  were  assigned  to  all  the  final  states  particle  in  the  decay  chain. 
Tbe  errors  for  charged  tracks  were  suiurued  linearly  and  the  result  was  added. 
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Figure  6.1  Mbc  distributions  for  the  B — * modes. 

In  (a)  the  -»  D+e*  in  the  D+  decay  mode  D+  -*  . In  (b) 

through  (d)  the  B~  -*  D®ir”  in  the  Cfi  decay  modes  it'tr'*'.  and 

K^ir+ir'ir*  respectively.  In  (e)  through  (g)  the  S"  -t  in  the  D°  decay 

Diodes  K“ir^,  and  K~7r'*ir‘3r*.  In  (h)  through  (j)  the 

in  the  Hfi  decay  modes  JC”ir'*'.  and 


Mbc(GcV) 

Figure  6.2  Msc  distributions  for  the  B — * modes. 

In  (a)  the  B®  — r D*eT  in  the  O'*  decay  mode  D"*"  -»  K“7r'*’rr'*’  . In  (b) 
through  (d)  the  B"  -t  D®p"  in  the  Zfi  decay  modes  K”!?"*",  and 

K“ir‘''rr“7r*  respectiveiy.  la  (e)  through  (g)  the  B~  -*  D*®p“  in  the  D®  decay 
modes  K”jr'*'ir^  and  In  (h)  through  Q)  the  — t D"*"p“ 

in  the  D’’  decay  modes  A'" it* , K"ir"''a°  and  K'Tr+ir'ir'*’. 
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Mji-CGeV) 

Figure  6.3  Mgc  distributions  for  tbp  B -*  modes. 

In  (a)  the  B°  -*  D*of  in  the  D+  decay  mode  D*  -*  K'-s+»*  . In  (h) 
through  {dj  the  B~  -t  in  the  D®  decay  modes  K~’r'*’,  and 

K'x'*'x”x'*'  respectively.  In  (e)  through  (g)  the  B~  in  the  D®  decay 

modes  K”x'*',  Jt'x'^x®  and  K”x'*'x”x'*'.  In  (h)  through  (j)  the  — ♦ D“7‘a]‘ 

in  the  ifi  decay  modes  K"x+,  K'x+x®  and 
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Itible  6.1 

Table  of  braiichiag  fractioDS  B — » decays 


B Mod: 

D Mode 

Yield 

eme(%) 

B(%) 

fl(%)  average 

D*ir- 

K"ir+x+ 

170 

27.8 

0.308±0.02S 

0.308  * 0.026  * 0.028  * 0.031 

K-X+ 

52 

31.0 

0.287  ±0.041 

D'+ff- 

K-x+x® 

S5 

12.3 

0.314±0.047 

0.304  * 0.024  * 0.025  * 0.027 

K-X+X-X+ 

46 

12-2 

D.313±0.050 

K-X+ 

162 

38.3 

0.496  * 0.042 

D«ir- 

K-x’-x® 

306 

17.4 

0.544*0.052 

0.534  * 0.025  * 0.033  ±0.047 

K-x*x-x* 

134 

19.1 

0.580*0.052 

If-x-^ 

33 

U.l 

0.429*  0.081 

D'°ir- 

K"x+x® 

73 

7.5 

0.475*  0.079 

0.497*  0.044  * 0.048  * 0.057 

K--X+X-X+ 

SI 

7,0 

0.654*0.101 

Tkble  6.2 

Table  of  branching  Craciions  B -*  decays 

B Mode 

D Mode 

Yield 

Cme(% 

B{%) 

B(%)  average 

D*r 

K-X+X+ 

171 

9.9 

0.861  *0.078 

0.861  * 0.078  *0.109  * 0.086 

K-X+ 

51 

lo.a 

0.816*0.125 

D‘*p- 

K-x+x® 

84 

4.2 

0.897  * 0.110 

0.844  * 0.071  * 0.096  * 0.076 

K-x+x'x-^ 

44 

4.7 

0.794  * 0.137 

K-X+ 

123 

14.1 

1.081*0.116 

DOp- 

K-x+x® 

206 

6.8 

0.934  * 0.099 

1.022  * 0.067  * 0.109  * 0.086 

K-x+x- X* 

126 

6.9 

1.119*0.146 

K-X+ 

46 

5.6 

1.508*0.240 

D'^p- 

K-x+x® 

76 

2.5 

1.498*0.204 

1.444*0.134*0.188*0.161 

K-X+X-X+ 

33 

2.4 

1.275  * 0.265 
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I^ble  6.3 

Table  of  braodiiog  fractions  B -*  decays 


B Mode 


D Mode 


37 

27 

47 


g(%) 

0.757±0.063 
1.569  * 0.270 
0.921  ±0.1S2 
1.461*0.309 
0.812*0.215 
0.641*0.198 
0.740  * 0.241 
1.685*0.453 
2.132*0.454 
1-874*0-487 


S{%)  average 

0-757  * 0.065  * 0.116*  0.076 


1.205  * 0.140  * 0.138*  0.09S 


0.812  * 0.215  * 0.151  * 0.022 


1.898*  0.268  * 0.236  * 0.221 


in  quadrature,  to  the  total  meuirai  pion  uncertainty,  the  second  systematic 
error  includes  the  uncertainties  in  the  D and  D'  branching  fractions. 

In  Figure  6.4  the  branching  fractions  are  grouped  to  show  the  measurement 
taiiatiODS  within  each  B meson  decay  containing  a D“.  The  variations  between 
the  measurements  all  fall  well  within  2.0  standard  deviations  of  the  statistical 

In  all  three  fl”  -+  modes,  (Figure  6.2  b-d)  the  background  levels  in 
the  Mjjc:  distributions  are  considerable.  Tb  insure  that  Che  large  background 
did  not  significantly  afiect  the  measurements,  we  recalculated  the  branching 
fractions  using  only  the  half  of  the  p~  helicity  range  which  corresponds  to 
the  fast  p~  decay  daughter.  This  cut  removed  a significant  amount  of 

background  including  the  background  associated  with  B~  — » with 

jyO  which  forms  a broad  peak  in  the  signal  region  (see  Section  5,1,1). 

The  B~  — » Mgc  distribution.^  with  the  slow  tt^  veto  are  shown  in  Fig- 

ure 6.5.  The  branching  fractions  obtained  by  fitting  these  distributions  gave 
values  consistent  with  the  values  obtained  with  the  full  helicity  range.  The 
values  were  found  to  differ  by  less  than  0.4o,  where  c denotes  the  statistical 
error  (see  Figure  6.4  (e)). 

For  the  B~  -*  D®aj"  decays,  only  the  D®  -*  K~7r*  mode  is  used  in  deter- 
mining the  branching  fraction.  The  statistical  significance  of  D®  — r 
submode  is  4.7o  compared  to  less  than  3.0ir  for  the  D®  — * K’jr'*'?^  and 
D®  — * K”jr'^jr”jr'*'  channels.  In  addition,  the  large  background  levels  found  in 
the  O®  -»  and  -►  modes  would  lead  to  large  uncer- 

tainties in  the  shape  of  the  background  function  and  would  further  reduce  Che 
significance  of  the  signal. 
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Figure  €.4  Plots  of  tbe  d branching  fractions  with  all  three  ^^ubmodee. 
Jn  (a)  the  S'  -4  Efiit'  modes,  in  (h)  the  B~  -t  D'^x"  modes,  in  (c)  the 
modes,  in  (d)  the  B~  — t D^p'  modes  with  the  full  helicity 
(filled  triangles),  and  half  the  helicity  (open  triangles),  in  (e)  the  S'*  -•  D'*p~ 
modes,  in  (g-i)  the  0'  -t  D®af , the  B"  -*  and  the  8“  -t  D'+af 

respectively.  The  enor  bars  conespood  to  the  la  statistical  error 
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Figure  6.5  Mgc  spectrum  for  B~  -♦  D^p~  after  appUcatiou  of  bebcity  tui. 
Id  (a)  the  D®-fK-ir+.  In  (b)  the  D°  -i  Jv-tt+tt'’  and  in  (c)  the  D° 

6.1.2  Decays  to  Color-Suimressed  Modes 

Clear  signals  were  not  observed  for  any  of  the  color-suppressed  distri- 
butions investigated.  The  lack  of  a prominent  enhancement  pcevented  us  from 
extracting  brandling  ftactions  and  thus  only  upper  limits  were  calculated  for 
each  colot-suppressed  B decay  channel.  To  simplify  the  analysis  and  to  insure 
that  the  upper  limits  obtained  do  not  undetestimate  Che  branching  fractious, 
the  yields  were  measured  for  the  sum  of  the  three  D^'  subdecays  in  each  B 


To  find  the  upper  limit  on  the  number  of  events,  the  Mgc  distributions 
were  fit  to  a Gaussian  of  fixed  width  plus  the  background  function  described 
in  Section  5.1.5.  but  with  additional  requirements  imposed  on  the  function 
parameters.  The  mean  of  the  Gaussian  was  fixed  to  5.28  GeV  and  cbe  slope 
of  the  CLEO  background  function  was  fixed  to  zero.  Also,  the  offset  of  the 

limit  on  the  number  of  events  obtained  tbc  most  conservative  value.  The  upper 
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Figure  6.6  Mjc  spectra  for  color-suppressed  B decays,  D®  modes  com- 

In  (a)  the  B°  ->  D®ir®  mode,  in  (c)  the  fi®  -*  in  (c)  the  B°  V°r,. 

in  {d)  the  e®  D-®».  in  (e)  the  fl°  -r  D®n',  in  (f)  the  fi®  -r  J>*®p',  in  (g) 
the  S®  -»  D®p®,  in  (h)  the  B®  D‘®p®,  in  (i)  the  B®  D®u  and  in  0)  the 

8®  -t  D’®oi. 


Figure  6.T  Mgc  for  color-suppressed  B 

Each  row  corresponds  to  a different  B decay  mode,  In  (a-c)  the  three  B®  -t 
D”jr^  modes,  in  (d-f)  the  S®  -•  modes,  in  (g-i)  the  S°  -»  D®i>'  decay 
modes,  in  (j-1)  the  fl®  -*  D®p®  and  in  (m-o)  the  B®  -»  D®w.  The  columns 
correspond  to  different  I?®  decay  modes,  from  left  to  right  the  D®  — » 

D®  -t  A-tr+rr®,  and  the  D®  K-ir+ir-ir’’. 
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Figure  6.8  Mgc  spectra  for  color-suppressed  B — » 

Each  row  corresponds  to  a different  B decay  mode.  In  (a-c)  the  three  S®  -* 
D'^ir^  modes,  in  (d-f)  the  -*  D‘®i?  modes,  in  (g-i)  the  -*  D"®!/  decay 
modes,  in  (j-1)  the  6°  -•  D'®p®  and  in  (m-o)  the  fl°  -+  D'®u).  The  columns 
correspond  to  different  D"®decay  modes,  from  left  to  right  the  X)®  — * K”a'*‘, 
0®  K-n*ir°,  and  the  D®  K-ir+ir-ir+. 


Table  of  braoching  (tactions  for  color-suppressed  decays 


Decay  Mode 

Nob. 

Cmc{%) 

a (%)  at  90%  C-L. 

SO-rD^ir® 

<33.3 

29.2,15.3.10.8 

< 0.033 

<9-4 

26.6,9.7, 13.1 

< 0.033 

<2.3 

13.6, 6.9.8.4 

<0.029 

g0_j£>0^ 

<33-7 

19.0, 7.0,10.9 

<0.060 

^ D°ui 

< 13-0 

11.0.3.2,3.9 

< 0.057 

S“  -» D-'’tr° 

< 14-6 

12.2,6.5.6.6 

< 0.055 

-»  D-Of) 

<3.6 

9.6, 4.1. 4.9 

< 0.050 

fl«  -» 

< 24 

5.6, 2-3,2.5 

<0.13 

< 19.1 

7.4,4.0,3.6 

<0.12 

fiO  -t  D’°<j 

< 11.8 

5-1. 2.6, 2.7 

<0.12 

(or  each  Cfi  submode  and  the  upper  limits  on  the  branching  fraction  are  listed 
in  Thble  6.4.  The  Mgc  distribution  for  the  color-suppressed  modes,  with  all 
three  Lfi  submodes  combined  are  shown  Figure  6.6.  In  each  plot  both  the  fit 
(solid  line)  to  the  Mgc  distribution  and  the  fit  with  the  upper  limit,  at  the 
90%  confidence  level,  (dashed  line)  are  superimposed.  In  Figures  6.7  and  6.8 
the  distribution  o(  each  color-suppressed  B decay  mode  are  separated  into  the 
three  Efi  subdecays. 

The  upper  limit  for  the  branching  fractions  were  determined  by  a modified 
version  of  Equation  (6.1), 
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where  the  value  of  Aoi,j  is  now  the  upper  limit  on  the  number  of  events  summed 
over  the  three  subdecays  and  the  Monte  Carlo  efficiency  term  is  replaced 
by  the  sum  of  the  efficiencies  limes  D°  branching  tractions.  The  upper  limits 
were  determined  by  the  method  described  in  section  III  of  the  Particle  Data 
Group  [2]. 

While  no  clear  signals  were  obtained  Ibr  any  of  the  Mgc  distributions. 
Figure  6.7  (a)  shows  a signiBcant  enhancement  in  the  signal  region  of  the  D®  -* 
submode  in  the  B®  -t  D®ir®  channel.  Since  signals  were  not  obtained 
in  cither  the  D®  -*  K~-n*  nor  the  D®-*  modes  we  scrutiniaed 

the  events  in  all  three  modes  for  possible  errors  in  the  event  reconstruction 
procedure.  No  obvious  errors  were  found,  however  we  did  find  that  in  the 
D®  — t K”*"*"  mode  approximately  16%  of  the  events  in  the  signal  region  were 
double  counted  events.  Double  counting  occurs  when  the  daughters  of  the 
D®.  the  K~T*.  are  misldentified  as  ir~  K"*"  and  reconstruct  to  a D®.  The 
misidentified  D®  Is  then  combined  with  the  same  x®  used  to  reconstnictiou 
the  fl®  and  thus  the  event  is  counted  twice.  The  double  counting  does  not 
significantly  alter  the  final  answer  since  the  Monte  Carlo  efficiencies  are  affected 
in  exactly  the  same  way.  If  the  double  counted  events  are  removed  from  the 
plot  there  is  still  a significant  excess  in  the  signal  region  which  cannot  be 
accounted  for  by  either  the  continuum  background  Figure  5.25  (a)  nor  the 
known  BB  background  shown  in  Figure  5.23  (a).  Further  work  on  this  mode 
is  in  progress. 


6.2  Measuremeut  of  Polarisation 

Fully  recoDstructed  B decays  were  used  to  measure  the  polaiiaations  of 
the  final  states  in  B~  -*  D‘®p”  and  6®  -♦  decays.  The  fractional 


polarization  of  the  D’*p  in  S®  — » D’*p~  decays  provides  information  in 
tests  of  factorization.  The  polarization  of  D*^p”  iri  B~  —*  provides  for 

information  on  the  relative  sign  and  magnitude  of  o^/ot  [50|. 

The  longitudinal  polarization  measurements  tvere  obtained  by  fitting  the 
cosine  of  the  helicity  distributions  with  a functional  form  derived  from  the 
diSerential  decay  widths.  The  helicity  angle  (see  Section  4.6)  is  the  angle  one 
of  the  daugbtere  of  the  vector  pareot  makes  with  the  direction  of  the  vector 
particle  in  the  B's  cest-fi-ame  as  measiued  in  the  vector's  real-frame.  The 
differential  decay  width  expressed  as  a function  of  helicity  angles  Bjp*  and 
and  after  integrating  over  x is 

(6.3) 

where  x is  the  angle  between  the  normal  to  the  D‘  and  p decay  planes  and 
Ji±l  and  Ho  ate  the  transverse  and  loogitudinal  helicity  amplitudes  [46|. 

Completely  polarization,  in  the  longitudinal  direction,  would  result  in  a 
cos^B^j  distributions  for  both  the  D'  and  the  p~.  Complete  transverse  polar- 
ization would  result  in  sin-  8/,  distributions.  Ib  allow  for  both  possibilities  we 
fit  the  cos  8fy  distributions  to 

/^{“s*e/,  + |f(l-3cos2e*))  (6.4) 

where  N and  Ti/F  are  free  parameters  in  the  fit. 

Tb  measure  the  polarization  states  of  the  decay  products  of  the  B mesoo. 
the  cosine  of  helicity  distributions  were  formed  by  selecting  events  in  the  signal 
region  of  the  spectrum  defined  to  be  betvreen  5.272  GeV  and  5.286  GeV. 
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Figure  6.9  DeteclloQ  efficiencies  verses  cos(0;,). 

These  distributions  were  then  background  subtracted  and  effideocy  corrected 
to  insure  that  the  cosine  ofhelidty  distribution  was  represeuCative  of  Cruesignai 
events.  The  background  events  were  taken  from  the  Mgr;;  side  band  dehned 
to  be  between  5.2  GeV  and  5.26  GeV.  The  background  distribution  was  then 
scaled  and  subtracted  from  the  signal  distribution.  The  effidency  correction 
was  performed  to  compensate  for  the  differences  in  detector  acceptance  for 
final  states  which  occupy  different  regions  of  the  allowed  phase  space. 

To  correct  for  the  difference  in  detector  acceptances  we  measured  the  effi- 
ciencies as  a function  of  beiicity  angle  in  Monte  Carlo  generated  with  the  final 
states  uapolarized.  Figure  6.9  shows  the  effidencies  as  a function  of  cos  6^  for 
Che  three  types  of  final  states  considered.  For  the  the  efficiency  drop  in 


the  region  of  tbe  plot  near  cos  6/j  = (dot-daslied  line)  is  due  to  tbe  momen- 
tum dependent  detector  losses  for  low  momentum  tracks.  For  the  p",  the  drop 
in  the  efficiency  curve  (dashed  line)  in  the  region  near  cos  8/,  = can  also  be 
ascribed  to  detector  losses  for  charged  tracks  u'lth  Ion*  momentum.  The  drop 
in  the  p~  efficiency,  in  the  region  near  cos  0/,  = +1  is  due  to  the  photon  energy 
cut-off  which  reduces  the  acceptance  for  low  energy  ir®s.  Neutral  pions  from 
D*^  decays,  unlike  those  produced  in  p~~  decays,  are  produced  with  a relatively 
small  momentum  spread  which  explains  the  nearly  flat  efficiency  curve  (soUd 
line)  of  the  D"®. 

The  efficiency,  background  subtracted  helicity  distribution  was  fit  to  Equa- 
tion (6.4).  The  results  ate  shown  in  Figure  6.10  for  both  the  B~  -»  and 

the  t D’*p~  modes.  AU  three  Cfi  subdecays  are  combined  in  these  plots. 
The  measured  fraction  of  longitudlnaUy  polarized  p~  and  in  the  decay 
B~  -i  D'V  dec!^  are  Fi/F  = 86-7±6.7±4.4%  and  Fg/F  = 81-2±7.4±4.1% 
respectively.  In  tbe  D'*p~  mode  the  fraction  of  longitudinally  polar- 

ized p-  and  D'+  is  Fi/F  = 91.5  ± 4.9  ± 4.6%  and  Fg/F  = 88.2  ± 5.4  ± 4.4% 
respectively.  In  order  to  reduced  the  statistical  errors  a binned  two  dimen- 
sional fit  to  the  joint  (cosdjy-.cosSpj  distributions  was  performed.  This 
methods  yields  Fg/F  = 84.2  ± 5.1  ± 4.2%  for  the  B”  ->  mode  and 

F g/F  w 90.0  ± 3.7  ± 4.5%  for  the  — » D’'^p~  mode.  The  systematic  error 

was  estimated  by  varying  the  efficiency  correction  and  background  subtraction 

To  verify  the  efficiency  correction  and  background  subtraction  scheme,  the 
same  method  was  used  to  extract  the  polarization  of  the  D'  in  B”  -*  D'°ir" 
and  B®  -»  D'‘‘‘ir”  decays.  If  the  procedure  is  reliable  it  will  find  100%  of  the 


Figure  6.10  Plot*  of  polarizatioo  in  fl  -*  D‘*p  modes. 

Tbe  background  subtracted  efficiency  corrected  cos  Qi,  distributions  of  (a)  the 
and  (b)  tbe  p~,  in  B~  -*  decays.  In  (c)  the  ET'*'  and  in  (d)  tbe 

p~  distribution  in  the  -t  D'*p~  dec^. 


cos<8,^)  cos(,8g.) 

Figure  6.11  Plots  of  polarization  \a  B —*  D‘^p~  modes. 

The  background  subtracted  efficiency  corrected  cos  distributions  of  (a)  the 
i)*®and  (b)  D"'*’.  in  the  B~  — ♦ D’®tr”  and  — + D’^tr”  decay  modes,  respec- 

t's loogitudiDally  polarized  as  required  by  conservation  of  angular  momentum. 
The  fitted  distributions  are  shown  in  Figure  6.11.  The  fits  yield  r^/F  = 102.3± 
6.0%  for  the  B*  D'^tr-  and  Fi/F  = 99.8±4.6%  for  the  S®  -»  mode. 

Both  results  are  consistent  with  100%  polarization. 

6.3  Tests  of  Factorization 

In  Chapter  2 we  introduced  the  factorization  hypothesis  which  asserts  that 
amplitudes  in  decays  that  proceed  through  spectator  diagrams  can  be  factor- 
ized into  two  independent  matrix  elements,  each  governing  a separate  transi- 
tion. In  class  1 decays  the  amplitude  factorizes  into  a B — r i)  transition  and  a 
transition  from  the  vacuum  to  a light  meson.  To  test  the  factorization  hypoth- 
esis nonleptouic  class  I decays  were  compared  directly  to  semlleptonic  decays, 


Table  6.5 


Table  of  coostauts  used  for  test  of  [actoTlsatioo 


Constants 
1.12  ±0.10 
131.74  ±0.15  MeV 
224  ± 1 MeV 
298  ± 8 MeV 


Model 

BSW 

ISGW 

KS 


dT(B 


0.0023 

0.0020 

0.0024 


D'lu)ldq^  (GeV-7) 


at  the  appropriate  Similarly,  a comparison  between  tbc  polarization  in 
nonleptonic  class  I decays  and  scmileptonic  decays,  at  Exed  were  used  as  a 
further  test  of  factorization. 

6.3.1  Branchiae  Ratio  That 

If  factorizatioD  bolds  tbeo  the  efiective  Hamiltonian  described  in  Chapter  2 
con  be  u.sed  to  write  down  the  amplitude  for  a nonleptonic  decay  process  as  the 
product  of  two  independent  matrix  elements.  In  the  factorization  hypothesis 
the  amplitude  is  given  by 

A'>  = Bl(A-|d7«(l  - 75)u|0)(D-+|c7"(1  - tsWS®)  (6.5) 

where  the  operators  come  horn  the  factorized  current-current  effective  Hamil- 
tonian (Equation  (2.15)).  The  h~  here  represents  a light  meson,  either  a ir,p 


To  test  factorization  we  need  to  experimentally  verify  that  the  relation 


Aihyud\- 


(6.5) 


is  ssiisBeid.  Here  is  the  mass  squared  of  the  lepton-neutrino  s>‘stem  aud  is 
set  to  the  mass  of  the  of  the  light  mesons  to  insuie  tiiat  kinematics  properties 
in  both  systems  match. 

In  order  to  veri^  Equation  (S.G)  we  need  to  evaluate  the  expression  on  the 
right-hand  side  [Rih)-  The  ingredients  that  go  into  Rt>,  are  listed  in  Table  6.S. 
The  QCD  constant  23  can  be  obtained,  in  the  factorization  approximation, 
from  perturbative  QCD  (see  Equation  (2.16), (2.14)  and  (2.12)).  The  error  in 
Sj  refiecis  the  uncertainty  in  the  moss  scale  at  which  the  coefficients  Cj  and 
C2  were  evaluated.  The  value  of  = 0.975  ± 0.001  was  used. 

The  values  of  the  decay  constants  come  from  a variety  of  sources-  The 
pioD  decay  constant  is  obtained  from  pion  decay,  and  together  with  the  values 
for  the  Je  and  /«,,  are  listed  in  Table  6.5.  The  decay  constant  can  be 
determined  from  either  e^e"  /fi  or  from  r decays.  The  first  method  leads 
10  fp  s 215  ^ 4 MeV.  The  second  method  involves  extracting  the  value  of 
Jp  from  the  rate  T(r"  Vrp”).  The  value  of  /q,  can  be  similarly  obtained 
through  the  decay  process  r"  — t vroj".  A determination  of  both  fp  and 
has  been  recently  performed  in  Ref.  [51].  The  new  results  are  more  precise 
than  previous  values  and  are  used  in  this  analysis  [4Sj. 

To  evaluate  the  expression  on  the  right-hand  aide  of  Equation  (6.6).  iRexp)- 
we  need  to  find  the  diferential  branching  fraction  of  the  semileptonic  decay. 
To  do  so  we  must  interpolate  the  results  from  the  observed  diEerential  q^ 
distribution.  In  Figure  6.12  the  distribution  of  ^ is  shown. 

The  data  are  weighted  averages  of  CLEO  and  ARGIIS  results  and  the  fits  are 
to  the  three  theoretical  models  listed  in  the  plot  [52].  Since  the  distribution  is 
poorly  measured  In  the  q^  range  of  interest  (0  < < 1.6  GeV^)  we  must  rely 
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Figure  6.12  The  distribution  for  —r  U'*/ 


on  the  theoretical  models  to  determine  the  semileptonic  branching  fractions  at 
= 771^.  Trip  and  Tn^^,  The  spread  in  the  theoretical  models  is  small  so  we  use 
the  prediction  of  WSB  in  oixr  comparison.  The  values  obtained  at  each  are 
also  listed  in  Table  6.5. 

We  now  have  all  the  necessary  ingredients  to  perform  the  test  of  factoriza- 
tion. In  Table  6.6  the  comparison  of  the  dau  in  the  three  available  regions 
are  shown.  The  error  in  Rnp  is  tahen  solely  from  the  error  in  the  ^ 
branching  fraction.  The  error  in  Rn,  is  dominated  by  theoretical  uncertmnties 
in  Qj.  From  these  measurements  we  find  that  at  present  levels  of  precision 
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Table  6.6 

Table  of  constants  used  for  test  of  factoriaation 


region 

flfip  (GeV*) 

ft,h  (GeV2) 

B-^D*+x- 
8 -» D'*p- 

S-*D’+ar 

1.3  ±0.1  ±0.2 

3.4  ±0.3  ±0.5 
3.8  ±0.4  ±0.5 

1.2  ±0.2 
3.3  ±0.5 
3.0  ±0.5 

there  is  good  agreement  between  the  experimental  results  and  the  expectation 
from  factoriiation  in  the  range  0 < j-  < . 


As  suggested  by  Kornerand  Goldstein  [53],  more  subtle  test  of  factorization 
can  be  performed  by  comparing  the  final  state  polarization  in  hadronic  two 
body  decays  to  polarization  in  semlleptonic  decays  at  fixed  For  instance,  the 
ratio  of  longitudinal  to  transverse  polarization  (Vi/rx)  in  should 

equal  the  polarization  ratio  in  the  semlleptonic  decay  at  s 0.6  GeV^, 

D-+P-)  = ^(8“  (6.7) 

This  method  has  the  added  advantage  that  the  QCD  coefficient  ai  does  not 
enter  the  cnmpatison  [45|.  This  coefficient  contributes  a sizeable  uncertunty 
to  the  branching  fraction  test  as  shown  in  the  previous  section. 

There  are  simple  physical  arguments  which  describe  the  polarization  of  the 
final  states  in  semileptouic  decays.  At  low  rp  the  longitudinal  polarization 
dominates.  Here  the  IT  has  the  maxknunj  possible  momentum  and  the  lepton 
neutrino  pair  travel  collinearly  in  the  opposite  directiou  with  their  spins  aligned 
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Figure  6.13  Polariiation  componeals  of  distribution  —* 

to  conserve  aaguJar  momeutunj.  This  causes  the  D”  to  be  loogitudlaally  po- 
larized. As  7^  — t 7^^  we  expect  that  ail  three  polarizations  are  equally  likely 
since  the  D*  is  almost  at  rest  and  its  small  residua]  momentum  will  equaUy 
populate  the  three  polarization  states.  Here  we  expect  r£/rj  = 1/2  and  thus 
transverse  polarization  dominates  as  7^  — » 7^^. 

Fbom  these  simple  arguments  Rosner  predicts  86%  longitudinal  polarization 
for  S'’  -*  D’*ri/  at  7^  = mp  (54].  Neubert  predicts  similar  results  [55). 
Figure  6.13  shows  Neubert's  results  for  the  production  of  transversely  and 
longitudinally  polarized  D"*"  in  decays  for  the  whole  7^  range. 

The  dashed  solid  curve  is  the  6t  to  data,  the  dashed  curve  is  the  component 
due  to  transversely  polarized  D“s  and  the  dot-dashed  curve  is  the  prediction 


for  longitudinally  polarized  D‘s.  The  data  are  weighted  averages  of  CLEO  and 


l^ble  6.7 


Table  of  ratios  of  color-suppressed  braacbing  fractioas. 


Modes 

Ratio  of  B 

0(0“  -»  A-V)/B{D«  -►  /f-rr-*-) 

0.57±0.13 

8(D°  A—rr-^) 

0.47*0.23 

B{D®  -t  X-X+) 

0.77  ± 0.25 

B(D*  -»  A’'''A'+)/B(D,  ->  «Sir+) 

0.95*0.10 

B(D*  -*  K<>K*)/B(D,  -*  <>ir+) 

1.01*0.16 

ARGUS  results.  The  solid  line,  at  = mj,  gives  F^/r  = 85%  which  is  in 
good  agreeiDeiit  with  the  measured  result. 

ri/r  = 90±3.7±4.5% 

from  Section  6.2.  This  also  supports  the  factorization  hypothesis  since  the 
comparison  reveals  similar  results  for  both  the  nonleptonic  and  semileptonic 
decays. 


6.4  Determination  lui  l.la*rl  and  the  Relative  Sign  of  u^/oi 
Importam  infortuation  on  factorization  and  the  contribution  to  the  am- 
plitude from  non-factorizable  terms  can  be  obtained  by  finding  the  values  of 
the  parameters  ni  and  02.  In  the  factorization  hypothesis  the  amplitude  for  a 
particular  spectator  decay  process  is  written  in  terms  of  these  parameters  and 
factorized  hadronic  matrix  elements.  The  parameters  ax  and  02  are  defined  as 

where  i — and  O]  and  02  are  related  to  the  scale  dependent  Wilson  coeS- 
cienis  calculable  from  perturbative  QCD. 
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In  cbarm  data  the  factorizaCiou  approach  outlined  in  Chapter  2 does  not 
explain  the  bulk  of  exclusive  non-leptonic  decays.  For  example,  the  ratio  of 
color-suppressed  to  non  color-suppressed  branching  fractions  listed  in  Table  6.7 
are  in  dramatic  disagreement  with  the  predicted  value  [66].  Also  the  results 
obtained  from  a fit  to  a iarge  sample  of  charmed  decays  gives  a value  of  ( = 0 
which  implies  Nc  = co.  This  corresponds  to  dropping  the  contribution  from 
the  color  mismatched  4-quark  operator  altogether  and  seems  co  suggest  that 
the  non-factorisable  contributions  cancel  the  contributions  from  the 

i-Cs(DAIOilB) 

term.  In  charm  decays,  however,  the  situation  is  complicated  by  final  state 
intoractioDS  ^ce  the  low  momentum  transfers  allow  for  the  possibility  of  final 
state  rescattering  before  the  hadrons  have  moved  significantly  far  away  from  the 
strong  interaction  region.  Measurements  of  these  parameters  in  the  B system, 
where  final  state  interaction  ate  believed  to  play  a less  significant  role  could 
shed  some  light  as  Co  the  role  played  by  the  nou-faccorisabie  contributions  to 
heavy  quark  decays  [67]. 


Tb  determine  the  the  values  of  the  parameter  ai<  the  magnitude  and  sign 
ofnz/ar  we  use  the  branching  fraction  measurements  in  Tables  6.1,  6.2  and  the 
theoretical  predictions  from  Neubert  el  al.  and  Deandiea  et  of.  [10,23].  The 
theoretical  predictions,  listed  in  Table  6.8,  have  been  rescaled  co  accommodate 
common  values  of  /p  = 220  MeV,  1]^  = 0.041  and  rg  s 1.44  ps  [68],  In  the 
Neubert  et  of,  approach  the  heavy  quark  symmetries,  with  mass  corrections. 


employed  in  determining  the  form  factors  in  heavy  to  heavy  (B  —*  D) 


'Kible  6.8 


Table  of  theoretical  predictions  for  B branching  fractions 


B Mode 

5,fc  (%)  BSW  11 

6, ft  (%)  CDDFGN 

D+ir- 

0.264o2 

0.278af 

D*p- 

0.621«2 

0.717a? 

D’  + 7T- 

0.2540= 

0.278a? 

D’*p- 

0.702OJ 

0.949a? 

Efi-r:- 

0.265(oi  1.250o2)= 

0-278(oi  + l. 12702)= 

[fip- 

0.662(oi  + 0.662o2)= 

0.717(oi+0.458o2)= 

D'®7r- 

0.255(o,  1.292o2)= 

0.278(01  + 1.52402)= 

D‘°p~ 

0-703(a?  + 0.635a?  + 1.487ai02)  0.949(o?  + 0.53o?  + I.31ai02) 

DOffO 

0.201a| 

0.502a? 

D»n 

0.147a? 

rflr/ 

0.030a? 

D®/ 

0.136a| 

0.077a? 

D®u! 

0.264a? 

0.077n? 

D*®n® 

0.201a? 

0.405a? 

D*®>) 

0.264a? 

0.097n| 

D'%' 

0.030a? 

0.136a? 

0.270a? 

D*®u; 

0.264a? 

0.270a? 
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transitions.  For  heavy  to  light  (fl  -*  /i]  transitions,  they  use  a slightly  modified 
version  of  the  original  BSW  model  where  the  form  factors  Ff,  aud  4]  remain 
monopole  types  and  the  F\,  .4o,  and  V form  factors  have  dipole  type  g~ 
dependence.  Dcaodrea  el  al.  use  exact  heavy  quark  symmetry  to  extract  the 
form  factors  in  heavy  to  heavy  transitions  by  extrapolating  the  Isgur-Wlse 
function  horn  the  symmetry  point  using  an  unproved  form  of  the  relativistic 
oscillator  model.  The  heavy  to  light  form  factors  are  estimated  by  use  of 
heavy  quark  and  chiral  symmetries  with  mass  corrections  and  monopole  type 
dependence  for  all  form  factors.  These  models  will  henceforth  be  referred 
to  as  BSW  II  and  CDDFGN  respectively  (57,50). 

In  clans  I decays,  the  branching  fractions  are  proportional  to  the  BSW 
parameter  s^.  By  performing  a leant  squares  fit  of  the  branching  fractions, 
B®  -»  D*p~  and  D'*p~  to  the  theoretical  prediction  of  BSW 

II  and  CDDFGN  we  obtained 

|tti!  =1.18±  0,025  ± 0.023±  0.198  BSWI! 

(6.8) 

(oil  =1.09  ± 0.023  ± 0.021  ± 0.183  CDDFGN. 

The  first  error  Is  statistical,  the  second  error  is  the  systenutic  error  obtained 
in  the  least  sqaures  fit  and  includes  both  branching  fraction  systematic  errors, 
the  third  error  is  due  to  the  uncertainty  in  the  B meson  production  fractions 
aud  lifetimes 

= 1.14  ±0.14  ±0.13(15). 

The  theoretical  uncertainties  were  not  included  in  the  systematic  enore. 
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The  relative  sign  and  value  of  a2/ui  can  be  extracted  fay  comparing  daae  lit 
decays  to  class  1 decays  by  formiag  the  following  ratios  of  brancbing  fractions 


The  numerical  constants  and  multiplying  the  aj/ui  terms  are  functions 
of  the  various  form  factors  and  decays  constants  which  enter  in  the  covariant 
decomposition  of  the  amplitudes  (see  Chapter  2).  The  theoretical  predictions 
from  the  BSW  II  and  CDDFGN  models  are  also  listed  in  Table  6.8.  The  values 
were  used  in  separate  least  squares  fits  to  Equations  (6.9)  and  give 


In  the  fit  we  ignored  the  theoretical  uncertainties  in  ili  through  /24.  The 
systematic  errors  were  first  estimated  for  each  ratio  separately  and  then  used 
as  weights  in  the  least  squares  fit.  The  results  both  obtain  a positive  value  for 
the  relative  sign  of  02/01  and  are  in  good  agreement  in  maguitude. 

The  model  predictions  for  the  ration  ihi  are  believed  to  be  unrebable  [58]. 
To  account  for  this  we  performed  the  fit  a second  time  using  only  the  ratio 
Rl  - flj.  The  results, 

— (Rl  - fij)  = +0.13  ± 0.034  ± 0.081  ± 0.022  BSWII 

ai  (6-11) 

— (Rl  -R3)  = +0.15  ±0.032±0.057±  0.023  CDDFGN 


(6.9) 


— = + O.U±0.032±D.043±0.025  BSWII 
“1 

— = + 0.16±0.031±  0.041  ±0.027  CDDFGN 


(6.10) 


1S3 

do  not  differ  signiffcanlly  from  Cbe  result  obtaioed  with  all  four  ratios. 

The  fits  to  data  described  above  do  oot  include  theoretical  uncertainties 
in  calculating  the  form  factors  which  are  quite  large  (on  the  order  of  o0%  in 
some  modes).  Abo.  the  expected  variation  between  naive  factorization  and  the 
rule  for  discarding  the  1/Nc  term  in  class  I amplitudes  may  not  apply  for  class 
II  amplitudes.  These  uncertainties  renders  the  global  fit  procedure  somewhat 
questionable  since  it  relies  on  ratios  of  class  I decays  to  claas  III  decays  cbe 
latter  containing  both  types  of  amplitudes.  Tb  examine  Che  rebtive  sign  of 
02  and  O)  Gourdin  el  al.ase  a method  which  relies  on  class  III  decays  only 
[37].  Their  method  consists  of  plotting  the  allowed  region  io  the  ni,  03  plane 
and  then  look  for  a noo*empty  intersection  of  the  allowed  domain.  In  their 
paper  they  compare  three  theoretical  modeb,  including  BSW  I]  and  CDDFGN. 
already  discussed  and  the  original  BSW  model  which  they  refer  to  as  BSW  L 
Their  plots  were  recreated  with  the  branching  fraction  results  and  numerical 
predictions  listed  in  Thble  6.S  and  in  Ref.  [57]. 

The  plots  were  made  by  plotting  the  lo  standard  deviation  results  for  the 
branching  ratios  of  class  III  modes  and  The  union 
of  all  four  bands  gives  Cbe  allowed  region  for  the  values  of  ox  and  03.  Both 
the  BSW  II  and  CDDFGN  modeb  have  common  Intersection  domains  which 
intersect  in  the  positive  region  of  03  and  ax . 

Tb  further  consCrain  Che  allowed  region  of  ax  and  03  we  superimposed  the 
ax  bands  from  class  I decays.  The  theoretical  prediction  for  the  three  modeb 
considered  are  abo  obtained  from  Ref.  [57]  and  from  the  values  in  Thble  6.8. 
The  union  of  all  bands  gives  a graphical  indication  of  how  well  Che  data  are  fit 


by  cbe  cbeorecicai  modeb. 
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Table  6.9 


6.4.2  TlpterminatiDn  FVom  B~  — * D"^n  PnlAriy.atinn 

Recently,  Keutn  has  pointed  to  another  method  of  detcrminiiig  the  relatii 
sign  of  02/01  [60].  His  method  consists  of  relating  the  polarizacion  state  i 
class  III  vector-vector  decays  to  the  decay  amplitude  which  depends  on  bol 
oi  and  02.  In  particular,  the  analysis  uses  the  B~  -*  decay  mode. 

The  lougitudinal  polarisation  is  related  to  the  amplitudes  by 


Pi  = 


(6.12) 


where  ^ w 02/ai,  .4^^  and  are  the  amplitudes  for  the  color-matched  and 
color-mixed  amplitudes  in  a particular  helicity  state.  FVom  the  expression 
above  several  prediction  for  the  longitudinal  polarizatiou  were  made.  In  ob- 
taining the  predications,  factorization  is  assumed,  the  effects  from  final  state 
interactions  are  ignored  and  the  calculations  were  done  using  both  zero  and 
finite  width  approximation  for  the  p~  maas.  The  results  obtained,  in  the  finite 


p width  approximation,  are  listed  In  Table  6.9. 


Figure  6.15  The  lo 


g-  D’^p~  decays. 


The  amplitudes  used  in  evaJuating  Equation  (6.12)  are  computed  with 
hadronic  form  (actors  as  derived  in  several  theoretical  models.  In  calculating 
the  heavy  to  heavy  5 — » D*  form  factois.  the  set  of  predictions  labeled  HQET 1 
(as  in  CDDFGN)  are  associated  with  exact  heavy  quark  symmetry.  The  HQET 
U predictions  include  mass  corrections  (as  in  BSW II).  In  computing  the  heavy 
to  light  B-*  p transitions.  BSW  11  and  the  CDDFGN  are  the  models  already 
described.  In  the  P6  model,  monopole  and  constant  terms  are  assumed  for 
the  dependence  of  the  factors  factors  [50]. 

The  measurement  of  the  longitudinal  polarization  in  fl’  — t D"®p”  was 
discussed  in  Section  0.2.1  and  the  value  obtained  was 
Pi  = ^ = S4.2  ± 5.08  ± 4.2%. 


This  value  is  compared  with  the  theoretical  predicatioo  in  Figure  6.15.  The 
bauds  represent  the  minimum  and  maximum  predictions  taken  from  Table  6.S 
for  positive  and  negative  ^ = 0.20  values.  Our  measurement  is  not  sul&ciently 
precise  to  distiugnish  between  positive  and  negative  vaiues  of  C>  The  error  bars 
represents  the  Iff  statistical  error. 


Two  types  of  weak  decay  diagranrs  dominate  the  decay  processes  of  heavy 
Savor  mesons.  The  external  diagram  mediates  class  I decays  and  is  described  by 
the  production  of  a Ught  meson  from  the  hadronlaation  of  quark  pairs  produced 
at  the  W vertex.  The  internal  diagram  mediates  class  U decays  and  the  quarks 
produced  at  the  W vertex  badronize  into  separate  mesons  which  imposes  color 
matching  requirements  on  the  color  fields.  The  color  matching  requirement 
implies  chat  class  II  decays  are  suppressed  relative  to  class  I decays.  From 
simple  color  counting  we  expect  a suppression  rate  of 


where  the  additional  factor  of  1/2  comes  from  the  wavefunction. 

The  naive  expectation  is  modified  by  hard  gluon  corrections.  U we  assume 
strict  factorization  we  can  ignore  the  non-fsctorizable  terms  In  the  effective 
Hamiltonians.  Equations  [ i,i.ef!Halii]  and  ( i,i,elfHa2ii).  and  write  the  ratios 
of  partial  widths  in  terms  of  matrix  elemenu  and  QCD  coefficients  s]  and  02 


r(fi®  1 

r(B«-*D+x-)  “ 2 


(6.13) 


If  we  use  the  LLA  results 


T^ble  6.10 


Tablp  of  ratio  of  class  II  to  class  I B decays 


Ratio  of  bcanchisg  fractious 

UL  (90%  C.L.)  BSWII 

CDDFGN 

< 0.11 

0.0129 

0.0406 

5{S‘>-»dO>j)/B(SO-.D*x-) 

< 0.11 

0.0119 

S(S<>  -»  D+ir-) 

<0.09 

0.0019 

B(S®  ->dV)/B(B°  -*  D*p-) 

<0.07 

0.0037 

0.0024 

B(B“  ^ D“w)/B(B®  -»  D*p-) 

<0.07 

0.0072 

0.0024 

6(8“  -t  D*+x-) 

<0.18 

0.0164 

0.0328 

S(B°  ^ D'*n-) 

<0.16 

0.0176 

0.0079 

B(B"  -t  D*V)/B(B‘'  -*  D'+ir-) 

<0.« 

0.0020 

<0.14 

0.0033 

0.0064 

B{B<>  D'“ai)/B(fi‘’  -*  D'*p~) 

<0.14 

0.0064 

0.0064 

Note:  Tbe  values  of  aj/oi  used  are  from  Equation  (G.IO) 


we  got  a suppression  factor  of  1/200.  This  results  is  extremely  sensitive  to  tbe 
precise  values  of  the  QCD  coefficients  due  to  the  cancelation  of  tbe  coefficients 
in  the  nunicraCor.  For  example,  the  NNLA  values  gives  a suppression  factor 
of  10“^.  Other  factors  such  as  tbe  precise  value  of  tbe  fbriti  factors  and  the 
effects  of  final-state  interactions  influence  the  suppression  although  these  tend 
to  be  small  as  can  be  seen  by  inspecting  the  theoretical  prediction  in  Table  6.6. 

We  set  upper  limits  on  color-suppression  by  taking  ratios  of  class  II  to  class 
I decays.  In  Table  6.10  Che  upper  limits  at  tbe  60%  confidence  level  are  shown. 
The  comparison  is  made  between  modes  with  siraibu-  particle  content  to  reduce 
the  systematic  errors,  some  of  which  cancel  in  the  ratio.  Also  listed  in  Table 
6.1Q  are  the  predictions  of  Neubert  et  of.  and  Deandrea  e(  ol.  where  tbe  ratios 


180 

of  the  BSW  parameters  03/01  from  Equation  (6,111  were  used  instead  of  the 
QCD  parameters  used  in  the  estimate  abote.  The  upper  limits  listed  in  Table 
C.lOare  consistent  with  coior-suppresslon,  however  the  largo  uncertainties  in 
both  the  experimental  parameters  03/01  and  the  theoretical  models  makes  it 
difficult  to  estimate  the  expected  rate.  Both  models  are  consistent  with  our 


CHAPTER? 

SUMMARY  AND  CONCLUSIONS 
Using  2,040  pb"'  of  data  collected  with  the  CLEO  II  detector,  we  haw 
reconstructed  twenty-two  decay  modes  of  the  B meson,  exclusively,  in  three 
efi  and  one  D*  decay  mode.  The  large  sampled  allowed  the  measurement  of 
branching  fractions  for  twelve  of  these  inodes  to  an  accuracy  limited,  in  some 
inodes,  only  by  the  systematic  errois  associated  with  the  branching  fraction 
measurements  on  charm  decays.  Two  6 decay  mode  branching  fractions,  the 
B~  -»  D®oj"  and  the  -»  C'aJ"  were  measured  for  the  first  time. 

Nine  of  the  twelve  non  color-suppressed  decay  modes  measured  contained 
a Cf  in  the  decay  chain.  For  each  of  these,  a new  measurement  was  provided 
by  the  —*  K"Tt'*’5T^  subdecay  mode  for  the  first  time.  All  of  the  branchiog 
fractiou  measurements  are  listed  in  Tables  6.1  through  6.3.  They  form  the 
most  precise  measurements  so  far  for  this  important  class  of  decays. 

The  precision  measurements  were  used  to  test  the  factorization  hypothesis 
by  comparing  class  I decays,  fl®  -t  D'+fi",  to  semileptonic  decays  at  appropri- 

90.2  ± 3.7  ± 4.5%  of  the  -»  D'*p~  decays  were  found  to  be  longitudinally 

polarized  also  in  agreement  with  the  factorization  hypothesis. 

Tb  branching  fraction  measurements  were  also  used  to  determine  the  BSW 
parameters  [oil  and  Che  magnitude  and  sigu  of02/ax.  The  values  were  obtained 
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by  Buini  to  the  bruchiog  frartions 


tbeoretical  models-  These  sulues 


are  listed  below, 


BSW  II 

|ai|  1.18  ± 0.025  ± 0.023  ± 0.198 
02/01  +0.15  ±0.032±  0.043  ±0.025 


CDDFGN 

1.09  ±0.023  ±0-021  ±0-183 
+0.16  ± 0.031  ± 0.041  ± 0-027 


The  positive  value  for  os/tti  dlBet  from  Che  expectation  extrapolated  from 
measurements  performed  on  charmed  decays.  In  charm  decays  the  Bt  to  data 
imply  that  Che  noo'factocirable  terms  contribute  in  such  a way  as  to  cancel 
the  effect  from  the  1/jV;  term.  Extrapolating  this  to  6 decays  we  expect  a 
negative  value  of  02/01  which  is  not  obtained  by  the  fit  to  data. 

The  large  sample  of  events  was  also  u-sed  to  search  for  tec  color- 
suppressed  decays  of  the  B mesons  to  a single  charmed  plus  a light  neutral 
meson.  The  modes  examined  were 

8“  ePii°  8“  8°  -*  D°ti'  8®  -»  8®  -*  D®+ 

8®-*D’®ir®  fl®-.D'®ti  8®-*D‘®p®  8®->D‘®+. 

While  no  signals  were  obtained  nut  results  showed  that  color-suppression  is 
active  in  B mesons  according  Co  theoretical  expectations.  The  upper  limits  set 
are  about  a factor  of  20  higher  than  Che  expectation  from  naive  factorization. 

In  Che  near  future,  the  large  data  sample  will  be  used  to  further  investigate 
the  properties  of  the  8 mesons.  This  includes  additionaJ  modes  such  Cabibbo 
suppressed  decays,  double  D decays  and  perhaps  the  addition  of  other  D me- 
son subchannels.  Currently,  none  of  the  colot-suppressed  modes  have  beeu 
observed  but  an  IntTiguing  eubancement  is  seen  in  one  of  the  D®  submodes  in 
the  8®  — * D®x®  channel.  This  mode  is  curreotly  under  investigation. 


APPENDIX  A 

MEASUREMENT  OF  IMPACT 
PARAMETER  SPREAD 


Introduction 

The  sipresd  in  Che  impaoi  parameter  varies  as  a function  of  both  momen- 
tum and  azimuthal  angle  0 (FICD).  The  impact  parameter  (DBCD]  of  a track 
is  defined  as  its  point  of  closest  approach  to  the  beam  spot  (see  Figure  4.4). 
The  sign  of  DBCD  gives  the  direction  of  the  track  with  respect  to  the  displace- 
ment from  the  beam  spot.  Usually.  In  CLEO  analyses  a five  or  six  millimeter 
cut  is  placed  on  a track's  impact  parameter,  independent  of  its  momentum  or 
azimuthal  angle.  We  can  improve  this  selection  criteria  I7  Impo^ng  a DBCD 
cut  which  depends  on  both  the  track  momentum  and  angle  since  the  spread 
in  impact  parameter  varies  considerably  over  both  these  quantities.  In  order 

sure  the  dependence  of  the  spread  in  data.  !n  this  appendix  we  describe  the 
analysis  procedure  used  to  measure  the  spread  and  also  discuss  the  shift  in 
the  mean  of  DBCD  for  tracks  which  is  also  a function  of  momentum.  After 
the  measurements  were  performed  a set  FORTRAN  routines  were  created  that 
return  the  value  of  the  width  of  the  spread  and  the  shift  in  the  mean  of  DBCD 
as  a function  of  momentum  and  Crack  azimuthal  angle  «.  These  routines  were 
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UKd  ill  tb^  B analysis  to  select  tracks  trom  the  origin  in  lieu  of  the  standard 

The  moTnentum  and  0 dependence  of  the  spread  in  the  impact  parameter 
can  be  expressed  as 

‘'Leo  (f'  = ‘'t(p)  + ^ + ‘'w  <’■  (A-'l 

Here  the  angular  dependence  of  the  overall  DBCD  spread  is  assigned  to  the 
beam  spread  with  the  terms  and  the  momentum  dependence  is  assigned 

to  the  term.  All  three  eomponents,  added  in  quadraturei  describe  the 
functional  form  of  the  overall  spread. 

The  momentum  dependence  of  the  Impact  parameter  is  dominated 
by  multiple  scattering  for  tracks  that  ori^nate  at  the  primary  vertex.  As  a 
particle  traverses  detector  material,  its  trajectory  is  deSected  by  the  influence 
of  the  nuclei  it  encounters.  The  sum  of  these  deflections,  for  small  angles, 
follows  roughly  a Gaussian  distribution  where  the  rout-mean-square  deflection 
is  given  by 

, 2iMev  rr 

here  Srms  Is  deflection  angle  from  the  normal  to  the  plane  of  incidence,  A'q  Is 
the  radiation  lengths  and  t is  the  thickness  of  the  media  being  traversed.  The 
momentum  dependence  nf  the  r.m.s  deflection  thus  influences  the  momentum 
dependence  of  the  impact  parameter  spread  since  the  spread  is  the  r.m.s  dis- 
tance between  the  track's  point  of  closest  approach  and  the  beam  spot.  Other 
factors,  such  as  poorly  reconstructed  tracks  and  non-primary  vertex  decays 


also  contribute  to  the  impacr  parameter  spread.  To  allow  for  these  factors  the 


momeatum  dep«ndeDcc  of  the  spread  was  not  taken  directly  from  the  multiple 
scattering  form,  Instead  the  function 

<’TiP)  = -^  + Bp  + C (A-3) 

The  functional  form  of  Equation  (A.3)  was  found  to  adequately  fit  the  data. 

A.l  Data  Samnle  and  Analysis  Procedure 
The  data  used  consisted  of  4s2  through  4s8  hadronic  sample  used  in  the  B 
meson  analysis.  To  limit  the  number  of  events,  only  1000  events  were  taken 
froia  each  run.  If  ack  selection  requirements  were  also  kept  minimal.  The  only 
selection  requirement  imposed  was  that  the  tracks  hit  at  least  10  layers  io  the 
drift  chamber.  In  addition  a 40.000  B9  event  sample  of  production  Monte 
Carlo  events  used  to  verify  chat  the  Monte  Carlo  simulated  the  data  well. 

The  momentum  dependence  of  the  spread  in  DBCD  was  determined  by 
measuring  the  standard  deviation  of  DBCD  as  a function  of  momentum  at  a 
fixed  azimuthal  angle  |cos$|  < 0.96.  The  distributions  for  a selection  of  the 
4s2  data  are  shown  in  Figure  A.l  [69].  The  distributions  were  fitted  to  a dou- 
ble Gaussian  plus  a Cbebychev  polynomial  where  the  width  of  the  secondary 
Gaussian  was  constrained  with  respect  to  the  primary.  The  standard  devia- 
tions for  each  of  the  DBCD  distributions  were  then  plotted  as  a function  of 
momenta  and  fit  to  Equation  (A-3).  As  a coneisCeocy  check  the  4s2  impact 
parameter  spreads  were  refit  with  a single  Gaussian  plus  a higher  order  Cheby- 
chev  polynomial.  The  widths  obtained  did  not  change  significantly  but  the 


The  impact  parameter  spread  verses  momentum  distributions  determined 
by  the  double  Gaussian  fits  are  shown  in  Figure  4.5  for  each  of  the  7 data 
sets  plus  the  432-433  BB  Monte  Carlo.  The  value  of  the  fit  parameters  can 
be  (bund  in  Ref.  (59).  To  insure  that  the  momentum  dependence  was  well 
modeled  for  low  momentum  tracks  the  lowest  momentum  bin  in  the  4s2  sample 
was  divided  into  20  MeV  bins  centered  at  70  MeV.  90  MeV.  110  MeV  and  170 
MeV.  The  DBCD  spread  verses  momentum  distributions  were  then  refit  to 
Equation  (A.3).  The  results  obtained  did  not  change  appreciably  from  the 
large  low  momentum  bio  results.  Figure  4.5  shows  the  variation  in  the  spread 
of  from  data  set  to  data  set.  The  maximum  variation  is  approximately  60 
microns  at  1.0  GeV  between  the  4s5  and  4s3  data  samples.  The  variations 
between  data  sets  ate  automaticly  included  in  the  FORTRAN  routines. 

To  measure  the  azimuthal  dependence  of  the  impact  parameter  spread 
the  momentum  was  fixed  between  1.0  GeV  and  1.2  GeV.  This  momentum 
interval  corresponds  to  a slowly  varying  region  of  the  momentum  spectrum. 
The  standard  deviation  of  the  spread  was  then  measured  in  cosp  bios  0.1 
units  wide.  The  combined  4s2-4s4  plus  the  Monte  Carlo  data  were  fit  as  before 
and  the  standard  deviations  were  plotted  as  a function  of  cos  0.  The  two  plots. 
Figure  A.2  (a)  and  (b)  were  fit  to  the  square-root  of  equation  (A.l) 

I’oaco  = \/<r?  + iri,(l  - cos^b)  -I-  cos^  0 (A.4) 

with  the  value  of  Oy,  fixed  to  the  known  beam  spread  in  p.*  The  value  of  the 
spread  in  x.  the  fit  parameter  ax,,  is  the  z component  of  the  overall  spread 

f The  30  micron  value  used  in  the  fit  to  Monte  Carlo  data  is  the  default 
QQ  value.  The  5.6  mlcrou  value  used  in  the  fit  to  data  is  the  CESR  beam 
dimension  in  y divided  by  V2. 
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COS(  ®dbcd) 


Figure  A.2  The  azimuihal  dependeuce  of  the  spread  in  DBCD. 
The  fitted  distribution  in  (a)  Monte  Carlo  and  in  (b)  4s2<4s4  data. 
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due  to  tbe  spread  of  the  colliding  beams  in  tbe  horizontal  direction  (<r>|).  The 
value  obtained  was 

!ti.  = 336.1±11.3±14.0. 

As  a check,  the  horizontal  beam  spread  obtained  from  Monte  Carlo  was  com- 
pared with  the  known  QQ  default  value  of  SCO  microns.  The  value  obtained 
from  the  fit  to  Monte  Carlo  data  was  520.9  microns.  The  difference  between 
the  fit  result  and  tbe  default  QQ  value  was  used  as  the  systematic  error. 

A.2  Shift  in  Mean  of  DBCD 

Mike  Zoeller  performed  a similar  analysis  on  the  momentum  dependence  of 
DBCD  [60].  Ke  found  a track  sign  and  momentum  dependent  shift  of  the  mean 
of  DBCD  which  was  clearly  visible  for  low  momentum  tracks.  In  Figure  A.3 
(a)  and  (b)  tbe  DBCD  distributions  for  po^tive  and  oegative  low  momentum 
tracks,  in  4s2  data,  are  shown  together  with  the  Monte  Carlo  predictions  for 
tbe  momentum  dependence  of  the  shift.  Figure  A.3  (c). 

The  track-sign  dependence  of  the  shifts  can  be  understood  if  wc  both  the 
sign  convention  of  DBCD  and  the  effect  of  energy  Jesses  on  charged  track 
curvatures.  The  energy  loss  effect  Is  shown  in  Figure  A.4.  Without  energy 
losses,  both  of  tbe  tracks  in  Figure  A.4  have  equal  but  opposite  values  of 
DBCD.  Energy  losses  increase  the  curvature  of  the  tracks  affecting  tbe  track 
on  the  left-hand  differently  than  the  track  on  tbe  right.  The  track  on  the 
left  side  of  the  beam-spot,  which  b defined  to  be  positive,  moves  towards  the 
beam-spot.  The  track  whose  point  of  closest  approach  lies  on  tbe  other  side 
is  sbifted  away  from  the  beam-spot.  Thus  tracks  with  the  same  sign  but  with 


Figure  A. 3 Shift  ie  mrae  of  DBCD  for  pceitive  and  negative  tracks. 

The  DBCD  distributions  for  (a)  positive  and  (b)  negative  low  momentum  pions. 
In  (e)  the  Monte  Carlo  predictions  for  the  momentum  dependence  of  the  shift 
in  the  mean  of  DBCD  for  low  toomentum  <r±  [60j. 


Figure  A.4  The  effect  of  energy  losses  on  primary  tracks. 

points  of  closest  approach  on  opposite  sides  of  the  beam-spot  are  shifted  from 
zero  in  one  direction  or  the  other  depending  on  the  sign  of  the  track. 

The  shifts  in  the  mean  were  measured  in  4s2  data  using  a sample  of  low 

momentum  tracks  can  curl  multiple  times  in  the  large  drift  chamber  volume, 
the  possibility  exists  that  the  tracking  software  will  pick  up  the  same  track 
more  than  once,  each  time  assigning  it  an  opposite  charge.  These  multiply 

distributions  of  the  spread  in  DBCD  (see  Figure  A.3  (c)).  We  verified  this 
aria  where  similar  backgrounds  were  eliminated  by  requiring  that 


in  Monte  Ca 


th&t  che  charge  of  the  recooslructed  track  match  the  charge  of  the  generated 


APPENDIX  B 
EFFICIENCY  CHECKS 


Various  studies  were  performed  to  verify  the  accuracy  of  the  Monte  Cario 

reconstruction  procedure.  The  accuracy  of  the  Monte  Simulations  is  important 
in  extracting  absolute  branching  fractions  since  Monte  Carlo  efSciencies  were 
used  to  correct  the  event  yields.  The  tracking  errors  used  in  the  systematic 
error  estimates  were  determined  in  these  studies.  In  this  appendix  we  outline 
the  analyses  used  to  verified  the  Monte  Carlo  accuracy  with  references  provided 
to  the  work  performed  by  various  members  of  the  CLEO  collaboration.  Also, 
we  tested  the  rcsponce  of  the  Monte  Carlo  effeciencies  and  under  the  particular 
particle  ID  consistency  requirement  used  in  the  B analysis.  Finally,  we  present 
Che  analysis  used  to  select  the  momentum  and  azimuthal  dependent  impact 
parameter  cut,  testing  the  performance  and  verifying  the  response  of  the  Monte 
Carlo. 

B.l  Monte  Carlo  Verification  StiiHips 
The  accuracy  of  detection  efficiencies  for  charged  and  neutral  tracks  were 
checked  in  several  ways,  by  various  members  of  the  CLEO  collaboration.  The 
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results  are  used  in  all  CLEO  analyses  as  an  esiiiuace  uf  the  accuraty  of  tbe 
Monte  Carlo  efficiencies.  Work  on  reducing  these  tracking  errors  is  an  ongoing 
and  important  process  since  some  analysis  are  already  dominated  by  systematic 
errors.  In  fact,  this  is  case  for  some  of  the  higher  multiplicity  modes  in  the  B 
analysis. 

The  uncertainty  in  detection  efficiencies  for  tracks  above  225  MeV  was  de- 
termined by  measuring  tbe  probability  of  not  finding  two  tracks  in  radiative 

This  is  then  repeated  in  Monte  Carlo  and  the  results  are  compared  [62].  For 
tracks  below  225  MeV  the  expectation  that  tbe  decety  angle  distributions  in 
Dd)r+  decays  are  symmetric  about  aero  wa.s  used  to  estimate  the  un- 
certainty in  detection  efficiency-  The  relative  corrections  between  Monte  Carlo 
and  data  were  compared  to  arrive  at  the  5.0%  tracking  error  used  it  the  fl 


to  the  average  PDG  value.  The  combined  systematic  and  statistical  errors, 
together  with  the  assumption  that  the  ^ulation  of  ir°  -t  77  and  i)  77  are 
similar  enough  so  that  the  ratio  of  their  efficiencies  cancel,  gives  an  uncertainty 
of  5.0%  in  detection  efficiency  for  single  neutral  pions  |62).  Several  ocher  con- 
sistency checks  including  comparisons  between  full  and  partially  reconstructed 
D®  -*  fC'tt+ff'’  event,  other  decay  angle  distributions  tests,  other  ratio  of 
branching  fraction  measurements  and  comparison  of  inclusive  D*"*"  with 

ciencies.  A synopsis  of  these  analysis  is  presented,  with  further  references,  in 


Ref.  [63]. 


lible  B.l 


dE/dx  PID  efficiency  for  data  and  Monte  Carlo 


2.0  - 2.5 
2.5 -3.0 

3.0  - 3.5 


97.8  ±2.9 

09.0  ±3.5 

98.0  ±3.5 


97.6  ± 4.3 

98.7  ±3.9 
96.3  ±3.8 


B.2  EfficienciL-CbeclcB  of  Particle  IdentlBcation  b\'  dEldx  Losses 
To  insure  that  cuts  based  dEldx  particle  ID  information  are  accurately 
modeled  in  Monte  Carlo  a check  on  the  efficiency  of  applying  these  cuts  to 
both  data  and  Monte  Carlo  events  were  compared.  The  particle  ID  based 
cuts  were  imposed  on  a sample  of  high  monieutmn  continuum  D**^s  where  the 
kinematics  of  the  low  decay  D’*  -»  Efix*  can  be  exploited  to  yield  clean 
signals  before  and  after  application  of  the  dE/dx  particle  ID  cuts.  In  Table 
B.  1 the  efficiencies  obtained  by  applying  the  Za  consistency  requirement  used 
in  the  B meson  analysis  for  Monte  Carlo  and  data  are  listed  in  three  different 
CP  momentum  regious.  The  tr  here  implies  a one  standard  deviation  from  the 
measured  specific  energy  loss  for  a particular  particle  hypothesis.  The  results 
for  data  and  Monte  Cario  differ  by  less  than  half  of  a percent  for  D**s  in  the 
momentum  region  of  relevance  to  the  B meson  analysis.  This  indicats  that 
the  effect  of  a 3o  dE/dx  particle  ID  cut  is  well  modeled  by  our  Monte  Carlo 
simulations. 
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B.3  EfBcienfv  Checks  in  D‘*  -»  D°x'*~  Analysis 
Tb  test  tbe  performance  of  the  momentum  and  azimuthal  dependent  im- 
pact parameter  cut  based  on  tbe  analysis  described  in  Appendix  B we  tested 
various  selection  criteria  on  a sample  of  high  momentum  D‘*s.  The  was 
reconstructed  in  the  D'*  -t  D°>rJo^  with  -t  K'n*  . The  performance 
was  tested  by  comparing  tbe  signal  and  background  yields  of  various  cuts  a 
single  track  in  Che  decay  chain  to  the  yield  obtained  with  no  impact  param- 

requirements. 

The  data  used  were  continuum  4s2-4s4  D'*s  (ip  > 0.5)  and  a O’*  contin- 
uum 4s2  Monte  Carlo  sample.  Tbe  energy  losses  of  tracks  used  in  reconstruct- 
ing the  Cfi  candidates  were  required  to  be  within  So  of  the  K/ir  hypothesis- 
In  addition  events  with  Imp.*  - mflg  - nu,+ 1 < 0.0025  were  selected-  These 
selection  requirements  alone  provided  a very  clean  signal  in  the  mp  and  mp. 
invariant  mass  distributions  and  allowed  accurate  yield  measurements. 


B.3.1  Efficiency  Measurements  on  Pions  From  D° 

The  efficiency  measurements  were  made  by  comparing  the  LP  invariant 
mass  yields  with  a cut  on  the  impact  parameter  (DBCD)  of  tbe  ir'*'  from  tbe 
D®  decay  to  the  yield  obtained  without  the  DBCD  cut  on  any  track  in  the  event. 
The  absolute  value  of  the  impact  parameter  was  required  to  within  1,2,3, 4 and 
5 standard  deviations  of  the  measured  spread  in  DBCD,  Tbe  spread  varies 
as  a function  of  momentum,  track  azimuthal  angle  at  tbe  origin  (<p)  and  run 


• Chebychev  with  all  Ga 


Figure  B.l  Fitted  mass  spectra  for  different  DBCD  cuts. 

In  (a)  the  invariant  mass  distribution  without  a DBCD  cut.  In  (b-f)  Che  invari- 
ant mass  distribution  with  a DBCD  cut  of  5, 4, 3, 2 and  respectively. 


Ibble  B.: 


Efficiencies  of  DBCD  cuts  on  the  tt*  in  C®  mass  yields 


rP  mass  distributions  are  shown  in  Figure  B.l  with  the  Gve  different  DBCD 
requirements  (b-f).  The  plot  in  Figure  B.l  (a)  is  the  normalization  plot.  The 
double  Gaussian  Sts  are  superimposed  on  each  histogram. 

The  efficiency  results  for  the  five  different  DBCD  cuts  together  with  a 
constant  6 mm  cut  are  listed  in  Table  B.2.  Also  listed  in  the  table  are  the 
results  obtained  for  the  Monte  Carlo  data,  the  background  efficiencies  and  the 
signal  to  noise  ratio  S/'JS  + B-  The  background  efficiencies  were  determined 
by  integrating  the  polynomial  functions  under  the  mass  peaks.  The  efficiencies 

the  signal  to  noise  ratio  in  the  momentum  regions  above  250  MeV  over  the  Bat 
6mm  cut  commonly  used.  With  a 4<r  cut  the  background  is  reduced  by  11.6% 
while  the  signal  only  drops  by  0.8%  compared  to  tbe  6mm  cut  were  only  5% 


sly  a slight  gain  in  signal  efficiency. 
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Figure  B.2  Fitted  D'*  mass  spectra  for  different  DBCD  cuts 

In  (a)  the  D‘*  invariant  mass  distribution  without  a DBCD  cut.  In  (b-f)  the 
invariant  mass  distribution  with  a DBCD  cut  of  7,5, 3, 2 and  respec- 
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Table  B.3 

Efficiencies  of  tbe  DBCD  cuts  on  the  in  D'*  mass  yields 


DBCD  Cut 

tmc  {%) 

‘r-sncl  m 

«WiS  (^) 

6mm 

96.8 

94.1 

73.2 

78.1 

7ir 

96.5 

95.3 

78.0 

77.8 

60- 

95.7 

94.8 

77.2 

77.7 

5er 

94.8 

93.9 

75.3 

77.5 

4tr 

93.1 

92.4 

73.3 

77.1 

3o 

90.4 

89.5 

69.0 

76.3 

2o 

83.1 

82.2 

61.5 

73.5 

Itr 

59.5 

59.4 

43.2 

62.8 

B.3.2  Effieienev  Measurements  on  Slow  Pions  FVom  D'* 

To  test  the  DBCD  cut  in  a different  momentum  region  the  same  procedure 
was  repeated  on  the  D"*’  invariant  mass  distribution  with  the  DBCD  cut  now 
imposed  exclusively  on  the  bachelor  pion  defined  to  be  below  225  MeV. 

A selection  of  the  invariant  mass  distributions  with  various  DBCD  cuts  are 
shown  in  Figure  B.2.  Since  the  DBCD  spread  for  low  momentum  tracks  is 
large,  an  5 mm  DBCD  cutoff  was  Imposed  In  addition  to  the  momentum  and 
$ dependent  cut.  The  results  for  tbe  various  cuts  on  data  and  Monte  Carlo  are 
listed  in  Table  B.3.  The  reduction  in  efficiency  over  that  expected  for  normal 
distributions  is  due  to  tbe  uon-Caussian  distiibutiou  of  tbe  spread  in  DBCD 
for  low  momentum  tracks  (see  Figure  A.l  (a)).  Imposing  a DBCD  cut  of  7<r 
reduces  tbe  background  I?  22%  while  retaining  97%  of  tbe  signal. 
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